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13  ABSTRACT 

The  SWEM  Allocation  Model  is  a  mathematical  optimization  model  which  allocates 
strategic  weapons  in  a  two-sided  nuclear  exchange  and  provides  damage  results  and 
specific  allocations.  The  model  is  useful  in  defining  levels  oi  parity  between  .he 
nuclear  arsenals  of  two  nations.  Version  II  of  this  model  expands  the  allocation 
matrix  to  include  attack  on  defense  installations  and  provides  the  optimal  tradeoff 
between  attack  of  defense  installations,  weapon  targets  (silos)  and  cities. 

This  document  describes  a  mathematical  optimization  model  which  provides  allo¬ 
cations  and  damage  levels  in  a  two-strike  nuclear  war.  It  is  an  extension  of  an 
allocation  model  (Strategic  Weapons  Exchange  Allocation  Model-Version  I)  which  nas 
been  documented  in  manuals  mentioned  in  the  introduction.  The  new  allocation  model 
allows  attack  on  defense  installations  as  well  as  cities  and  weapon  targets,  and  pro¬ 
vides  optimal  trade-offs  in  allocation  between  defense  installations,  cities,  and 
weapon  targets.  This  document  describes  the  capabilities  and  limitations  of  the 
model,  gives  the  format  of  the  input  and  output  of  the  model,  provides  example  problems 
and  contains  the  annotated  Fortran  listing  of  the  program. 
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FOREWORD 


The  model  described  in  this  document  is  an  extension  of  one  of 
the  Strategic  Weapons  Exchange  Models  (SWEM)  ■which  hare  been  programmed 
to  aid  in  decision  making  concerning  the  procurement,  deployment,  and 
allocation  of  strategic  weapons*  As  detailed  further  in  the  Introduc¬ 
tion,  these  models  were  documented  in  a  recent  series  of  documents  by 
the  Research  Analysis  Corporation.  The  document  provides  the  basic 
information  necessary  to  run  the  model,  but  is  not  intended  to  be  used 
independently  of  the  other  documentation  of  the  SWEM  models. 

This  model,  SWEM  Allocation  Model,  Version  II,  is  an  extension 
of  the  allocation  optimization  model  to  include  defense  installations 
as  targets.  The  model  will  be  an  aid  in  the  decision  making  process 
regarding  the  allocation  optimization  of  offensive  strategic  nuclear 
weapons  against  defense  installations.  It  is  useful  in  defining  levels 
of  parity  between  the  nuclear  arsenals  of  two  nations  and  provides  the 
optimum  tradeoff  between  attack  of  defense  installations,  weapon  targets 
(silos),  and  cities. 


J.  ROSS  HEVERLY 
Vice  President 
Technological  Systems  Group 
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Chapter  1 


INTRODUCTION 


The  model  described  in  this  document  (SWEM  Allocation  Model  - 
Version  II )  is  another  member  of  the  set  of  Strategic  Weapons  Exchange 
Models  (SWEM)  which  have  been  produced  and  documented  by  the  Research 
Analysis  Corporation.  This  model  is  an  extension  of  the  SWEM  Allocation 
Model  which  was  used  in  the  SWEM  Second  Strike  Model  (Version  II )  and 
documented  in  the  manual,  A  Planners  Guide  and  Users /Programmer  Manual 
for  Strategic  Weapons  Exchange  (SWEM)  Second  Strike  Model  -  Version  II, 
RAC-CR-Ug.  The  model  user  should  also  be  familiar  with  the  SWEM  Alloca¬ 
tion  Model  -  Version  I  which  is  documented  :.n  the  manual.  The  Strategic 
Weapons  Exchange  Models  Users  Manual,  RAC-CR-lU,  and  the  general 
description  of  the  SWEM  models  which  is  contained  in  the  document, 
Strategic  Weapons  Exchange  Model  Planners  Guide,  RAC-CR-36.  The  model 
uses  the  same  nonlinear  mathematical  algorithm  for  solution,  the 
Sequential  Unconstrained  Minimization  Technique  (SUMT),  as  other  SWEM 
models.  This  algorithm  is  documented  in  the  manual,  A  SWEM  System/ 
Programmer 1 s  Guide  to  SUMP  -  Version  IV,  RAC-CR-3^*  The  present  manual 
is  meant  to  serve  as  an  annex  to  the  above  documents  and  is  not  self- 
contained. 

This  model  extends  the  capability  of  the  SWEM  Allocation  Model  - 
Version  I  in  several  ways.  The  major  change  is  the  inclusion  of  defense 
installations  as  target  classes.  The  model  now  gives  optimal  allocation 
of  offensive  strategic  weapons  against  city  targets,  weapon  targets 
( I CBM  silos),  and  the  defense  installation  of  each  of  the  city  and  weapon 
targets.  In  addition,  this  new  allocation  model  which  was  included  in 
the  SWEM  second  strike  model  has  several  features  which  were  not  avail¬ 
able  in  the  earlier  version  of  the  SWEM  Allocation  Model.  These  features 
include: 
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.  More  than  a  single  type  of  each  level  of  defensive  weapon:  area 
missile  defender,  terminal  missile  defender,  "bomber  interceptor 
and  SAM,  can  be  input  into  the  model. 

.  Options  for  the  employment  of  defense  against  the  attack  can 
span  the  range  from  a  randomly  committed  defense  to  a 
proportional  defense  which  attacks  incoming  missiles  on  the 
basis  of  damage  capability.  Thus  a  range  of  defense  effectiveness 
can  be  investigated  with  the  model. 

.  First  strike  offensive  weapons  axe  designated  as  countervalue 
only  or  countervalue  and  counterforce.  This  feature  reduces 
the  number  of  variables  in  the  model  and  allows  more 
geographical  areas  or  target  classes  to  be  included. 

.  Bombers  are  allocated  optimally  instead  oi  on  the  basis  of 
population . 

This  manual  documents  the  model  changes  for  the  strategic  analyst, 
user,  and  programmer.  The  first  section  contains  a  description  of  model 
characteristics,  assumptions,  and  uses  for  the  strategic  analyst  interested 
in  applying  the  model  to  strategic  problems.  The  mathematical  formulation 
of  the  model  is  also  provided  in  App.  A.  The  second  section  contains  for 
the  user  a  detailed  description  of  the  model  input,  a  sample  force  table 
with  all  information  needed  to  run  the  model,  and  a  sample  Fortran  code 
sheet  containing  the  model  input  from  the  force  table.  Output  from  the 
model  for  a  series  of  six  sample  problems  using  the  sample  input  data 
is  given  in  the  fourth  section. 

For  the  programmer,  a  Fortran  listing  of  the  program  that  contains 
extensive  comment  cards  is  provided  in  App.  B,  together  with  a  partial 
Fortran  glossary  which  references  the  most  important  Fortran  variables 
with  their  counterpart  in  the  mathematical  notation  in  App.  D. 
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Chapter  2 

MODEL  CHARACTERISTICS,  ASSUMPTIONS,  AND  USES 


Two  types  of  problems  arise  when  considering  attacks  on  defense 
installations.  The  problem  of  the  offensive  planner  is  to  determine  if 
defense  installation  should  be  attacked  and  if  so,  which  installations 
should  be  attacked  and  at  what  level  in  order  to  achieve  his  strategic 
objectives.  The  second  problem  is  for  the  defense  planner  who  must 
deploy  defenses  geographically  and  determine  the  level  of  defenders  at 
each  installation  in  order  to  achieve  his  strategic  objectives.  This 
model  is  designed  to  study  some  of  the  tradeoffs  involved  in  decisions 
of  this  type  at  the  strategic  rather  than  the  tactical  level. 

The  current  scenario  of  the  SWEM  Allocation  Model  does  not  allow 
•for  a  direct  attack  by  the  first  striker  on  the  defensive  installations 
of  the  second  striker.  Such  attack  scenarios  are  feasible  since  an 
attack  on  defense  installations  might  result  in  more  damage  to  value 
targets.  This  additional  damage  is  caused  by  the  remaining  attackers 
which  survive  the  degradated  defenses. 

The  new  model  is  an  optimization  model  which  considers  the 
optimal  allocation  of  offensive  strategic  weapons  in  a  two-strike 
nuclear  war  against  city  targets,  weapon  targets,  and  defense  installa¬ 
tions.  Side  one  is  assumed  to  have  an  inventory  of  different  types  of 
offensive  weapons  which  attack  target  classes  of  side  two  which  can  be 
weapon  classes  such  as  ICEM  silos,  geographical  city  regions,  or  defense 
installations  of  each  of  the  weapon  or  target  classes.  Side  two  is  then 
assumed  to  strike  back  with  all  of  his  remaining  weapons  on  side  one  cities. 
Side  one  weapons  can  be  optimally  allocated  against  side  two  targets  in 
the  model  according  to  a  range  of  criteria.  The  objective  function  of 
the  model  involves  the  difference  in  value  damage  to  the  two  nations 
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and  the  first  striker  objective  can  range  from  maximing  damge  to  side 
two,  to  maximing  the  difference  in  damage,  to  minimizing  his  own 
damage.  Model  output  includes  the  damage  to  both  nations,  the  optimal 
set  of  allocations  again, it  value  targets  and  defense  installations, 
and  many  other  detailed  characteristics  of  the  exchange. 

The  model  takes  into  account  two  main  effects  of  the  attack  on 
defenses. 

1.  The  attrition  of  defenders  due  to  attack  on  their  def  nse 
installation,  theieby  allowing  more  attackers  through  the 
defenses. 

2.  The  missiles  actually  attacking  the  defense  installation 
are  deleted  from  attack  on  other  targets. 

The  new  model  feature  allows  the  user  to  select  up  to  four  side 
one  weapon  types  which  can  attack  specific  types  of  defense  installations 
of  side  two.  Thus  a  side  one  missile  type  can  be  chosen  to  attack  side 
two  SAM  bases  prior  to  the  attack  on  cities,  or  a  side  one  missile  can 
be  specified  to  attack  hard  site  defense  of  ICBMs  prior  to  an  attack. 

The  model  will  determine  if  the  attack  on  defense  will  achieve  the 
strategic  objectives  of  side  one  and  what  the  optimal  level  of  attack 
would  be  to  aichieve  these  objectives.  The  model  will  simultaneously 
optimize  attack  on  all  typ^s  of  defense  installations  for  each  city 
class  or  weapon  target  and  the  attack  on  cities  and  weapons.  The  model 
also  can  handle  several  different  defense  installations  defending  a 
certain  region  as  well  as  multiple  defense  installations  of  the  same 
type  in  the  same  region. 

The  attack  on  defense  installations  is  formulated  on  an  expected 
value  basis.  Damage  curves  as  a  function  of  the  level  of  attack  for 
each  defense  installation  are  calculated  based  on  calculating  the 
probability  that  a  single  missile  will  be  reliable,  penetrate  the 
defense,  and  hit  the  target  for  each  level  of  attack.  The  number  of 
surviving  defenders  free  to  counteract  the  attack  on  values  is  just 
the  total  original  number  at  the  defense  installation  times  the 
probability  that  a  single  missile  will  impact  to  destroy  the  targets. 

Of  course,  the  probability  of  impacting  on  target  increases  the 
greater  the  level  of  attack  on  the  defense  installation.  These  damage 
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curves  are  generated  inside  the  program  and  then  an  analytic  function 
is  fit  in  order  to  provide  the  main  optimization  routine  with  the 
defense  damage  curves. 

The  methodology  for  including  attack  on  defense  installation 
contains  the  following  limitations  and  assumptions. 

I-  Single  type  of  offensive  weapon  can  he  targeted  against  a 
given  type  of  defense  installation.  The  model  allows  either  missiles 
or  bombers  to  attack  any  type  of  defense  installation,  but  only  a  sing!  o 
type  of  weapon  can  attack  a  given  type  of  installation. 

2.  The  first  impacting  missile  on  the  target  defense  installation 
will  destroy  the  defenders  associated  with  that  installation. 

3*  The  geographical  regions  chosen  are  non -overlapping  and  must 
be  the  same  for  each  defense  installation  type. 

4.  The  attack  on  defenses  is  one  sided  in  that  only  the  defenses 
of  the  second  striker  are  subject  to  attack. 

5.  The  solutions  obtained  by  the  model  are  not  guaranteed  to  be 
globally  optimal  because  of  the  non -convex  nature  of  the  problem,  thus 
different  standing  points  may  have  to  be  used  to  obtain  the  globally 
optimal  solution. 

The  first  assumption  is  made  to  reduce  the  number  of  allocation 
variables  and  does  not  represent  a  serious  limitation.  In  general  a 
single  type  of  weapon  with  high  kill  probability  against  a  particular 
defense  installation  will  be  chosen  to  attack.  The  second  assumption 
indicates  that  an  impacting  warhead  will  destroy  either  all  defenders 
or  the  command  and  control  system  necessary  to  target  defenses.  This 
assumption  is  probably  valid  for  all  defenses  except  aircraft  with  AWACS. 
In  this  case  the  aircraft  would  simply  not  be  input  into  the  model  as 
being  destroyable  at  the  defense  installation.  The  third  assumption  is 
the  most  stringent  in  the  current  model.  The  geographical  regions  for 
the  model  must  be  the  same  for  each  defense  type.  That  is,  for  each 
defense  type  the  defense  is  assumed  to  operate  over  the  same  geographical 
region.  This  will  usually  entail  certain  compromises  with  realism  in 
the  choice  of  geographic  regions  and  footprints  of  the  actual  weapon 
system.  Attacks  against  only  one  type  of  defense  installation,  for 
instance,  airbases,  can  be  realistically  modeled.  In  general  the  com¬ 
promise  will  be  defense  conservative  since  usually  defense  regions  will 


be  larger  than  they  actually  are.  The  damage  and  allocations  that  come 
out  of  the  model  will  generally  be  defense  conservative. 

The  fourth  assumption  can  be  overcome  to  a  certain  extent  by  a 
series  of  runs  where  side  one  and  side  two  are  interchanged  to  model  the 
situation  where  both  side  one  and  two  attack  defenses. 

The  fifth  item  is  not  a  serious  limitation  of  the  model.  It  means 
in  certain  instances,  primarily  for  the  cases  where  the  damage  difference 
is  being  optimized,  that  more  than  a  single  run  may  have  to  be  made  to 
insure  a  globally  optimal  solution.  This  is  usually  best  done  by  running 
different  values  of  the  K  parameter,  say  .1,  1.0,  10  instead  of  only  the 
value  of  1. 

The  model  can  be  used  to  address  either  the  problems  of  the 
offensive  planner  planning  an  attack  which  can  include  attack  on  defenses 
or  the  defensive  planner  who  must  deploy  defenses  so  they  are  not 
vulnerable  to  attack  by  the  offense.  Some  of  the  tradeoffs  which  occur 
in  this  type  of  analysis  and  which  can  be  studied  with  the  model  include: 

1.  The  degree  of  defense  dispersal  and  the  number  of  defenders 
at  each  installation,  i.e.,  the  number  of  missile  radars,  airfields 
necessary  to  moke  the  offense  indifferent  to  whether  he  does  or  does  not 
attack  the  defense. 

2.  The  proper  balance  of  different  types  of  defense,  i.e.,  bomber 
and  missile  defense,  to  defend  value  targets. 

3*  The  type  and  proper  level  of  allocation  of  strategic  weapons 
against  defended  target*  between  the  value  target  itself  and  the 
defense  installation. 

In  general  these  types  of  studies  can  be  carried  out  with  the  model 
by  parametrically  varying  the  follovring  model  parameters. 

•  1.  The  area  or  region  of  defense  effectiveness  for  defender  types. 

2.  The  number  of  defense  installations  for  value  targets  and  the 
number  of  defenders  at  different  installations. 

3.  The  overall  level  of  defense  against  different  types  of 
threats . 

One  of  the  limitations  of  the  SWEM  models  generally  is  a  noninteger 
weapon  allocation  to  individual  silos.  Thus  less  than  a  single  missile 
can  be  allocated  to  a  single  silo  and  cause  damage.  The  damage  to  silos 
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from  an  attack  is  accurate  when  survival  probabilities  are  large  enough 
to  cause  at  least  one  missile  allocated  per  silo.  However  with  very 
small  survival  probabilities,  less  than  one  on  one  allocation  will  result 

t 

and  the  damage  will  be  overestimated.  The  current  model  limits  survival 
probabilities  to  be  >  .01  to  partially  overcome  this  problem.  Thus  in 
the  calculation  of  Sij  which  is  done  in  subroutine:  KEADIN,  if  j3ij  comes 
out  to  be  less  than  .01,  it  is  set  to  .01  +  Sij.  This  value  of  .01 
could  be  adjusted  so  as  to  obtain  at  least  one  on  one  allocation  of 

•  •  i 

offensive  missiles  on  silos. 
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Chapter  3 


MODEL  INPUT 

INRJT  DATA  FOR  THE  ALLOCATION  MODEL 

1.  First  card  -  SUMT  parameter  card 
'  2.  Next  set  of  cards  -  Starting  point  cards 

3.  Next  set  of  cards  -  Data  for  the  allocation  model 

4.  Last  three  cards  -  SUMT  option  and  additional  parameters 

Details  concerning  the  groups  of  cards  will  be  discussed  in  order.  The 
reader  may  find  it  convenient  to  examine  the  sample  input  data  contained 
on  the  coding  sheets  in  conjunction  with  the  discussion. 

1.  First  card  -  SUMT  parameter  card 


COLUMNS 

01-12 

12-24 

25-36 

37-48 

U9-60 

61-64 

65-68 

69-72 

FORMAT ; 

E12.0 

E12.0 

E12.0 

E12.0 

E12.0 

14 

14 

14 

NAME  OF  INPUT 

EPSI 

RHOIN 

THETAO 

RATIO 

TMMAX 

MC 

NV 

MZ 

SUGGESTED  VALUE 

.01 

.05 

.00001  2.0 

1800. 

.where  MC  =  the  number  of  inequality  constraints 
NV  =  the  number  of  variables 
MZ  =  0  or  blank 


2.  Next  set  of  cards  -  Starting  point  cards 

This  set  of  data  contains  [NV/6]  +  1  cards.  For  the  sample  problem 
NV  =  87  and  we  have  15  cards .  Hie  starting  point  cards  have  the  format 
6E12.0.  By  setting  NUL  =  1  (column  60)  on  the  first  card  after  the  last 

i 

starting  point  card  STARTB  will  select  a  feasible  starting  point,  thus 
enabling  the  user  to  free  himself  of  this  task.  A  set  of  [NV/6]  +  1 
blank  cards  may  then  be  substituted  for  the  starting  point  data. 

3-1  <®> 


mttas&mtum 


For  those  prefering  to  select  their  own  starting  points,  the 
information  is  read  in  the  following  order:  first  the  non  defense  allo¬ 
cations  are  read  in  by  columns  i.e.,  x(i,l)  (i  =  1,....II),  x(i,JJ) 

(i  =  I  +  l,..r,II),  next  the  values  of  the  first  striker  defense  alloca¬ 
tion  variables  are  read  in  by  rows  u(i,k)  (k  =1,...,  D),  u(ll,k)  (k  =  1, 

. . . ,  D) .  If  D  ^  0  then  there  are  no  allocations  of  side  one  weapons 
against  side  two  defenses. 

3.  Next  set  of  cards  -  Data  for  the  allocation  model 

The  input  data  for  the  allocation  model  is  read  in  once  per  problem. 
For  convenience  of  reference  the  input  data  is  discussed  in  the  order  it 
is  read  in  and  by  the  type  of  data  it  is. 

(3-1)  Problem  Size  and  Option  Inputs 


Format 

1315 


S1 11 

M  J 

JJ 

K  : 

1  D1 

D2 

where 

I 

=  total 

number 

of 

side 

II 

=  total 

number 

of 

side 

M 

=  total 

number 

of 

side 

J 

=  total 

number 

of 

side 

JJ 

=  total 

number 

of 

side 

K 

=  total 

number 

of 

side 

L 

=  total 

number 

of 

side 

N 

=  total 

number 

of 

side 

D 

=  total 

number 

of 

side 

(2nd  card) 


NDEFS  =  side  2  defense  allocation  doctrine 
(-  ^  1  is  proportional 
2  is  uniform 

as  3  the  user  must  read  in  his  own 
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IB1  =  side  one  defense  allocation  ooctrine 
r  -  1  is  proportiona.1 


2  is  uniform 


^  ^  3  the  reader  must  read  in  his  own 

NU1  =  If  NHL  =  1  a  starting  feasible  point  is  calculated  automatically, 
otherwise  the  starting  point  is  read  in 

NU3  =  If  NU3  =  1  a  printout  cf  the  first  strikers  allocations  occurs 
for  every  point  of  the  STMT  algorithm 

NU4  =  The  total  number  of  strategies  (different  ratios  of  k2/kl)  to 
be  considered  for  a  problem  set 

D1  -  The  maximum  number  of  types  of  side  2  defensive  weapon  types 
which  can  be  attacked  by  side  one  weapon  type  1 


0  ^  D1  *  k 

D2  =  The  maximum  number  of  types  of  side  2  defensive  weapon  types 
which  can  be  attacked  by  side  one  weapon  type  2 

0  ^  D2  ^  2  .D1  +  D2  ^  4  and  D1  *  D  2 

(3-2)  Damage  Curve  Fit  Parameters  and  Fractional  Populations 


Format 

10F8.0 


A(1),...,A(II-I) 

A(1),...,A(L) 

P(l),...,P(lI-l) 

P(1),...,P(L) 


each  line  is  one  card 


where 

Ai  =  fitting  constants  for  side  2  populations  (i  =  1, ...,II-l) 

A  =  fitting  constants  for  side  1  populations  ( l  =  1, ...,L) 

Xj 

=  fraction  of  side  2  populations  in  ith  area  (i  =  l,...,II-l) 
P  =  fraction  of  side  1  populations  in  £th  area  ( l  -  1,...,L) 

Xj 

(3-3)  Side  Two  Offensive  Systems  and  Single  Shot  Survival 
Probabilities  Against  Side  One  Defensive  Weapon  Types 


Format 

5E8.0,A8 


n  r  w2  PSI  ZN2  TITLE2 
m  m  m  m  m  m 


10F8.0  /  P  (n  =  1, ...,N) 

L  mn  ’  ’  ' 


First  I  side  2  offensive 
weapon  systems  (weapons 
in  hard  sites) 


n  r  w2  ZI'2  TITLE2 
m  m  m  m  m 


Format 

4e8.0,a6 

10F8.0  /  p  (n  =  1,...,N) 

(_  *mn  ’  ’ 

where 


Next  M-I  side  2  offensive 
weapon  systems 


n^  =  total  inventory  of  side  2  mth  type  offensive  weapon 

r  =  force  reliability  of  side  2  mth  type  offensive  weapon 
m 

w2  =  warhead  yield  (MT)  of  side  2  mth  type  offensive  weapon 
m 

ISI  =  hardness  (P.S.I.)  of  base  (or  side  2  mth  type  offensive 
weapon  (m  =  1, . • . , I) 

ZN2  =  number  of  warheads  delivered  by  side  2  mth  type  offensive 
m 

weapon 

TITLE2  =  name  of  side  2  mth  type  weapon 
m 

p  =  single  shot  survival  probability  of  side  2  mth  type  weapon 
10X1  engaged  by  side  1  nth  type  defender  (n  =  1,...,N) 

(3-*0  Side  One  Offensive  Systems  and  Single  Shot  Survival 
Probabilities  Against  Side  Two  Defensive  Weapon  Types 


Format 

7E8.0,A8 


10F8.C 


m.  r.  e.  wl.  CEP.  FPP.  ZN1.  TITLE1 . 

3  3  3  3  J  3  3  3 


0  [p  (k=l,...,K) 


First  J  weapons 
with  both  CV  and 
CF  capability 


Format  f  m.  r.  e.  wl.  ZN1.  TITLE1. 


5E8.0,A8 

10F8.0 


3 


Pjk  (k  = 


Next  JJ-J 
weapons  with 
CV  capability 
only 


where 

m.  =  total  inventory  of  side  1  jth  type  offensive  weapon 
r  =  force  reliability  of  side  1  jth  type  offensive  weapon 

J 

e.  =  number  of  independently  targetable  warheads  per  each  side  1 
^  attacker 

wl.  =  warhead  yield  per  side  1  weapon  type  j 
J 

CEP.  -  radius  of  circle  (N.  Mi)  of  error  probability  .5  for  side  1 
J  weapon  type  j  (j  =  1,...,  J  only) 

FPP.  =  number  of  cluster  warheads  for  side  1  weapon  typt  j  (j  =  1, ...3 
0 

3-^ 


J  only) 


ZN1.  =  number  of  warheads  per  side  1  weapon  type  j 
0 


TITLE1.  =  name  of  each  side  1  jth  type  weapon 
0 


p.,  =  single  shot  survival  probability  of  side  1  jth  type  attacker 
J  when  engaged  by  side  2  kth  type  defender  (k  =  1,...,K) 


(3-5)  Side  One  Offensive  Weapon  Inventory  Which  Can  Be  Used  For 

Other  Military  Targets  (i.e.,  held  back  from  the  first  strike ) 


Format 

10F8.0 


!o  j  uuj  =  i>  jj) 


where 


uu.  =  number  of  each  side  1  weapon  assigned  to  other  military  targets 


(3-6)  Strategies  to  be  Investigated  for  the  Problem  Set 
(Ratio  of  (k2/kl)  ) 


Format 

10F8.0 


|  RAT(i)  =  ratio  of  k2/kl  for  case  i,  i  =  1, ...,NU4 


(3-7)  Side  One  Reliable  Defensive  Weapons  Inventory  For  Each 
Side  One  City  Class  jj  (i  =  l?L) 


Foimat 

30F8.0 


!o  [  %  (n  =  l,.,.,N) 


where 


d  =  number  of  reliable  side  1  defenders  of  side  1  jfcth  target 
1X1  [a  =  l,...,L)(n  =  1,...,N) 


(3-8)  Side  Tvfo  Defensive  Weapons  Systems  Force  Reliability 


Format 

10F8.0 


rk  = 


where 


.  r£  =  reliability  of  kth  type  side  2  defender  (k  =  1,...,K) 


[3-9)  Side  Two  Defensive  Weapons  Inventory  for  Side  Two  Weapons 
Which  Are  Not  Attacked  by  Side  1 


Note:  If  D  £  0  there  is  not  an  attack  on  side  two  defensive 
weapons  and  all  side  two  defensive  weapons  inventories 
are  read  in  at  this  point. 


For  each  side  two  resource  i  (i  =  1, ...,II; 


(k  =  D  +  1,...,K) 


Format 

10F8.0 


where 


df  =  number  of  side  2  defensive  weapon  type  k  defending  side  2 
ith  type  resource 


FOR  ATTACK  ON  DEFENSES  ONLY 

(3-10)  Side  One  Offensive  Weapon  Survival  Probability  Against 
Side  Two  Base  Defenders 


Format 

10F8.0 


''  (k  = 


where 


S  is  the  survival  probability  of  the  side  one  offensive  weapon 
which  attacks  side  two  defenders  of  type  k 

(3-11)  Detailed  Side  Two  Defense  Structure  (for  defenses  which 
can  be  attacked)  ~~ 

For  each  side  two  defensive  weapon  type  k  (k  =  1, ...,D) 
(3-11-1) 

Format  f  Read  for  each  side  two  resource  type  i  (i  =  l,...,Il) 

(1315)  < 

[  NBASE,  IFIT,  LLO,  LL1 
where 


NBASE  =*  total  number  of  defense  installation  of  resource  type  i 

IFIT  =  1  if  user  wants  to  select  points  LLO  and/or  LL1  over  which 
least  squares  curve  fit  will  be  made 

LLO  =  starting  point  for  curve  fit  (l  unless  user  overrides) 

LL1  =  final  point  for  curve  fit  (the  last  value  for  which  the  expected 
number  of  surviving  side  2  defenders  is  less  than  1  unless 
user  overrides) 

Note:  If  we  chose  to  place  no  kth  type  defenders  of  a  particular 
resource  i  a  blank  card  is  inserted. 

(3-11-2)  Defense  Installation  Information  (where  applicable  NBASE  £  1) 

For  each  defense  installation  associated  with  defensive  weapon 

type  k  and  resource  i 

Read 


Format 

(10F8.0) 


HARDU)  DCAP(je)  ACAP(f) 


l  =  1,...,  NBASE 
1  card  per  each  l 


where 


HARD( l)  =  the  hardness  (P.S.I.)  of  defense  installation  i 


DCAP(  /,)  =  the  number  of  defenders  associated  with  defense  installation  l 
ACAP(jfc)  =  the  number  of  defenders  of  defense  installation  l 


4.  Last  three  cards  -  SUMT  option  and  additional  parameters 

CARD  1*  COLUMN  7  14  21  28  35  42  49  56  63  70  77 

Format 

(1117)  VALUE  3121211  1  321 

(no 

sensitivity)** 


CARD  2 

COLUMN 

01-12  13 

Format 

(6E12.0) 

VALUE 

.01  .00001 

CARD  3 

COLUMN 

7  14 

Format 

(1117) 

VALUE 

l  1 

*  See  RAC-CR-34  ior  discussion  of  SUMT  (version  4)  options. 

**  See  section  VIII  of  RAC-CR-14  for  discussion  of  sensitivity. 
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USER’S  OPTIONS  FOR  THE  ALIOCATION  MODEL 


The  options  available  to  the  user  can  be  classified  as  SUMT 
options  and  model  options.  The  SUMT  options  are  discussed  in  EAC-CR- 34. 
Special  SUMT  option  available  for  use  with  the  model  are  discussed  in 
RAC-CR-14  Section  VIII. 

The  model  options  available  to  the  user  are  as  follows: 

1.  Option  of  attacking  or  not  attacking  side  two  defenses. 

2.  Option  of  selecting  defense  allocation  doctrines  for  both  sides. 

3.  Option  of  selecting  strategies  (objective)  to  be  considered. 

4.  Option  of  having  program  compute  starting  point. 

Explicit  details  on  how  to  initiate  these  options  are  discussed 

in  the  sections  (3-1)  and  (3-6)  of  Input  Data  for  the  Allocation  Model. 

The  choice  of  strategies  is  made  by  selecting  the  ratio  1%/K of 
the  parameters  K^,  Kg  in  the  objective  function  K^|B  -  as  follows: 


Objective 

Ratio  Kg/l^  Function 


0 


0 


1  0-3 

>  107  -P 


Strategy 

First  striker  minimized 
damage  to  himself  (CF 
attack) 

First  striker  minimizes 
difference  in  damage 

First  striker  maximizes 
damage  to  side  two  cities 
(CV  attack) 


These  values  of  ohe  ratio  represent  extreme  points  of  interest 
when  developing  "fish-tail"  type  curves.  Ranges  from  .1  to  100  should 
be  considered  when  sorting  out  local  solutions  particularly  for  the  case 
when  ratio  =  1  since  the  objective  function  is  not  glcjally  strictly  convex. 
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SAMFIiE  INPUT  DATA 

Tables  1,  2,  3,  4,  and  5  contain  the  force  table  data  necessary 
for  input  into  the  model.  Following  these  tables  in  Table  6  are 
Fortran  coding  sheets  containing  the  correct  data  input  for  the  model 
using  the  data  in  the  force  tables. 


FORCE  TABLE  FOR  SIDE  TV/O  OFFER'S E 


1 


o  o  o 
o  o  o 
on  co  on 


00HH00OO0JOO.3- 
H  H  OJ  CM 


o  o  o 

ir\  in  in 

H  H  H 


888888 


COO-d-CMOCMCVICMCVICM 

t^co  i>-no  no  t>-  c — 
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ONONONONONO  O  O  O  O 


ON  ON  ON  ON  on  O  O  OOO 
ONGNONONONOOOOO 
ONONONONONOOOOOOOOOO 
ONONONONONOOOOO 


ON  on  ON  ON  On 
On  on  On  On  On 


O  Q  Q  OQ  On  ON  On  ON  On 

00  00  00  00  do  ON  ON  ON  ON  ON 

C\J  C\J  CM  CM  CM  ON  ON  ON  ON  ON 


s  888 


8QQOQQNONQNONON 
OOOOOnonOnonOn 

IfN  IfN  LfN  IfN  LTN  O'  On  ON  ON  ON 
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8QOOLfNOLTNQQO 
O  m  LfN  J-  CM  3-  IfN 
NO  Jt  H  m  H 


UOUhhOOO 

hiHHwwnncQ 


Jill 

H  H  H  C/3  C/2 
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FORCE  TABLE  FOR  SIDE  TWO  DEFENSE 
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No.  of  Defense  Installations  No.  of  Defense  Installations  No.  of  Defense  Installations 


LTN 

ON 


in 

(U 


FORCE  RELIABILITY  OF  SIDE  TWO  DEFENDER  TYPES 


FORCE  TABLE  FOR  SIDE  ONE  DEFENSE 


PUNCHED  INPUT  TO  THE  MODEL 
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Chapter  4 

MODEL  OUTPUT  AND  SAMPLE  PROBLEM 


The  output  for  the  allocation  optimization  model  can  be  classified 
as  follows: 

(1)  Total  damage  to  each  side 

(2)  Side  one  counterforce  strike  characteristics 

(3)  Side  one  countervalue  strike  characteristics 

(4)  Side  one  counterdefense  strike  characteristics  (when  applicable) 

(5)  Side  two  second  strike  characteristics 

The  output  is  fairly  self-explanatory  and  reference  to  the  output  given  in 
next  section  is  suggested. 

TOTAL  DAMAGE  TO  EACH  SIDE 

The  damage  summed  over  all  population  classes  of  each  side  is 
printed — first  side  one,  then  side  two. 

SIDE  ONE  COUNTEREORCE  STRIKE  CHARACTERISTICS 

The  attack  by  each  side  one  counterforce  missile  type  on  each  side 
two  ICBM  type  is  characterized  in  the  following  fashion: 

(a)  Number  of  each  side  one  missile  allocated  against  each  side 
two  ICBM. 

(b)  Number  of  each  side  one  missile  arriving  over  each  side  two 
ICBM  type. 

(c)  Number  of  warheads  from  each  side  one  missile  type  arriving 
over  each  side  two  ICBM  type. 


(d)  Number  of  warheads  from  each  side  one  missile  type  impacting 
on  each  side  two  ICMB  type. 

(e)  Number  of  warheads  from  each  side  one  missile  type  impacting 
on  each  ICBM  silo. 

(f)  Survival  probability  of  each  side  two  ICBM  attack  by  each 
side  one  missile  type. 

SIDE  ONE  COUiTERVALUE  STRIKE  CHARACTERISTICS 

The  attack  by  each  side  one  countervalue  weapon  type  on  each  side 
two  population  class  is  characterized  as  follows: 

(a)  Number  of  each  side  one  countervalue  weapon  allocated  against 
each  side  two  population  class. 

(b)  Number  of  each  side  one  countervalue  weapon  arriving  over 
each  side  two  population  class. 

(c)  Number  of  each  side  one  countervalue  weapon  impacting  on 
each  side  two  population  class. 

(d)  Number  of  1  MT  equivalents  from  each  side  one  countervalue 
weapon  impacting  on  each  side  two  population  class. 

SIDE  ONE  COUNTERDEFENSE  STRIKE  CHARACTERISTICS 

When  the  option  of  attacking  side  two's  defenses  is  used  the 
following  information  is  provided: 

(a)  Total  number  of  side  one  weapons  allocated  against  each  side 
two  defender  type  by  resource  defended. 

(b)  Total  number  of  side  one  weapons  arriving  over  each  side  two 
defender  type  by  resource  defended. 

(c)  Original  number  of  side  two  defenders  of  each  type  before 
side  one's  attack. 

(d)  Number  of  side  two  defenders  of  each  type  one  by  resource 
defended  surviving  side  one's  attack. 


SIDE  TWO  SECOND  STRIKE  CHARACTERISTICS 

The  second  strike  by  side  two  against  side  one's  population 
classes  is  characterized  as  follows: 

(a)  Number  of  side  two  offensive  weapons  of  each  type  before  i 
side  one's  counterforce  strike. 

(b)  Number  of  side  two  offensive  weapons  of  each  type  surviving 

side  one's  counterforce  strike.  ! 

(c)  Number  of  side  two  offensive  weapons  of  each  type  arriving 
over  side  one's  population  classes. 

(d)  Number  of  side  two  offensive  weapons  of  each  type  impacting 
on  each  of  side  one's  population  classes. 

(e)  Number  of  1  MT  equivalents  delivered  by  each  side  two  weapon 
type  on  each  of  sid^  one's  population  classes. 

SAMPLE  PROBLEMS 

.  •  I  .  ; 

The  following  set  of  problems  is  presented  both  as  illustrative  , 
of  the  model  capability  and  to  .familiarize  the  users  with  the  model  output 

i 

,  I 

The  nuclear  arsenal  and  level  of  defense  for  two  nations  is  given. 
The  problem  is  to  derive  .the  optimal  set  of  allocations  and  the  resultant 
damages  in  a  two-strike  nuclear  war  under  various  scenarios.  In  particu- 

t  \ 

lar,  a  comparison  is  to  be  made  of  the  different  damages  and  allocations 
under  two  sets  of  assumptions: 

a.  side  one  (first  strike)  is  allowed  to  attack  and  destroy  the 
defenses  of  side  two. 

' 

b.  side  one  is  not  allowed  to  attack  and  destroy  the  defenses 

of  side  two.  ! 

Scenario: 

A  two  strike  war  is  assumed  with  side  one  attacking  side  two  with 

a  specified  portion  of  his  input  arsenal,  and  side  two  striking  back 

1 

with  all  surviving  weapons.  The  target  structure  of  side  two  is  divided 
into  six  different  classes  of  targets:  three  city  classes  and  three  ICBM 


i 


weapon  classes.  The  city  classes  are  aggregated  geographically  into  an 
Eastern,  Midwestern,  and  Western  region  with  the  division  roughly  along 
the  85  and  105  axes.  The  three  weapon  targets  are  hardened  ICBM  silos. 
The  defense  of  each  target  is  assumed  independent  of  other  targets  --no 
overlap.  The  level  of  defense  of  target  class  is  given  below. 


Target  class 

Level  of 
missile 
defense 

No.  of  missile 
defense  sites 

Level  of 
bomber 
defense 

No .  of  bomber 
defense  sites 

ICBM3 

624 

2 

0 

0 

ICBM2 

408 

2 

0 

0 

ICBM1 

242 

2 

0 

0 

East 

958 

7 

550 

4 

Midwest 

1035 

13 

600 

4 

West 

370 

3 

425 

4 

Fqr  the  first  set  of  cases  run,  side  one  has  the  option  of  attacking 
the  defensive  sites  of  each  target  class  before  attacking  the  target  in 
each  case,  or  of  directly  attacking  the  targets  without  first  attacking 
the  defenses.  If  the  defenses  are  attacked,  he  must  pay  the  price  of 
the  number  of  missiles  required  to  destroy  the  defenses.  On  the  other 
hand,  if  the  defenses  are  not  destroyed,  part  of  his  force  is  attrited  by 
the  defense.  For  the  second  set  of  cases  run,  no  attack  on  defense  was 

allowed . 

,  1 

For  each  option,  three  different  cases  have  been  run  corresponding 
to  an  attack  by  side  one  maximizing  damage  on  the  cities  of  side  two 
(countepvalue),  an  attack  which  would  minimize  the  damage  to  side  one 
(couhterforce),  and  an  attack  which  minimized  the  difference  in  damage 
between  the  two  nations. 

The  entire  input  arsenal  and  characteristics  are  given  first,  then 
the  curve  fits  for  attack  on  defenses,  and  finally  the  output  for  the 
six  cases  in  the  following  order: 


Attack  on  defenses  option 

Objective 

Yes 

Max  damage  side  two 

Yes 

Min  damage  side  one 

Yes 

Max  difference  in  damage 

1-4 

No 

Max  damage  side  two 

No 

Min  damage  side  one 

No 

Max  difference  in  damage 

l  the  output  to  the  cases,  the  following  numbering  scheme  is  given 

targets  and  defender  types. 

Target  class 

Side  two  resource 

1 

ICBM3 

2 

ICBM2 

3 

ICBM1 

4 

City  class  (East) 

5 

City  class  (Midwest) 

6 

City  class  (West) 

Defender  number 

Defender  type 

1 

.Anti-missile 

2 

Anti-missile 

3 

Anti-bomber 

4 

Anti-bomber 

l*-5 
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attacking  the  defenses  we  have  substantially  decreased  damage  to 
side  one  and  substantially  increased  damage  to  side  two. 


STRATEGIC  WEAPONS  EXCHANGE . MCCEL S  (S HEM) 


CO 

0+ 

UJ 

H 

u 

X 

-J 

X 

►4 

o 

V4 

a 

to 

> 

to 

X 

w 

to 

X 

1-4 

o 

►- 

►4 

a 

a 

UJ 

»- 

X 

X 

Ul 

w 

o 

UJ 

u 

or 

o 

X 

« 

UJ 

»- 

<r 

a 

UJ 

►- 

< 

UJ 

o *- 

«x 

«x 

►- 

a 

o 

a 

k- 

UJ 

<t 

<r 

a 

a 

. 

X 

£1 

Ul 

z 

z 

X 

o 

♦- 

> 

c 

o 

c 

UJ  O 

a 

UJ 

*- 

o 

o 

© 

u 

to  a 

o 

z 

z 

►4 

z 

z 

►4 

»-4 

*•. 

_J  Of 

►4 

►4 

X 

«x  ►- 

C 

> 

u  u 

o 

Z  IO 

o  a 

►4 

<L 

«1 

a 

«X  us  O'  Q.  U  w 

>  a  x  x 

DUCHHU- 


_j  cr  o  o  c  *-< 

CKOUUVI 

<x  w  w  w  VI 


X 

JOO'O’O 

r> 

ni  rj  uj 

44  44 

cr 

• 

•  •  • 

PJ 

O' 

• 

•  .  . 

cr  * 

o 

H 

M  O'  H  (VI 

fO 

to 

O  W  •£> 

44  J- 

*-4 

J- 

r>  O'  r; 

• 

• 

C\J 

cvj  CJ  fy 

•  . 

t-4  H 

44 

“I 

, 

X 

tVJ 

2 

1 

7 

O' 

%D 

n- 

2 

0 

9 

N-  j' 

CJ 

. 

»  .  » 

U> 

. 

... 

n  irv 

o 

CJ 

j  ip  ip 

H 

in 

p- 

UJN  J 

J-  44 

►4 

ro 

CVl  j-  O' 

• 

if 

.T  P-  P- 

•  • 

X  UJ  , 

,  x  z  a 

a 

a 

UJ  z 

3  x  > 

>• 

> 

X 

z  o 

_J  z>  »- 

►- 

o 

O  -J 

<1 

UJ 

O  o  UJ 

Ul 

U 

u 

►4 

o 

wo  -J 

-J 

-J 

r 

vD 

© 

O'  P- 

** 

N. 

►4 

W  M 

►4 

►4 

u 

m 

• 

. 

»  . 

«0 

P- 

CO 

a*  to 

to 

(O 

o 

a 

vTJ 

s.  O 

44 

J’ 

a  uj  to  to  to 
>-  a  *-•»-<  *-<  • 
iiiiu 
>-  o 

u.  u_u-u.tr 

OiOOpH 

o 

i/i  to  to  \r 

U.  IU  u  o  u  c 

«J  «J  «X  «t  «X  > 

H  M  UJ  UJ  UJ  H 


►4  N  |f\  (O  r(  •  • 


oeoentr 
•  •  •  •  r>  u> 

O  O  C  (T  Cl  o 

-r  t»>  r>  •  • 


o  ©  o  o  o  o 

•  •  •  too 

o  o 


Ul 

O 

.  to  to  X 

X  X  > 

’ 

44 

• 

z 

X 

to  to  a 

a  a  a 

• 

o  , 

►- 

HH< 

«x  <  D 

• 

XXX 

X  X  VI 

j 

1 

UJ  ! 

UJ 

1  { 

•• 

1 

o  ; 

o 

t 

a 

» 

r 

* 

►4 

1 

X 

CJ  3COOO  1 

•“*00000 

o  o  o  o 

to  i 

to 

CD 

•  •••03 

•  •  •  too 

#  •  •  • 

JL.  ■  i 

X  epDiOcui  h  c  cr  ^  uj  .»  ir  if  cr 

O.'  ••••rlO  •  •  j  •  •  <T  •••• 

o  »■#  H  r-  o  -4  in  c  n  a  C  c  *J  O'  it  IP  o 

►4  ri  OOK  »  •  K  UJ  vD  M  •  •  vD  in  IT.  PJ 


z  « 

K) 

CSJ 

'  44 

M  UJ 

X 

X 

X 

; 

:  oj 

CD 

1  CJ 

o 

o 

o 

6=7  .CCO  .CCO  .000  .OCi  .041 

Q  42  - i.CCO . -  1.0C0  i.QQC  .995  .639 


c 

o 


•H 

-p 

-P 

to 

CD 

03 

£3 

43 

0) 

U 

W 

EH 

43 

O 

-p 

i — l 

cd 

rH 

43 

• 

o3 

<53 

-P 

-p 

G 

to 

<0 

CD 

O 

03 

43 

-P 

O 

G 

H 

O 

o 

o3 

•H 

0) 

-P 

-P 

« 

-P 

to 

• 

ft 

O 

<0 

CD 

o 

O 

£ 

o 

-P 

o 

-p 

cd 

rH 

<D 

(0 

cd 

i— 1 

43 

CD 

G 

03 

CO 

-P 

o3 

CD 

•H 

iH 

CO 

<D 

O 

to 

CD 

U 

-P 

o3 

O 

a 

O 

-P 

03 

+3 

G 

co 

CD 

o 

o 

CD 

cd 

-P 

bD 

rH 

cO 

OJ 

•H 

O 

CO 

o 

£ 

CO 

O 

• 

C 0 

<0 

to 

H 

+5 

03 

*tH 

rH 

o3 

-p 

S 

<0 

G 

0} 

d) 

cO 

,G 

CD 

CD 

o 

-P 

JG 

JH 

<D 

G 

-P 

CO 

G 

CD 

O 

cd 

N 

03 

to 

43 

•H 

G 

3h 

(D 

-p 

e 

<0 

CD 

03 

•H 

43 

•H 

to 

X 

CD 

S 

to 

0) 

co 

G 

O 

to 

£ 

rH 

P 

O 

G 

<0 

-P 

O 

> 

CD 

O 

<L) 

-P 

43 

(0 

<D 

-P 

on 

03 

to 

-P 

t*- 

O 

•H 

G 

<D 

• 

G 

i— 1 

CO 

CD 

O 

•H 

to 

•rH 

> 

O 

43 

•H 

43 

•H 

H 

-P 

rH 

(D 

O 

rH 

-P 

bO 

CD 

<0 

to 

cO 

*o 

CO 

P 

43 

to 

o 

co 

o 

•H 

o 

o3 

STPATEGTC  WEAPONS  EXCHANGE  HCOEL ?  (SWEH> 


[  ‘ 


u  UJ 

o  o 
c  z 


tr 

►>« 

(X 


or 

«3 

T 

O 


z  o 

X 

c  x 

*-* 

►->  iii 

-I  Or 

«i  *- 

«X  u. 

Z  V) 

rv  o 

UJ 

M 

cr  u/ 

r 

a  o 

*-•  *-i 

«  a. 

►-  */> 

c 

a  >- 

UJ  u 

C  -J 

Z  O' 

c  UJ 

z  z 

►" 

O  <a 
»-< 

fc’5 

►-  o 

»-  c 

«a  *-i 

to  o 

o  e 

C  UJ 

3  UJ 

u  l- 

UJ  z 

-J  vt 

z  c 

«  Ct 

►- 

UJ 

V* 

o 

</} 


CL 

C_' 


•9 

r 

UJ 

c 


X 

cr 

o 


o  o 
o  o 
o  o 


o  o  oo 
r  •  O  O 

o  o 

•  • 
■r* 


OO  oooo 

•  •  •  •  o  o 
o  o 
•  • 


oooo  o  o  ooo  ooo  oo  o  ooo 

•  •••OO  •••  •  ou  ••  •  •  O  <" 

o  o  o  o  o  o 

•  •  •  •  •  • 


OOOCJOC3  OOOOOCJ  DO  OOOO 


o  o 

o  o 
•  • . 


•  •  o  o 

.o  o 

•  • 


«  CD  O 
O  o 

•  • 


ooooo 


(j  O  C  O  (Vj 


z  o 

c  e> 
m  u  i 

c  *- 

HIU> 
ST  C  > 

*-«  u.  a 
a  u.  z> 

tr  a.’  vi 


Hue 

*~<  D  Z 

t)  H  W 
*- 

ft  o 

u  2  a 

>  c  a 

C  X 

o  *-< 
©  z 
z  ►-  — 
>-  z 

>  c  x 

►-  «  13 
ft  a  _» 
axe 
c  e 


cc  OOC  ON 
X  •  •  • 

C  o  o  © 

Qj  J  J  W 


o  OC  ^  N 


C  ^  ' 

«i  z  z  z  : 
cr  x  x  x  3 
«t  z>  r>  3  : 
x-j-j-j- 
o  c  c  c  c 
coy( 
ijj  ww  w  * 
XT 

MUUUJl 

or  ft  o.  ft  c 

c-i  t 

o  ft  U-  u.  I 
z  c  c  c  < 

OI/rlA  VH 

u,  z  z  z  : 

(TCCCi 
a  ft  ft  i 

c  «j  «  «i  < 

X  kJUUL  I 

I-  X  X  x 

UJCOC 

ox*  ? 


0j  o  a  c  o  n 

X  •  •  • 

C  O  o  © 

c.  i  -r  «Ni 


COD  wr 
•  ■  •  *  • 

*f\  IT  IT  O'  o 
N.  S.  CJ  ©  IT 

tr>  K>  cvj 


ODOCW 


pa  o  m  ©  cr 

IT  IT  C  fV  f\j 


|»f  luU'UJ*- 
O  C  C  O  X 
M  *-  *_  *- 

ir  v-  v.  ** 

I1U.I.UU. 

c  c  c  c  c 


X  OORJCN 

cr  •  •  •  •  • 

u  ooNiCi 

H  K.  to  C\J  *■> 


o  c  o  c  o 

:z  Z  z  Z  z 


DOQ  J  J 


O  oooee 

M  OONff'C' 


X  ©  O’©  *  K 

CD  •  •  •  •  • 

O  O  ©  ©  W 

M  O  O  ©  -J  *M 

©  ©  -T 


muWSauww. . 


L 


10 

03 

O 

X 

to 

•H 

g 

U 

03 

to 

0) 

P 

P 

0) 

to 

P 

rH 

03 

>> 

•H 

xf 

-P 

P 

rH 

to 

• 

CO 

03 

03 

CO 

P 

a 

P 

P 

•H 

C3 

0 

•H 

P 

03 

£ 

CO 

< — i 

P 

X 

£ 

ft 

03 

P 

0 

O 

£ 

G 

cO 

cO 

<13 

P 

O 

O 

P 

G 

<0 

0 

03 

P 

O 

03 

P 

o3 

to 

P 

X3 

P 

no 

O 

p 

P 

•H 

xj 

rH 

O 

CO 

G 

G 

•H 

3 

P 

•H 

•rt 

to 

w 

rH 

S* 

03 

03 

rH 

X3 

O 

O 

< 

0 

rH 

-P 

£0 

to 

H 

03 

d 

03 

0) 

P 

P 

• 

no 

H 

O 

X* 

0 

P 

CO 

*H 

03 

O 

£ 

P 

X 

CO 

G 

03 

• 

r« 

0 

XJ 

03 

£0 

>» 

O 

rH 

P 

p 

03 

u 

s 

P 

■rH 

P 

0 

p3 

*0 

H 

g 

-P 

H 

0 

•H 

0 

3h 

H 

P 

P 

•H 

0) 

0) 

O 

CO 

P 

N 

P 

xJ 

G 

ft 

• 

ft 

•H 

G 

g 

CO 

ITs 

O 

£ 

3 

to 

0 

CO 

•H 

0 

•H 

J- 

03 

c 

0 

CM 

03 

• 

P 

•H 

on 

O 

P 

6 

03 

3 

2 

CO 

to 

O 

P 

0 

M 

e 

5 

•H 

0 

P 

H 

M 

U 

<u 

to 

H 

P 

03 

no 

G 

to 

•H 

to 

03 

P 

CO 

03 

■H 

t* 

G 

rH 

G 

•H 

•H 

d 

CO 

to 

03 

P 

03 

p 

0 

xf 

03 

> 

0 

no 

a 

O 

•H 

03 

03 

C3 

P 

+3 

*\ 

03 

3 

to 

O 

-P 

to 

TZ3 

> 

3 

03 

03 

•H 

03 

CO 

•H 

•O 

>> 

p 

G 

fH 

P 

03 

P 

O 

•H 

O 

p 

O 

p 

e 

p 

03 

P 

u 

03 

P 

to 

*rl 

03 

P 

0 

P 

03 

p 

P 

H 

to 

(0 

XJ 

&H 

h- 51+ 


!  S 


m, 

X 

H\«00 

iMoa 

O 

UN  «o  o  o 

CD 

•  •  •  • 

•  •  « 

• 

•  *  •  • 

J 

v» 

N  IP 

K  * 

i 

; 

! 

N  UN 

;■  i 

1  ! 
j 

j  j 

| 

l  .•  i  . 

i  » 

U 

X 

V\  M  OO 

UN  (V  O 

o 

UN  o  O 

CD 

•  •  •  • 

•  •  • 

• 

•  •  •  • 

D 

M 

K  N 

K  M 

H  H 

K  M 
«■»  H 

o 

sc 


v> 

1A 

< 

-I 

o 

► 


o 

X 


o 

x 

CD 

o 


-»  .  M  CD 


* 

X  o 

X 

I 

o 

CD 

■ 

Of  z 

o 

I 

««  c 

n 

(A 

2 

I 

O 

(AO  U 

*-4 

(ACC  Z 

►- 

<  W  M 

Vi 

-1  >- 

►4 

DIN  O 

e* 

or 

(A  «X 

X 

ID 

>  «  (a  a 

CD 

i 

►- 

*-  _j  (a  r 

o 

t 

u 

MO«H 

►4 

< 

U  J 

Sr 

►  u- 

« 

»-  ►-  z 

X 

(a*-  ►  r 

o 

Z  O  »-  3 

1 

H  *-<  _J 

y 

£ 

Ul 

«ttt  OC 

X 

X 

O  U  O 

a. 

?? 

M 

2  - 

o 

ft 

C  C 

»- 

O  Ul 

(A 

WOOD 

1 

>-  z  z  ► 

H 

Ul 

«lHWh 

o>  ►- 

1 

-I 

C  M  o  z 

3 

-jft  ■*  c 

r 

1' 

> 

-J  ft  a  ft 

03 

ft 

<«a< 

o 

•* 

Ul 

m  u 

•4 

►- 

3 

z 

Z  2  — 

3 

r  i  z  x 

C 

3  3  X  C 

o 

-»  J  3  ft 

C  C  -J  U 

k-> 

Ul 

OO  C 

X 

z 

*»  —  o  • 

ft 

c 

^  (A 

o 

I 

Ul  Ul  O 

►4 

I 

Ul 

ft  ft  Ul  Ul 

I 

o 

►  >•  ft 

M 

ia 

►-  ►» 

U  U  U  X 
c  c  c 

v4 

. 

Vi  V  V 

ft 

z  z  z  U 

o 

c  c  c  C 
ft  c  a 

«  «  «t  • 
luluh-a 

M 

X  ■»  x  z 

k I 
£ 


(A 

•  o 

•  M 

Ul 

U  M 

O  ft 
»- 
z  z 

C  Ul 

X 

*-« 

r 

>-  u 

CD 

«  o 

o 

*- 

M 

Ul  *“ 
ft  Ul 
ft  (A 

fix 

V-  o 

z  * 

—  Ul 

J  . 


K«00 
•  •  •  • 

'  c  40 


C  iC  O  □  ClO  oo 

•  •••  •••• 

o«  «  « 


•  •  •  • 

O  K 

•  •  •  • 
OK 

•  •  •  • 
®  K 

i 

1 

Ktdoo 

SiSOO 

K  40  Cl  O 

m  m 

fA  «N| 

n  cni 

IACOO 

UN  40  O  O 

(0  40  O  O 

K  UN 

N  UN 

0.  UN 

o  o  o  o 

o  oo  o 

o  o  o  o 
•  •  •  • 

o  o  oo 

ooo  o 

a  o  o  o 

o  o  oo 
•  •  •  • 

O  O  O  O 
•  ,•  •  • 

O  O  O  O 
•  •  •.  • 

• 

'  ' '  ! 

it 

o  o  oo 
•  •  •  • 

o  o  o  o 

•  ##  •  • 

O  o  o  o 
•  •  •  • 

.  t 

*  t 

O  **  O  O 

o  o  o  o 
•  •  •  • 

•  1 

!  '  ( 

! 

o  o  o  o 
•  •  •  • 

’  1 

o  o  oo 

•  •  •  u 

1 

o  o  o  o 
•  •  •  • 

i 

o  o  o  o 

1 

•‘.j 

1 

,  1 

.  I 

? 

*4 

i 

1 

*r> 

28 

■H  ♦» 


c  § 

o  ND 

p. 

6J  O 

*: 
>*  V 
>-4  XI 
C  O 

.  o  « 


z 

c 

o 


mo  ao  «  *  o  o  s  «  ao 


OlOOO  OUV<4H 


ec 

X 

un  ft  o  e 

UN  ft  O  O 

UN  ft  ■  w4  «4 

m 

•  •  •  • 

•  •  •  • 

•  •  •  • 

>1 

(Nl  40 

ft  40 

(Ni  V0 

(A 

40  * 

40  •* 

40  *  ' 

I-  I 


U-56 


UN 

• 

CD 

O  O  O  O  O 

s: 

* 

o 

oOiO 

- 

CD 

UN  U>  40 

ur 

V  40 

•T 

4-4  UJ  to 

CD 

OOOON 

la  or  NU) 

X 

•  •  •  •  • 

vr  *-  4-  *-• 

o 

o  o  ® 

*-«  </>  i- 

cc 

j-  J’  eg 

or  o'-* 

4-  UJ  o 

40  O  •  UJ 

a  to  z  r 

UJOUiCO 

Oil.  H 

40 

atthui; 

c 

o  o  «c  O  CM 

CUIMOZ 

X 

U.  H  O  H  H 

D 

o  o  « 

co  z  v>*~ 

CD 

<T  <7  (VJ 

Uj  =5  Or  o 

h  OUJ  Z  < 

• 

z  o  >  c  a 

OCX 

c  C  o  *“t 

m  o  z  o  z 

oo 

O  M  Z  4-  ~ 

CD 

O  o  O  4-4  CJ 

*-<  ia  >  4-*  ►-  z 

X 

»-  or  w  >  c  i 

C 

IT  IT  a 

lrtC>H«  3 

CC 

<\J  CM  CT 

►«  u  or  er  a  _j 

H  H 

a  u»  o  a  x  c 

|t«  CU  40  <a  »-*  o 

4-  — 

O  ~  ■ - 

cl  Z  Z  Z  Z  UJ 

r< 

a  x  x  r  x  o 

cc 

oooon 

«DD33>- 

X 

I  JJ  J  Jt- 

c 

o  o  ® 

o  c-  c  c  c 

a. 

-7  -7  CM 

0  0  0  0  2 

v-  w  W  w  C 

V  0. 

4-i  la  UJ  UJ  UJ  <1 

a  a  a  a  a  uj 

»-  >■  v  >■  v  j 

ru 

co  4-  i-  4-  4- 

X 

a  a  CJ  to  -* 

T 

cc 

•  •  *  4  - 

o  u  u  u  u  o 

.J 

IT  IT  U“>  CM  r4 

2  C  C  C  C  «S 

V 

S  N.  CM  JJ  40 

C  UJ 

40  40  PJ 

O  CO  CO  CO  IO 

Ul  2  2  2  2  r 

CO  c  c  o  c  c 

o  a  o  a  or 

c  «S  «  C  «I  u 

v4 

J  Ul  UJ  UJ  U.' 

X 

O  O  O  UN  CT 

4-2222  • 

CL 

CO 

-J 

O  O  40  N.  K 

Ul  c  c  c  c  c 

CO 

Ifi  IT\  (T  *4  H 

C  2  2  2  2  la 

4-4  4-4 

4-4  4-4-4-4- 

so  • 

Ul  UJ  UJ  UJ  K 

cccc  r 

4-4-I4-4- 

rl 

V»  40  CO  10.  »4 

X 

a  «  40  fu  \D 

c 

tu  U  U  U  U 

o 

O  «4  *4 

c  c  c  c  c 

4-4 

40 

•  c  •  •  • 

.  c  c  c  c  c 

2  2  2  2  2 

* 

M 

r 

ojoonk 

CD 

•  •  •  •  • 

o 

O  *4  C'  *"4  *4 

/. 

■  4-4 

O  4-4 

UJ 

J’ 

U  M 

o  or 

4- 

Z  2 

O  UJ 

40 

4-4 

r 

o  o  o  a  \r 

4-  U 

CD 

•  •  •  •  • 

*r  o 

o 

ao  c  40  4 

4- 

4-4 

OOtCCK 

ta  4- 

%0  NO  «» 

or  uj 

O  40 

Cl 

4-  O 

2  « 

m  ia 

i 

), 

_  cr<-7 

<-4 

4-57 


I 


<5*5> 


co 

CD 

CO 

N 

a> 

•H 

CO 

OJ 

s 

<0; 

•H 

O 

C 

II 

•H 

H 

B 

CO 

w 

U 

TO 

CD 

> 

<u 

3 

<D 

CO 

CQ 

CO 

a 

O 

o 

•H 

• 

CO 

P 

CO 

•H 

ft 

G 

XJ 

O 

O 

•H 

o 

P 

co 

CO 

.  . 

c 

C 

O 

o 

r— ! 

<p 

O 

<u 

S 

TO 

TO 

P 

TO 

<D 

P 

CO 

O 

o 


o 

cO 

p 

p 

CO 


1+-58 


_ _ _ .i 

- - 


to  ICBM2  and  ICEM1  which  are  almost  completely  destroyed,  and  the 
countervalue  attack  is  to  the  West  Coast.  The  damages  are  .485 
to  side  one  and  .082  to  side  two. 
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Appendix  A 


MATHEMATICAL  FORMULATION 


MATHEMATICAL  STATEMENT  OF  THE  ALLOCATION  OPTIMIZATION 
MODEL  WITH  ATTACK  ON  DEFENSES  OPTION 

A.  General  Statement  of  Model 

The  mathematical  statement  of  the  allocation  optimization  model  is 


as  follows: 


Minimize  (Kj0  -  Kg 3) 


over  x, J5  ulk 


subject  to 


II  Dl-1 

Z  (x,  j+  Z  u,  J+lt)  s  r1ml  for  j=l  D1  s  D 
i=l  kaO 


2  (xi  j  ^  ^  ^2^2  ^r  j=2  j  •  •  «D-  (Dl-l)  if  Dl<  D  jpjp 

no  constraint  if  D1=D 

II  2 

2  (xj  j  +2  uU+k^  ^  r2m2  for  if  1)2=2 

i=l  k=l  and  Dl=2 


,L  xn  *  rjmj 


for  j=D-kL, . . . ,  J 


2  xn  < 
i=I+l 


for  j=J+l>  •  •  •  >  JJ 


XU>  uik  *  0 


where 


A-l 


ulk  =  number  of  reliable  side  1  weapon  type  (according  to  the  constraint 
conventions  above)  allocated  to  the  kth  type  side  1  defensive 
weapon  of  the  ith  type  resource 

Xj j  =  number  of  reliable  side  1  jth  type  offensive  weapons  allocated 
to  side  2  ith  type  resource 

m5  =  total  number  of  side  1  jth  type  offensive  weapons 
ri  =  force  reliability  of  side  1  offensive  weapon  type  j 
j3  =  fractioned  damage  to  side  2  population 

j3  =  fractioned  damage  to  side  1  population 

K1?K2=  input  parameters  determining  the  strategic  goals 
I  =  total  number  of  side  2  CF  targets  (iCBMs) 

II  =  total  number  of  side  2  targets  ( ICBMs  +  Cities) 

J  =  total  number  of  side  1  offensive  weapon  types  with  CF  and  CV 

capability 

JJ  =  total  number  of  s5.de  1  offensive  weapon  types 
D1  =5  maximum  number  of  side  2  defensive  weapon  types  that  can  be 
attacked  by  side  1  weapon  type  1 

D2  =  maximum  number  of  side  2  defensive  weapon  types  that  can  be 
attacked  by  side  1  weapon  type  2 
D  =  total  number  of  side  2  defensive  weapon  types  that  can  be 
attacked  by  side  1 

Note  if  D=£  there  is  no  attack  of  defenses,  i.e.,  no  ulkvariables. 

B.  Details  of  Model 

First  Strike1  Attack  on  Second  Striker  Defenses 

The  effects  of  the  first  striker* s  attack  on  the  second  striker* s 
defenses  are  modeled  in  the  following  way.  Let 

Llk  =  total  number  of  defense  installations  associated  with  defender 


type  k  of  resource  i 


=  total  number  of  defenders  of  type  k  of  resource  i  at  base  l 


Slki  =  probability  side  1  attacker  arrives  at  the  base  given 
it  encounters  side  2  base  defender 
r'k  =  force  reliability  of  side  2  base  defenders 


The  parameters  qlk  are  then  defined  by  the  Z  that  yields  the  minimum  of 


min 

Z 


2  (E„  -  F,  Z*)S  . 

W=1 


In  order  to  limit  the  number  of  variables, restrictions  have  been 
placed  on  the  number  of  different  types  of  side  1  weapons  that  can  attack 
side  2  defenses.  The  following  five  matrices  depict  the  permissible 
values  of  side  X  weapon  types  j  which  can  attack  side  2  defensive  types  k. 


D=4  D=4 
DlJ+  Dl=3 
D2=0  D2=d 


X 

X 

X 

Dl=2 

D2=l 


Villi. 


X 

X 

X 

_l 

1 — 

X 

D=4 

Dl=l 

D2=a 


V1 2  3  4 

iVl-TT 


X 

X 

X 

X 

D=4 
D1  =2 
D2=2 


The  first  type  side  1  weapon  is  the  only  type  which  can  attack  four 
side  2  defender  types.  This  is  the  case  for  the  sample  problem. 

Thus  the  admissible  index  value  of  j  in  the  equation 

L'k  , 

dtk  =  rk  dik  9tkUlkJ  where  djk  =  Z  hlk  (k=0-,...,D) 

;,=a 

(where  ulkJ  is  the  number  of  reliable  side  1  jth  type  attackers  allocated 
against  defender  type  k  of  the  ith  side  2  resource)  is  limited  to  the  above 
five  situations.  Once  D1  and  D2  (given  D=4)  are  given  the  j  index  value 


is  known  and  ulkj  is  designated  by  u,k.) 


Second  Striker  Defenses 


Let  div  =  number  of  reliable  side  2  kth  type  defenders  (area  or 

terminal)  of  side  2  ith  type  target  (geographic  area  of 
missile  type),  k  is  now  permitted  to  be  greater  than  1 
Pjk  =  single  shot  survival  probability  of  side  1  jth  type  attacker 
when  engaged  by  side  2  kth  type  defender 
Xj j  =  number  of  reliable  jide  1  jth  type  weapons  allocated  to 
side  2  ith  type  target,  (geographical  area,  missile  type) 

Jk  =  allocation  of  side  2  kth  type  defender  to  side  1  jth  type 
attacker  of  side  2  ith  type  target  (geographical  area, 
missile  type) 

JJ  =  total  number  of  side  1  attacker  types, 
jj 

Note  E  v1Jk  =  1. 

Different  defense  doctrines  correspond  to  different  rules  for  choosing 
the  Y1Jk. 

Example  1:  Consider  a  "proportional"  defense,  where  the  defenders 
are  assigned  to  attackers  in  numbers  proportional  to 
the  fraction  of  the  attack  represented  by  a  particular 
attacker  type  and  to  the  capability  of  a  defender  to 
bring  down  a  particular  attacker. 


Then 


Vi  ik 


j  j 

^  *i  j(i-Pjk) 
j=i 


Example  2:  Consider  a  "uniform"  defense,  where  the  defenders  are 
assigned  to  attackers  in  numbers  proportional  to  the 
fraction  of  the  attack  represented  by  a  particular 
attacker  only. 


*Yi  = 

I  xt 

j-i 


Then 


A  more  general  form  to  use  would  be 


xUa3k 

**  =  “ - 

Z  xt  3aJk 


The  aJk  are  constants  to  be  read  in  (or  calculated)  at  the  beginning  of 
the  program. 

The  fraction  of  side  1  jth  type  attackers  of  side  2  ith  type  targets 
surviving  the  defenses  is  then  given  by 


a«  =  (i-Pjk) 
in  example  1, 


in  example  2. 


In  +i\e  general  form  this  reduces  to 


F 


13 


j*dlk 


Jn 


xl3a 


») 
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where  K  =  total  number  of  side  2  defender  types. 

Retaliators  Surviving  the  First  Strike 

The  number  of  side  2  mth  type  weapons  retaliating  on  side  1  is 
y,  3  m  =  1, . . .  ,M, 

where 

rB  =  reliability  of  side  2  mth  type  retaliator, 

nB  =  number  of  side  2  mth  type  weapon  prior  to  side  1  CF  attack, 

SB  =  fraction  of  side  2  mth  type  weapon  surviving  side  1  CF  attack, 


where  yB  =  number  of  side  2  mth  type  weapons  retaliating  on  side  1. 

This  done,  constants  a,n  can  be  read  in,  as  was  done  for  the  side 
2  defenses  such  that 


yx»a»n 

M  _ 

•«1  1 


represents  the  performance  of  side  1  defenses,  whether  uniform  or 
proportional.  The  fraction  of  side  2  mth  type  retaliators  on  side  1 
ith  type  targets  is  then  given  by 


a»n<1i« 


*n  p.x 


^  »a»  ® 

R*  1 


> 


where  N  =  total  number  of  side  1  defender  types. 

In  summary  it  should  be  pointed  Out  that  the  single  shot  survival 
probabilities  (l)  pJk,  (2)  sl3,  and  (3)  P.n  are  instrumental  in  modeling 
the  interaction  of  (l)  the  first  striker  offense  with  the  second  striker 
defense,  (2)  the  first  striker  offense  with  the  second  striker  CV 
targets,  and  (3)  the  second  striker  offense  with  the  first  striker 
defense,  respectively.  Of  course  information  on  the  plk  and  pan 
may  then  be  or  not  be  used  in  determing  the  defense  allocation  doctrine 
constants  ajk  and  aan. 


Damage  Determination 

The  number  of  1MT  equivalents  delivered  on  side  two's  ith  type  CV 
target  is  given  by 

J  j  j 


EQi  Jxi 


+  E  Fj 
j=  J+i 


lX!  JkJ 


i  =1,...,  II-I 


A- 8 


MMfia 


.  , . .  ,  ... 


The  number  of  1MT  equivalents  delivered  on  side  one’s  £th  type  CV  target  is 
EQt  =  E  fguYiA. 

B  =  1 

As  noted  previously  the  fraction  damage  to  side  two*s  ith  type  CV  target  is 
3i  =  1  -  (l  +  Aj/EQj )  exp(-Alv/EQj ) 

The  fraction  damage  to  side  one’s  jgrfch  type  CV  target  is 

Pi  =  1  -  (1  +  exp  Miv^) 

The  total  fraction  population  loss  to  side  two  is 
i  II-1 

^  =  T05.T  ill  P‘Pl 

and  to  side  one 

P  =  — —  Z  PiPi 
TOTP  i=1 

where  TOTP  =  total  population  of  side  two, 

TOTP  =  total  populations  of  side  one, 

Pt  =  population  in  side  two’s  its  area, 

PA  =  population  in  side  one’s  £th  area. 

C.  Remarks  Concerning  the  First  and  Second  Partial 

Derivatives  With  Respect  to  Allocations  Against  the  Defenses 

The  first  and  second  partial  derivatives  of  the  objective  function 

with  respect  to  side  2  defensive  variable  dllt  have  been  evaluated  in  a 

prior  model  (see  RAC-CR-43) .  The  same  evaluation  is  made  in  the  present 


model  even  though  the  dlk  are  not  independent  variable s(  of  the  current  i 

! 

model.  Ihe  following  equations  apply  in  evaluating  derivatives  with  respect 

i 

to  allocations  against  the  defense  ulk  in  terms  of  known  derivatives  with 
respect  to  dlk.  These  equations  are  used  by  the  subroutines  GRAD1  and  MATRIX. 


,  /  .  /  Ui  k 

dtk  =  rk  dlk  qlk  i 

i 

=  r£  d{k  eb*ituik  where  qlk  =  eblk  1 


Mlk 

dGTT  =  blk  dlk 


a2d,  k 

-  =!bik2  dlk 

duIk2 
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SUBROUTINE  READIN 

-  THIS  SUBROUTINE  READS  IN  THE  INPUT  CONSTANTS  OF  THE  PROBLEM  -- 

COHMCN/SHAPE/  X(i%),  DEL(IOC),  A  (  iuC  ,  10  r  )  ,N  V,MC  ,HN  ,NPi  ,NM1 
C0MM0N/0NF/NI,NJ,NII,NJJ,NM,NK,NL,NI1, NJ 1 ,  NI  2  ,  N  J2  ,N  IIJ , IEX,I£D,IEN 
l,IEU,IrXi,TEyEN,Il,I2,I3,I4,I5,I6,I7,NDEFS,N00FFS,NN,NUl,NU2,NU3 
2,IB1,NP  ,NU4,NU5,NU6,NU7,NU8,NC1,N0? 

COMMON/TWC/AA (6)  ,AAB(2)  ,FP0P(6>  ,FP0PB(2)  ,01(9,7)  ,ENL (1G ) , RBCiO > , 

IE  KB  ( in )  ,TITLE2(1C)  ,P(15,7)  ,AB(10 ,15) ,  B(2,  lu  )  ,0E  (2,  It)  ,  EM  (15)  , 

2E  ( 15)  ,EK(15> ,TITLF1(15) ,PB(lu,15) ,AN(15,7) ,SS(5,6),IU(15) 

C OMMCN/ THPF E /PREFER, P RE FIN, PR F INC ,PRED FT  ,  B 1 ( 9 , 5 ) ,R0 (5) 

COMMCN/FOUR/  C(15,7) ,C6(1j,15) ,PLNS (5,6) 

COMMON/FI VE/XX(9,15) ,D( 9, 7) ,EN  <  1C  )  ,  0 (  9 , 5  >  ,  U(  9 , 7) 

C0MMnN/SIX/Z(9,7)  , ZB (2, 1C)  ,F(9,15)  ,S(10 )  ,  FB  <2 , 1C  )  ,  BET  A  (6),Y(6) , 
1BETAP(2) ,YB(2),BFTP,BrT 

COmhON/SFVEN/GFX (9,15,15) ,OFO(9,15, 7) , BE TP <6 ) ,OYX (6  ,15)  ,OBX (6 , 15)  , 
10 YO ( 6 , 7 ) ,OPO(6,7) ,PETPP (2) , EFTPP(6) ,BETBPP<?) 

COMMON/ EIGHT/ CSX (5 ,6) ,D SO (5 ,7)  ,  DF1V( IQ ,5 ,6)  ,0F2X( 10 ,5 , 6) , 

10 FID (10,5,7) ,CF20(1?,5,7)  ,  CYBX  (2 , 5 , 6)  ,0Y E0( 2 , 5, 7 ) 

CC^MON/NIMF/nPey (5,6) ,DRB0(5,7> 

COMMON  /OUTIN/  DnB  ( ?u  )  ,  7N 1  (  2  j  )  ,ZN2(?C)  ,RAT(1D 
COMMON  /ECAL/  H,mi,M7 

DIMENSION  H2(2C)  »PSI(2C)»  Wl(20),  CEP  (2Q ) ,FFP(23 ) , 

CRL ( 2G , 2  3  ) 

DIMENSION  CCAP(2C) ,ACAP<2  0)  ,EJN(?C)  ,QBAR(20)  ,SPRIME(20) ^eARClOO^  ) 
D I  ME  Nc I  ON  PARC (2r )  ,PSUR  (5  ) 

---DATA  FOR  THF  ALLOCATION  MODEL 

---SEE  DESCRIPTION  OF  INPUT  OATA  FOR  DETAILS  - 

WRITE  (6, FEB) 

WRITE  ( 6 , 66n ) 

WRITE(6,25C) 

50  FORMAT  ( / / 45X  ,  37 H* *  IN PU T  CONDITIONS  FOR  THIS  RUN  **  ) 

--(3-n-PROBLFM  SIZE  ANO  OPTION  CAROS  --- 

READ (5,1)  NI ,NII ,NM,NJ,NJJ,NK,NL,NN,NO,NOEFS ,TB1 ,NU1 ,NUT,KU4 
READ (5 , 1 )  NO  1 ,ND? 

.  USER  OPTIONS  ----- 

*****  IF  NU1= 1  PROGRAM  SELECTS  STARTING  FEASIBLE  PCINT  ******** 

*****  jc  fcUTd  A  PRINTOUT  OF  ALL  °OINTS  IS  GIVEN  DURING  THF  SOLUTION 
NU4  IS  THC  TOTAL  NUMBER  OF  STRATEGIES  CONSIOEREC  WHERE  RATIO 
***** (K2/K1)  EQUALS  RAT ( I )  ***** 

. NT  LESS  CR  ECIJAL  ZF^O  NC  ATTACK  ON  DEFENSES - 

IF(NP.GT.w)  WRITE(F,246) 

?46  F0RMAT(/U8X,31HATTACK  ON  OFFENSES  OPTION  USEC  ) 

IF(Nn.Lc.C)  WRIT  p ( 6 ,2  47 ) 

247  FOR MAT(/48X,?9H STANDARD  DEFENSE  OPTION  USED  ) 

WRITE(6,?5)N0EFS ,161,  NU1  ,NU3,MU4 
WRITF (6 , 3 )  NI,  Nil,  NM ,  NK,  NJ,  NL ,  NJJ,  NN 
WRlTF(6,7uC)NC,NCl,N02 
7^  F0RMAT(1HC ,5X  ,515) 

Nil  =  NI  4  1 
NI?  =  Nil  -  NI 


ooo  OOO  OOOO  ro  ooo 


NIIJ  =  Nil  *  NJ 
NJ1  =  NJ  ♦  1 
I F (N J J • EQ • NJ )  NJ 1~N J 
NJ2  =  NJJ  -  NJ 
I  EX  =  NI  *  NJ  ♦  HI 2  *  NJJ 
I EX1  =  IEX  ♦  1 
IF(NO.LT.C)  ND*n 
IED  =  IEX  ♦  Nil  ♦  NO 

11  =  NO 

IF(N02.NE.2.AND.ND*GT.CI  11=  NO  -  <  NOl  -  1  ) 
IF<N02.EQ.2.ANC.ND.GT.O>  11=  2 

12  =  NJ 

13  =  NJJ 
HC  =  13 
NV  =  IEO 
HZ  =  P 

—  (3-2>-0AMAC-E  CURVE  FIT  PARAMETERS  ANO  FRACTIONAL  POPULATIONS 

REA0(5,2>  (AA(I),I=1,NI?) 

WRITE (6,5) 

WRI TE ( 6 , 4)  (AA( I) ,1*1, NI?) 

RE AO (5 i 2 )  (AAB(I> ,I=1,NL) 

WRITE ( 6 ,  E  ) 

WRITE<6,4)  (AAB(I) ,I=1,NL> 

RE AO (5, 2)  (FPOP (I) ,1=1  ,NT2) 

WRITE (6,7) 

WRITE (6, 4)  (FFOP(I) , I=1,NI2) 

READ (5,2)  (FFOPB(I) , 1= 1, NL ) 

WRITE ( 6 , 8 ) 

WRITE(6,4>  (FPOPB(I) ,I=1,NL) 

WRITE  <6,F55) 

WRITE  (6,295) 

5  FORMAT  ( / / 4J X  , 53 H  FORCE  STRUCTURE  ANC  WFAFONS  CHARACTERIS T 

1ICS  ) 

WRITE (P ,910^ ) 

00  lw 9  1  =  1, NI 

-- (3-3 ) -SIDE  TWO  OFFENSIVE  SYSTEMS  AND  SINGLF  SHOT  SUPVIVAL  PR09S • 
- AGAINST  SIOE  ONE  DEFENSIVE  WEAPON  TYPES 

Re AO (5, 91)  £N(I),RB(I) ,W2( I) ,DSI(I),ZN2(I),TITLE2(I) 

RE  AO (5 , 2 )  (  PP(I,J), J=1,NN) 
i:-9  CONTINUE 

00  ill  I=NI1,NM 

-(3-3)  CCNTIKLFC 

RE AO (5, 1C)  EN(I) ,RE(I), W2(I),ZN2(I> , TITLE  2 (I) 

RE AO (5 , 2 )  (  PO(I,J), J=1,NN) 

111  CONTINUE 

00  776  M  =  i , N M 

. COMPUTE  CEFENSE  ALLOCATION  CONSTANTS  FOP  SICE  CNE  - 

IE ( T81-?) 751,733,735 

751  00  752  N= 1 , NN 

752  A  B  ( M , M)  =  1.  -  PP(M,N> 


B-2 
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GC  TO  736 

733  00  734  N=1,NN 

734  AP (M,M)  =  1. 

GO  TO  776 

735  REAj(5,2)  (AB(*,N),K’=1,NN) 

736  CONTINUE 
WRITE  (f>  ,  211) 

211  FORMATdHt ,4. X,17H0FFENSI VE  SYSTEMS) 

WRTTF(6,216) 

216  F0RMAT(1HG,3C'X,  8HSIDE  TWO) 

WRITE(6,2C1) 

23 1  F0RMAT(1HG,4X,76HTYPE  SYSTEM  INVENTORY  FORCE  REL •  WARHEAD 

1YTEL0  NO.  OF  WH  P.S.I.  ) 

20  6  FORMAT  (1H  ,  5X  ,12  ,4X  ,  A 6,  3X  ,F>’ . 2 , 6X  ,F 5. 4  ,9X  , F8 . 3 , 7X  ,F  4  .U  ,  6X  ,F6.2I 
00  2C7  1=1, NI 

20  7  WRITE (6, 2C 6) I  ,TITLE2(I)  ,EN(I) ,RB(I) ,W2 (I  I ,ZN2(I) ,PSI(I) 

00  212  I  =  M1,NN 

212  WRITE (6, 2 J6) I,TITLE2( I)  ,EN(I) , RB ( I) ,W? (I )  , ZN2 < I ) 

00  115  M=1,NM 
00  115  N=i,NN 

115  C8(M,N)  =  AB (M,N )  ♦  A  LOG ( PB (M, N) ) 

WRITE (6 ,221) 

221  FORM AT ( 1HC , 3°  X ,  8HSI0E  ONE) 

WRITE (6 ,226) 

226  FORMAT (  /5X,lu8HTYFE  SYSTEM  INVENTORY  FCRCE  REL.  WARHEAO  Y 

1IEL9  NO.  OF  WH  NO.  INO.  TAR.  WH  NO.  CLUSTER  WH  CEP) 

231  FORMAT (1H  , 5X ,12 ,4X , A 6, 3X ,F 8 . 2 , 6X ,F5 . 4 ,9X , F8 . 3 , 7 X , F4 .0 , 11 X ,F4 . u ,1 4 
1X,F4.C,9X,F5.2) 

00  131  1=1, NJ 

-- (3-4 ) -SIDE  CNE  OFFENSIVE  SYSTEMS  AND  SINGE  SHOT  SURVIVAL  FRORS. 
. AGAINST  SICE  TWO  DEFENSIVE  WEAPON  TYPES 

READ  (5,94)  EM  (I)  ,P99(I  ) ,  Ed)  ,W1  Cl)  ,CEP<  I)  ,FPP  (I)  ,2N1  (I)  ,  TITLEK  I) 

RE  AO  (5 , 2  )  (P(I,J)  ,J=1,NK) 

WRITE  (6, 231)  I, TITLEK  I)  ,EM(  I)  ,  RBBCI )  ,W1(  I )  ,ZM  ( I)  ,E  (I)  ,  FPP(I)  ,  CEP  ( 
II) 

131  CONTINUE 
00  132  I=NJ1,NJJ 

-(3-4)  CONTINUED 

READ  (5,91)  EMI)  ,REB(  I)  ,E(I),W1(I),ZN1(I)  ,  TI TLE1  (I) 

RE AD (5 , 2 )  (  P (I , J ) , J=1 ,NK) 

WRITE  (6, 2  31)  I, TITLEK  I)  ,EM(i)  ,  RRF  (I )  ,Wi(I)  ,ZN1(I)  ,c  (I) 

132  CONTINUE 
00  136  1=1, NJJ 

- COMPUTE  O’FENSF  ALLOCATION  CONSTANTS  FOR  SICE  TWO  - 

IF  (NOEFS  -  2)  151,133,135 

151  00  152  J=1,NK 
AN ( I , J )  =  1.  -  P(I,J) 

152  CONTINUE 
GO  TO  136 

133  DC  134  J= 1 , NK 
AN ( I , J )  =  1. 

134  CONTINUE 


o  o  o 
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60  TP  136 

REAn<5,2>  ( AN (I ,J) ,J=1 ,NK) 

WRITF(6,18> 

W°ITF ( 6,  4)  (AN( I, J) , J=1,NK) 

1?6  CCNTTMUC 

00  2X4  1  =  1, NI 
00  2-4  J= 1 , N J 

RL(I, J)=2.8*(kl( Jl  **? .3  33)* (PSI(I)-7.3  7)**(-w.35?) 
SS(I,J)=0.5**((FPP(J) >* (RL ( I , J ) /CFP ( J) )**2) 

C 

- CGMFUTE  SINGLE  SHOT  SURVIVAL  crops.  FOR  SIDE  TWC  HAROSITES  VS. 

. -SIDE  ONE  OFFENSE . — 

***  IF  SINGLE  SHOT  SURV.  PROB.  IS  LESS  THAN  .01  INCREASF  IT  BY  .01  ** 

IF(SS(I,J) .LT.l.CE-2)  SS(I,J)  =  SS(I,J)  ♦  .ul 

2*' 4  CONTINUc 

395  FORMAT  ( / /2l X , 7* HS ING Lr  SHOT  SURVIVAL  PROBABILITIES  (SIDE  TWO  ICBM 
IS  VS  SIDE  ONE  MISSILES)  * ** /6X , 14HS IDE  TWO  ICBMc ,26X,26HSIQE  ONE  M 
2ISSILFS  (ACROSS ) /9X ,6H ( DOWN ) ) 

400  FORMAT  ( 7X , A6 ,12 (4X ,F6. 4) ) 

525  FORMAT  (  /15X ,12 (4X , A6)  ) 

WRITE  (6,655) 

W  R I T  E  ( 6 , 9  C  w  0  ) 

WRITE  ( 6 , 3Q5 ) 

WRITE  (6,525)  (TITLE1 (J),J=1,NJ) 

DO  2l 5  1=1, NI 

205  WRITE  (6,4-j)  IT  TL  E  2  ( I)  ,  (  SS  (I ,  J  )  ,  J=  1 ,  N  J) 

795  FORMAT  (/ /21 X , 76 HS INGLE  SHOT  SURVIVAL  PROBABILITIES  (SIDE  TWO  OFFS 
1.  VS  SITE  CNF  HEFORS.  )  * ** /SX , 14HS TOE  TWO  0 F F S .  ,  26 X  ,  26FS IOE  ONE  0 
2EFDRS  •  (ACROSS) /9X,6H( DOWN)) 

WRITE (6 , 795) 

WRITE(6,625) (N,N=1,NN) 

DO  5E  5  M  =  1 ,NM 

8;- 5  WRITt(6,4Cw)TITLE2(M)  ,(FB(M,N)  ,M-1,NN) 

DO  141  J=1,NJJ 

I F  (Wi ( J) .LE.i.u)  EK(J)=ZN1( J)*W1( J) **2  .6  67 
IF(WKJ).GT.l.C)  FK(J)=7N1( J)*M1( J) **.5C 
141  CONTINUE 

DO  142  J=1  ,NM 

IF(W2(J)  .GT.l.L)  EKB( J) =7N2(J) *W2C J)**C. 53 
IF  (W2(J)  .LE.l.C)  Exn( j) =  7N2 (J) *N2(J>** j.*67 
14?  CONTINUE 
WRITE(6,49C) 

495  FORMAT  ( //21X  ,76HSINGLE  SHOT  SURVIVAL  FROPAPILIT IES  (SICE  ONE  ICBM 
IS  VS  SI Dc  TWO  OEFDRS.  )  * **/5X , 14HS IDE  ONE  I CBMS , 26X ,26HS IDF  TWO  9 
2EFPRS.  (ACRCSS) /9X ,6H( OOWN) ) 

625  FORMAT  (  /12X ,12 (4X , I  6) ) 

WRITC (6 , 625) (K,K=1,NK) 
no  6^6  J=1 ,N J J 

6:6  WRITE ( 6 , 4  C  j )  TITLE1 (J ) , (P  (J ,K) , K=l, NK) 

00  155  I=NI1 , NM 
155  S ( I >  =  R T ( I ) 

00  16.  1=1, NI 
DO  16*  J  =  1 , N  J 

FLNS(I,J)  =  E ( J)  *  AL  OG (SS ( I  ,  J ) ) 

16n  CCNTINIJF 

DO  165  J=1 , N J J 


ooo  »-*•  oooo  ooo  ooo 


165  C(J,K)  =  AN(J,K)  ♦  AL0G(P(J,K>> 

WRITE(6,9lu()) 

--(3-5)-  SIDE  ONE  OFFENSIVE  INVENTORY  USEO  FOR  OTHER  TARGETS 

RFA0(5,2)  JUKI), 1=1, NJJ) 

WRITE ( 6 , 20 ) 

WRITE(6,4)  (UU(I)  ,1=1, NJJ) 

00  143  J=1,NJJ 

IF  t  UU(J) .LT.EM(  J)  )  EM ( J )  =  EM(J)  -  UUCJ) 

143  CONTINUE 

00  147  J  =  1 , NJJ 
147  EM<j)=R9B(J)»EM( J) 

WRIT  E (6 , 9C  c  n ) 

NU8  =  1 
WRITF (6 , 797) 

--(3-6>-STPATEGIES  TO  EE  INVESTIGATED  (RATIO  OF  K2/K1) 

RE AQ ( 5 , 9 )  (RAT  (I)  ,1  =  1,10 
00  799  1=1, NU^ 

799  WRITE(6,798) I ,RAT(I) 

797  FORMAT  (  // 10  x* , ?8  HRA TI  OS  SELECTED  FOR  STRATEGY/1FX,21HCASE 
1E(K2/K1>  ) 

798  F ORMAT  (  16X , 12 ,5X , £12 .4  ) 

I F (  R AT ( 1 ) . GT • 1 • E+7 ) PRFBET  =  loO. 

I F (  RAT (1) #GT.1.E+7)PPEFER=C . 

I F (  RAT(l) .GT.l.E+7)  GO  TO  937 
PREFER  =  lwu. 

PRESET  =  PREFER  *  RAT (1 ) 

537  CONTINUE 

WRITE  (6,655) 

WRITF ( 6 , 9 C l u  ) 

241  FORMAT (1HC ,4jX,17HCEFENSIVE  SYSTEMS) 

WRITF (0, 241 ) 

WRITF (6 , 221) 

WRITE(6,15> 

OC  13 1  1=1, NL 

-- (3- / ) -SIP 6  CNE  RELIABLE  DEFENSIVE  WEAPONS  SYSTEMS  INVENTORY  FOR 
. POPULATION  CLASS  L . 

RE  AO (5 , 2)  ( OP ( I , J) ,J=1 , NN) 

WRITE  (6,t*)  (  C  E  ( I  ,J)  ,J=1  ,NN) 

CONTINUE 
WRITE(6,216) 

WPITE(6,24) 

—  (3-3) -SIDE  TWO  OEFFNSIVE  SYSTEMS  FORCE  RELIABILITY  — 

R  r  A  9 ( 5 , 2 )  (ROCK)  ,K  =  1,NK) 

WRITE(e,4)  ( R  L  i  K)  ,  K  =  1 , NK ) 

IF(NO.FC.NK)  GO  TO  25^ 

IF(NP.GT. O.ANC.NO.LT.NK)  N03  =  NO  ♦  1 
IF(NO.LE.G)  NO 3  =  1 
90  262  1  =  1, Nil 

C  -- (3-9 ) -SI DE  TWC  CEF^NSIVE  WEAPONS  INVENTORY  FCR  SICE  TWC  WEAPONS 


B-5  (jiAch 


VALU 


-  -  ■"  - . —  —  --Tit-  rami-i 


ooo  oooooo  ooo  oooo 


c 

c 


NOT  ATTACKFC  BY  SIOE  CNE 


263  REA0<5,2) <0L (I,K > , K=ND3 ,NK> 

IF(NO.LE.Q)  GC  TO  1*17 

— (3-1CI-SI0E  ONE  OFFENSIVE  WEAPON  SURVIVAL  PROBABILITY  VS.  SIOE  TWO 
. DEFENSE  INSTALLATION  DEFENDERS - 

253  RE AD (5 , 2 )  (PSUR ( K) ,K= 1, NO ) 

WRITE (6,971) 

971  F0PMAT(/13X,49HSURVIVAL  FROB.  CF  SIOE  ONE  OFFS.  VS.  BASE  CEFORS.) 
WRITE (6, 4) CPSUP(K) ,K=l,NO> 

COMPUTATION  CF  MAXIMUM  DAMAGE  CURVE  *********** *** 

Cl  =  (  SORT  ( 5  •)  -  1.)  /  2. 

C2  =  1.  -  Cl 
DO  5LC  K  =  1 ,  NO 
DO  5uC  1  =  1, Nil 
LC  =1 

WRITE (6 ,80  3) I  ,K 

—  (3-11-D-SICE  TWO  RESOURCF  I  NO  OF  DEFENSE  INSTALLATIONS 
REA0<5,1)  NBASE,IFIT,LLC,LL1 

.  USER  FITS  CURVE  OVER  INTERVAL  (LLO ,LL1 >  IF  IFIT  =  1 

IF(NBASE.LF.O)  Q(I,K)  =  1. 

IFCNPASE.LE.QJGL  (I,K)  =  .Cl 
IF (NBASE  .LE.  * )  GO  TO  5C r 
NIK  =  NBASE 

WRITE(6,256)  K,  I,  NIK 

256  F0RMAT(1HC ,4X,15HQEF.  WEAP.  TYPE,I3,31H  RpG  ION  CR  CITY  CLASS  OEFE 
1N0ED, 17 , 38H  NO.  CF  BASES  DEFENDED  BY  WEAP.  TYPE, 13) 

IF(IFIT.EC.l)  L0 -LLO 

IF(TFIT.EQ.l)  NF IN  -  LL 1 

IF ( IFIT .EC  .1) WRITE (6, 97  7) LL  }  ,L L 1 

977  FORMAT (/1CX,47HUSER  CHOOSES  TO  FIT  CURVE  USING  STARTING  FOINT  ,17 
1  ,15H  AND  END  POINT  ,I4> 

SUM  =  C. 

WRITE (6,261) 

261  F0RMAT(1HC,9X,116HBASE  NO.  H ARONESS ( PS  I )  NO.  OF  CEFORS.  OF  BAS 
IE  NO.  OEFCRS.  AT  BASE  OFFENDED  BASE  UNOEFFNDE C  BASE/54X,31H 
2SURVIVAL  FROB.  SURVIVAL  PROB.) 

DO  21  wt  J  =  1  ,N IK 

--(3-ll~2)-SIOE  TWO  DEFENSE  INSTALLATION  INFORMATION 
READ (5, 2)  HACC(J),CCAF( J) ,ACAP(J) 

TERM  =  2.0  ♦  <W1(K)**.333)  *  (  HARD ( J ) -7 . 77)  **  (-0.352) 

SPRIME(J)  =  .5  **  ( (F  FP  (K ) )  ♦  (  TERM  /  CEP(K))  *«2> 

OPAR(J)  -  Cl.  -  SPRIME(J))  *  (  RD(K)  *  (PSUR(K)  -1.)  +  1.) 
TF(ACAP(J)  .LT.l.)  CB A R ( J)  =  1.  -  SPRIME(J) 

EJN(J)  =  DC AP  ( J) 

V  AL=1-CB  AR ( J ) 

WRITE (6,266)  J,H6P0(J), ACAP(J)  ,OCAP(J)  ,V A L  ,SRRIM E ( J > 

266  FORMAT (1H  , 1 5X ,  I  2 , 1CX  ,F 6. 2 , 12X , F5 .1 ,16 X, F 6 .1 , 14X , F5 .4 , 1 2X  ,F5. 4 ) 


o  o  o 


210  SUM  -  SUM  ♦  EJN(J) 

EBARP  =  SUM 
OL(I,K>  *  SUM 
00  23l  N  =  l,l  Vj 
BIG  = 

00  22C  J=1,NIK 

TERM  =  EJN(J)  *  OBAR(J) 

IF(TERM.GT.eiG)  JMAX  -  J 
IF(TERM.GT.BIG)  PIG  =  TERM 
220  CONTINUE 

FBAR(N)  =  SUM  -  EIG 
EJMJMAX)  =  E JN ( JMAX)  -  BIG 
SUM  r  EBAR (N) 

ACA°  (JMAX)  =  ACAP(JMAX)  -  1. 

IF(ACAP(JMAX) .EQ.O.)QBAR( JMAX)  si.-  SPRI*F(JMAX) 

IFdFIT.EG.il  GO  TO  23C 
IF (SUM. LT • 1. I  NF IN  =  N 
If (SUM.LT.l.)  GO  TC  24C 
230  CONTINUE 

NFIN  =  lGjr, 

IF(IFIT.EG.l)  NF IN=LL  1 

LCAST  SQUARES  FIT  TO  MAX.  UM .  CURVE  USICG  FIBONACCI  SEARCH  ♦♦♦♦♦ 

240  TU  =  1. 

TL  s  r. 

TB  s  Cl 
TA  =  C2 
I  SET  =  1 
V ALB  s*l. 

00  3u5  N=LC,NFIN 

305  VALB  =  VALB  ♦  (  EBARO  ♦  TB**N  -  EBAR(N))**2 
00  411  L=l,6" 

I F  (ISET -2)  34G, 310,310 
31  j  VALB  =  *. 

00  33C  N=LO , NFIN 

33  0  VALO  =  VALB  ♦  (  EBARC  *  T  B*¥N  -  EBAR(N))**2 
IF (IS FT -2) 34"  ,36C , 34C 
340  VALA  -  0. 

00  35C  N=LO , NFIN 

35^  VALA  =  V ALA  ♦  (  EBAR?  *  TA*»N  -  EBAR'N>)**2 
3F0  I F (VALA  -  VALB)  37 0,391,380 
37*  TU  =  TB 
TB  =  TA 

TA  =  TL  ♦  C2  *  (TU  -  TL ) 

VALO  =  VALA 
ISET  =  1 
GO  TO  410 
38*1  TL  =  TA 
TA  =  TB 

TB  =  TL  ♦  Cl  *  (TU  -  TL) 

VALA  =  VALB 
ISET  s  2 
GO  v0  41" 

390  TL  =  TA 
TU  =  TB 

IF  ( ABS (TL  /  TU  -  1.) .LE. l.E-3)  GO  TO  45" 

TA  =  TL  ♦  C2  *  (TU  -  TU 


N  O'  vn 


TB  =  TL  4  Cl  *  (TU  -  TL> 

ISET  =  3 
41;  CONTINUE 
42  ?  F0RMAT(1HC,6E16.  8) 

45:  Q  Cl  =  (TL  4  TU)/  2. 

WRITER, 271) 

271  F0ON1AT(1HC  , 38X, 5 1H EXP  EC TEO  NO.  OF  SURVIVING  OEFFNDERS  VS.  ATTACK  S 

1IZE/38X,4SHNC.  ATTACKERS  EXP.  NO,  SURV.  CURVE  FIT  VALUE) 

00  4 7b  N=  1  * NFIN 
VAL  =  EBARt  *  C(I,K)**N 
470  WRITF(6,276)N,FBAP(N) ,VAL 
276  FORMAT ( 1H  , 4 3 X , 1 3 , 9X , F9 . 2 , 9X , F 9 . ? ) 

5u^  CONTINUE 
1C  7  WRITF<6,9) 

WRITE (6  »  625) (K,K=1,NK) 

00  8C6  1=1, Nil 

8C6  WRITE (6,825) I, COL <I,K), K-1,NK) 

OP  ICC  1=1, Nil 

DC  lit  K  = 1 , N K 

OL ( I , K )  =  PC ( K)  *  0L( I, K) 

100  0(1, K)  =  0 L ( I » K ) 

IF (NO.Lc.w)  GC  to  lJOn 
no  105  1  =  1, Nil 
DO  1C  5  K  =  1 , NO 
lu  5  B1 ( I ,K)  =  ALOG (Q (I , K)  ) 

WRITE(6,193) 

193  F0RMATC1HC,  4  1HCURVE  FIT  VALUES  FOR  Q(I,K),  K(ACCROSS)  ,  I(OOWN)) 
WRITE (6 ,625) (K,K=1,ND) 

825  F0RMAT(  7X  ,  16 , 1? (4X,F^, 1) ) 

925  FORMAT(  7X  ,  16 ,12  (4X ,F6. 4) ) 

00  841  1  =  1, Nil 

841  WRITE(6, 925)1, (  G ( I ,K ) , K= 1 , NO) 

WRITE(6,9cl?) 

IOCw  IF  (NU1.EC.1)  CALL  STARTS 
RF  TURN 

1  F QRMAT (1415) 

2  F  CRM AT ( 1 0  F8  •  n ) 

3  FORMAT (IHl ,4X,9HSI0E  2  --,I7,11H  CF  T  ARGETS  ,1C X  ,17  ,8H  TARGETS, 

1  9X , I  7 ,13H  WEAPON  TYPES, 17, 15H  OE  F  END  FR  TYPES  // 

2  5 X , 9HS 10 F  1  --,I7,21H  CF  4  C V  WEAFON  TYPES, 17, 

3  17F  GEOGR  AFFr  C  AREAS, 17, 13H  WEAPON  TYPES, 17, 

4  15H  0  EEENOF  R  TYPES) 

4  FCRmaT(  8X , 1 jF12. 4) 

FCRMAT(1H?,4X,39HFITTING  CONSTANTS  for  SIDE  2  POPULATION) 

FORMAT (lHu ,4X,39HFITTING  CONSTANTS  FOR.SIOp  1  POPULATION) 

FORMAT (IHl , 4X ,4 1HFPA CT ION  OF  SIOE  2  POPULATION  IN  ITH  AREA) 

8  FCRNAT(1H:,4X,41HFRACTI0N  OF  SIOE  1  POPULATION  IN  LTH  AREA) 

9  FORMAT (1H- ,4X,47HNUMBER  CF  K  TH  TYPE  SIOE  2  OEFDRS •  (ACCROSS)  OF, 

123HI  TH  TYPE  TARGFT  (OOWN)) 

1j  F  CRM AT ( 4F 8  •  p  ,  AE  ) 

11  FO°MAT ( 1HG , 4X  ,11HSI0E  2  TYPE,I2,22H  WEAPON  --  NUMBER  ,F7.2,14H 

1  RELIABILITY, F7. 4, 1SH  WARHEAn  YIELD,  F7.2,  18H 

2  NO.  OF  WARHEADS, F7. 2/  5X , 7HN  AME  -  ,A6) 

12  FORMAT (lHr ,4X,44HSINGL E  SHOT  SUPVTVAL  PROBABILITY  ENGAGEO  BY  , 

1  24HSI0E  1  NTH  TYPE  DEFENDER) 

13  FORMAT (lHu ,4X,35HSI0E  1  DFFFNSE  ALLOCATION  CONSTANTS) 

14  FORMAT (lHu,4X,37HSI0E  2  OFFFNSIVE  ALLOCATION  CONSTANTS) 

15  FORMAT (1HI ,4X,44HN UMBER  OF  RELIABLE  SIOE  1  CEFENOEPS  OF  SI9F  , 


w*m 


1  12H1  LTH  TARGET) 

16  FORMATdHC  ,4X,liHSI0E  1  TYPE,I2,24H  WEAPON  --NO.  OF  WEAPONS ,F7.2, 
114H  RELIABILITY, F7. 4, 34H  NO.  OF  INO.  TARGETAPLE  WARHE AOS, F7 .2 , 
216H  WARHEAD  YI ELO ,F 7. 2/5X ,7HNAME  -  ,A6,18H  NC.  CF  W ARHEAQS ,F7  . 
32) 

17  F0RMAT(1H:,4X,48HSINGLE  SHOT  SURVIVAL  PROB.  WHEN  FNGAGFO  BY  SIDE 
1  , 19H2  KTH  TYPE  DEFENDER) 

16  FORMAT dHC,4X,35HSIDE  2  OEFFNSE  ALLOCATION  CONSTANTS) 

19  FCRMAT (1H0,4X,46HSINGLE  SHOT  SURVIVAL  PROB.  OF  SIDE  2  MTH  TYPE 

1  4  7HRET AL IATOR  ATTACKED  BY  SIDE  1  JTH  TYPE  ATTACKER) 

20  FCRMAT  (1HC , 4X,47HNUMBER  OF  EACH  SIDE  1  WEAPCN  ASSIGNED  TO  OTHER  , 

1  16HNILIT ARY  TARGETS) 

21  FCRMAT (1HG,4X,6HPREFER,F7.2  ,7H  PREFIN,F7.2  ,7H  FREINC,F7.2  , 

1  4H  PHI  ,F  7. 4  ,7H  PREBET, F7. 2  , ’’H  TL,F7.2  ,4H  EPS, 

2F7.4) 

23  F0RMAT(1HC,4X,32H  NO.  SIDE  2  OFFS.  TYPE  I  WEAPONS) 

24  FCRMAT(1HG,4X,4CHPFLIABILITY  OF  KTH  TYPE  SIOE  2  DEFENDER) 

25  FORMAT(1HO,4X,5HNDEFS,I2,9H  IB1=,I2,6H  NU1=,I2,6H  Nll3=,I2,6H 

1  NU4=,I2) 

655  FORMAT  (1H1) 

660  FORMAT  ( 29X , 62H* *********  STRATEGIC  WEAPONS  EXCHANGE  MODELS  (SWEM) 
±* ********** 9  / ALLOCATION  OPTIMIZATION  MODEL  ) 

91  FORMAT ( 5E  8  •  0  ,A8) 

94  FORMAT (7E8.0  ,A8) 

92  FORMATdHC  ,4X,11HSI0E  2  TYPF|I2f22H  WEAPON  --  NUMBER  ,F7.2,14H 

1  RELIABILITY, F7. 4, 16H  WARHEAD  Y IELO, F7 •? , 9H  p.S.I  ., F7.2, 18H 

2  NO.  OF  WARHEADS, F7. 2/  5X , 7HN AMF  -  ,A6) 

93  FORMATdHC,  4X,11HSI0F  1  TYPE,I2,24H  WEAPON  --NO.  OF  WEA  F0NS,F7 .2, 

114H  RELIABILITY, F7. 4, 34H  NO.  OF  INO.  TARGETAELE  WARHE ADS, F 7. 2 , 

216H  WARHEAO  YI FL0,F7. 2/5X ,7HNAME  -  , A6 , 9H  C  .€  .P  .  ,F7  .2  ,?6H  NO 
3.  OF  CLUSTER  WARHEADS ,F 7. 2, 18H  NO.  OF  WARHEADS  ,F?  .2) 

803  FORMAT (  lHC,t*,14I5) 

9G01  FCRMAT(  //132H  *************************************************** 

^»**********««4  *+*4*+***+** *++*9* ********************* ************ 

2**************) 

END 


ooo  ooo  oooo 


SUBROUTINE  RESTNKMT,  VAL) 


—  THIS  SUBROUTINE  GIVES  THE  VALUES  OF  THE  CONSTRAINTS  ANG  C8JECT- 
- IVE  FUNCTION  - 

CQMMON/SHARE/  X(1j0),  OEL(IOG),  A ( 1  CO , 10 u ) ,N V , MC , MN ,NP1 ,NMi 
C0M(0N/0NE/NI,NJ,NII,NJJ,NM,NK,NL,NIi,NJl,NI2,NJ2,NIIJ,IEX,IE0,IEN 
lflEUtICXI, IEXFN, IifI?,I3t 14,15 ,T6 ,17, NDFFS,NC0FFS,NK,NU1,NU2,NU3 
2, IBi,ND  ,NU4,NU5,NU6.NU7,NU8,N0l,Nn2 
COHMON/T  WO/ A A (6) ,AA3(2> ,FP0P(6) ,FP0PB(2) ,01(9,7) ,FNL(10 ) , RB(13 ) , 
lEKB(lC) , TITLE 2(14) ,P( 15,7)  ,A0(U  ,15) ,6(2, 10) ,DE(2,1G) , EM (15»  , 

2E  ( 15)  ,  EK  (15)  ,  TITLrl  (1 5)  ,0  0 ( 13 , 15)  ,AN(15,  71  ,SS  (5 , 6) ,  LU  ( 1 5) 
CCMM0N/THFEF/FPEFER,PREFIN,PREINC,PREB£T,B1(9,5)  ,R0(5) 

COHMCN/FOUP/  C(15,7) ,C0 (10,15) ,ELNS(5,6) 

C0HM0N/FIVE/XX(9,15) ,0(9, 7) ,EN ( U ) , 0(9 ,5) ,U( 9 , 7 ) 

COMMCN/SIX/Z (9,7)  , 79 (2 ,1  Cl ,F ( 9 ,15) , S( 10 ) ,FB (2 , 1 C ) , BET A (6) , Y<6 ) , 
1RET  AB ( 2 ) ,YP(?),BETo,RET 

COMMON/SEVEN /CFX (9,15,15) , CFO (9,1 5, 7) , BETP (6 ) , OYX (6 , 15) ,0  EX  (6, 15) , 
10YQ(6,7) ,0BD(6,7) ,BET2P(2) , BET  FP ( 6) ,RETBPP(2) 

COHMON/PIGHT/CSX  (5  ,6) ,0  SO  (5 , 7)  ,  DFlX(lu,5,6>  ,CF2X ( 10 ,5 , 6) , 

10Fn(lu,5f7)  ,CF2P(13,5,7>  f  DY3X  (2,5,6)  ,CYB0(  2,5,7) 

C0MMCN/NINE/C88X (5,6) ,0Be0(5,7) 

IF ( MT • EG • 1)  CALL  ICEMTV 
IF (MT-l)9CC,l,i 

1  IF(Ii.EC.C)  GC  TO  2 
IF  (MT-I1)  li  0 ,11  C  ,2 

2  IF(Il.EG.NJ)  GO  TO  3 
IF(MT-I2)2iu,2un  ,3 

3  IMI2.F0.NJJ)  GO  TC  lTu* 

IF(^T-I3)3L),7CC,lGwn 

-  CONSTRAINTS  ON  SIOr  ONE  ALLOCATIONS  TO  SIDE  TWC  OEFENSES  - 


1C*  J  =  MT 
S  U  M  =  ?  . 

00  lie  1=1, Nil 
ilu  SUM  =  SUM  4  XX(I,J) 

IF(J.EO.l)  GC  TO  15^ 

IF(N02.NE.2)  GO  TO  13? 

00  120  1  =  1, Nil 

12"  SUM  =  SUM  f  U(I,J*1)  ♦  U ( I , J*2  ) 

GO  TO  19? 

13-:  00  14;  1  =  1, MI 

lU-»  SUM  r  SUM  4  U  (I  ,  J  +  NB1  -1  ) 

GO  TC  19, 

1*.  no  IRC  1=1  , Nil 
00  16)  X  =  1 ,  N  0 1 
16)  SUM  r  SUM  4  U  ( I , X ) 

GO  TC  1 9 C 

19"*  VAL  =  EM(J)  -  SUM 
GO  TC  lC;r 

.  CONSTRAINTS  ON  SIOF  ONE  ALLOCATIONS  TO  SIDE  TKC  C V  ♦  CF  TARGETS 


SUM  =  (.  . 

00  21C  I  =  1 » N 1 1 
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21:  ru*=sum4>xx<i, j) 

VAL  =  E>'(J)-SU* 
GO  TO  11CC 


- CCNSTRAIMR  CM  SIDE  ONE  ALLOCATIONS  TO  SIDE  TfcC  CV  ONLY  TARGETS 

2u"  J  =  *T 
SUN  = 

00  TIC  I  =  M1,NII 
31  **  SUW=SUM*XX  <I,J) 

VAL=EM ( J) -SUM 
GO  TC  1 J w  C 

- OBJECTIVE  FUNCTION  (K1  3FT  A8AR  -  *?.  RETA)  - 

cu'  CALL  FPACTS 

VAL  =  PRFFER*PETp-PRE9ET>  BET 
1 :  j  RETURN 
ENO 


B-ll 
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GO  TO  l?6f 

-  CONSTRAINTS  ON  SI0F  ON?  ALLOCATIONS  TO  SIOF  TWC  CV  ♦  CF  TARGETS 

20?  J=NT 

00  21t  1=1, Nil 
N=(J-1)*NII*I 
21 j  OEL(N)  =  -1. 

GO  TO  liwO 


CONSTRAINTS  ON  SI0E  ON?  ALLOCATIONS  TO  SIOF  TWC  CV  ONLY  TARGETS 


2u0  J=MT-I2 

00  3ir  1=1, NI2 
L  =  NIIJMJ-1>  *KI2*I 

310  OEL (L) =  -l • 

GO  TO  1060 

-  OBJECTIVE  FUNCTION  (K1  BFTABAR  -  K2  BETA) 

9^6  CALL  IDENTV 
CALL  FRACTS 
WRITF(F,442)BET,BFTB 
CALL  DERB 
CALL  DFRPG 

442  FORMAT  (  /I o'  ,  4HBET= ,  E 11  *4  ,  8H  BETR=,E11.4) 

IF(NU3.NE.l)  GO  TO  4C  9 
WRITF(6,438) 

433  FORMAT (  31 X , 11HA LLCCA TI ONS) 

00  AC  1  I=1,NI 

4C1  WRITE (6  >427)  (  XX  (I  ,  J)  , J  =  i , N J) 

00  412  IsKIlfNII 

402  WRITE (6,437)  C  XX < I  ,  J) ,J  =  1 , N JJ) 

IF(ND.LE.C)  GC  TO  409 
00  413  1=1, Nil 

4C3  WRITE(6,477MU(I  ,K),K=1,ND) 

409  DO  *16  J=1,NJ 
00  416  1=1, NI 
N=(J-1)*NII*I 

410  0?L(N)  =  PREFER  *  CB9X(I,J) 

DO  43C  1=1, NI2 

00  426  J=1  ,N J 
N= (J-l) *NII*I*NI 

42 J  DEL (N)  =-FR*e?T  *  C9X(I,J) 

DO  436  J=1,NJ2 
N=NIIJ+(J-1)*NI2+I 

430  DEL (N)  =-PRFBET  *  0RX(I,NJ«-J) 

-  IF  NO  ATTACK  CN  OFFENSES  BYPASS  THIS  SECTION  — 

IF(NC.LE.C)  GO  TO  100  G 
DO  45G  K=l,NO 
00  441  1=1, NI 
N=  IEX  ♦  (1-1)  *  ND  «■  K 

44?  OEL  CN)  =  PREFER  *  0BBD(I,K)  *  B1(I, K)  *  0(1, K) 
DO  45 j  I=NT1 ,  MI 
N=IEX  ♦  (1-1)  *  ND  ♦  K 

450  OEL (N)  =-PREBFT»  DBO(I-NI,K)  ♦  B1(I,K)  ♦  0(1, K) 
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SU3R0UTINE  GRAC1  (MT) 


:  ,  :  !  I  ! 

-  THIS  SUBROUTINE  CALCUl AT ES  '  THE  GRADIENT  VECTOR  CF  THE  CONSTRAINTS 

- ANO  OBJECTIVE  FUNCTION - ; 

i 

CO^MCN/SFARE/  X(1G;),  DEL(lj5>,  A  (1 00  ?  l'j  C  )  ,N  V  ,*C  *  HN  ,NP1  ,N*1  1 

COMMON/ONE/NI,NJ,NII,NJJ,NM,NK-,NL,NIl,NJl,NI2,NJ2,NIIJ,IEXi,IED,IEN 
1,IEU,IFX1,IEXEN,I1,I2,T3,I4,T5,I6,I7,NOEFS,NOOFFS,NN,NU1,NU2,NU3 
2,IQl,N0-,NU4,NU5,NU6,NU7,NU8,Nm,N0? 

C0MM0N/TW0/AM6)  , A A n ( 2) ,FPCP (6 ) , FPQPR ( 2) ,0L(9,7) , ENL (1C ) , RB(ir) ,  j 
1FK0(1C> ,TITLE2UC) ,P(15,7)  j  A9  (10 ,1 5)  >  B  ( 2 , 1G  ) ,00(2,10 ,FM(15) , 

2E (15) ,EK(15)  i  TITLE1  (1 5)  sFPdD*i5>  ,ANC1?,7,)  ,SS<5 ,6)  -  LUC  15) 

COMHON/ THF EP /pc>EFER,PREFIN,Fc»E INC, PRE9ET,  Bl(o,5)  ,R0(5)  s 

CO^MCN/FCUR/  C  (1  5 , 7)  , CB  (It  ,15)  ,  El  Nc  (5,6) 

COMMON/ FIVE/XX (9, 15) ,0(9,7) ,EN < 10 ) , 0(9 ,5) ,U(9,7) 

CCMMON/SIX/Z (9,7 )  ,ZP(2,1C) , FT  9,15)  ,S(1C) , FP (2 , 11 ) , EFT A  (6 ) ,Y ( 6 1  , 
1RETA9(2) ,YB(?) ,8ETR,BET 

C0My0N/SEVEN/CFX(9,15,15)  ,0*0(9,15,  7) , BETP (6 ) , CYX (6 ,15)  ,0eX(6,15)  , 
10X0(6,7) ,CEO  (6,7) ,6FT8P(2) ,RETPP(6) ,8ETBPP(2)  1 

COFMON/EIGFT/DSX (5,6) ,0SD(5,7)  ,  DF1X<10,5,6) ,CF2X( 10,5, 6) , 

lOFl'Hlf  ,9,7)  ,CF?P(l'3,S7)  ,  >  0,Y8X  <2 ,5 ,6)  ,CY8C(  2,5,7) 

COMHON/NINE/08BX  (5,6)  ,0BB0(5,7) 

COMMON/ EL FVFN/CRXX (9C ,15)  ,00X0 (9u ,7) ,0000(30  ,5) 
C0MM0N/TWrLV£/DSXX(3':  ,6)  ,OSXD(  TO, 7^  ,OSOO(6,5) 
CONMON/FCURTN/CRBXX(3J,30 ) ,008X0(30 ,35 ) , O80DC ( 31 ,3? ),02YB(2)  i 
437  FORMATt  KF13.3)  ,  , 

00  1C  1  =  1, I** 

\>  DEL ( I )  =  I • . 

IF(MT-l)Cr;,l,l  ! 

1  IF  ( 1 1 • EO • w)  GO  TO  2  ' 

IF(MT-Il) ICO  ,lcC  ,2 

2  IF(U.EO.NJ)  GO  TO  7  ,  1 

IF(MT-I2) 2  J  ,  2wn  ,3 

3  IF(I?.FO.NJJ)  GO  TC  lali 

IF ( MT-I 3 ) 3Cj,3uu  ,  1°  uC  > 

i  i  |  i 

.  CONSTRAINTS  ON  SIOE  ONE  ALLOCATIONS  TO  SIOE  TWO  OFFENSES  - - 

dfl  J  =  M  , 

00  lie  1=1, Nil  ‘  1 

L=(J-l)*NIl4l  , 

113  0EL(L)=-1. 

IF(J.EO.l)  GO  TO  15^  s 

IF(ND2.NE.2)  GO  to  13" 

00  12C  1  =  1, MI  , 

00  12i  K= 1 ,NC2  | 

M=  I  EX  4  (1-1)  ♦  NC  +  U  J 

120  PEL(M)  =  -1.  1 

GO  TO  l‘i *  1 

11-  no  1 4i  1=1, Nil 

M=I EX  4  (1-1)  *  NO, ♦  J  4  NHi  -  1 

14C  OEL(M)  =  -1. 

60  TO  i:j  1 

153  00  IE  T  1  =  1, Nil 
00  16 J  K  =  1 ,N°1 
M=  I  EX  4  (1-1)  ♦  NC  4  K 

165  DEL(M)  r  -l.  1 


,  B-13  Cwac> 


1 


GO  TC  loaf 
-0C  RETURN 
ENO 
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o  o  o  o  olo  o  o  o  o  o  o 


COMMON/SEVFN/OFX  (9,15  ,15)  ,0^0(9,15,7)  , BETP(fc)  ,nYX(6,15>  ,0?X(6,15) 

10Vn(6,7)  ,rfi3(6,7),g|:TBP  <2),9£TPp(6)  ,BETBPP(2) _ 

COM'iCf /cTGHT/0SX  (5,6)  ,DSD  (5  ,7)  ,  OF IX ( 1  j  ,  5 , 6)  , CF2 X ( IP  ,5 , 6)  , 

lQFn(iu,5,7)  ,  0F2P(1%5,7)  , _ _ DYEX  (2,5,6)  ,PYBP(2,5,7) 

C0MMGN/NINF/nE3X  (5  ,6)  ,0380(5,7 ) 

CO^MOt'/FLFVEN/DRVX(9o  ,15)  ,PBXO(9ut 7)j  3  BOPJJmjl 5 > _ _ 

COMMCN/TWeiVE/CSXX  (3u  ,6  )  ,  DSX9  (  3C~,  7)  ,DS09  (  6  ,~5  ) 
CQMH0N/F0URTN/0QPXX(3g, 31 ) ,QP 9X0(30 ,35) , 08900 (3u » 30 ) ,02 YQ (2) 
IF(MT-i)9' v ,1,1 
1  IF(Il,tQ.L)  GO  TC  2 _ 

if(^t-ii) ir j,i:: ,? 

?  IF ( 1 1 , £Q , N J)  GC  TO  3  _ 

IF(MT-I2)2i  j,2^,3 
3  IF  (t?«EQ«NJJ)  GO  TO  lu 
IFCMT-T3) 3U  ,3CC  ,1wJL 

. CONSTRAINTS  CN  S’l9E~0NC  ALLOCATIONS  TO  SIDE  TWO  oefenses - 

T__.--.__ 

GO  TO  1H? _ 

.  CONSTRAINTS  ON  SIHE  ON^  ALLOCATIONS  TO  SIDE  ThC  CV  +  CF  TARGETS 

2  C:*  IK  =  1 _ _ _  _ _ 

go  to  i-  -wr 

. CONSTRAINTS  ON  SI HE  0Nr  ALLOCATIONS' TO  SIDE  TfcG  CV  ONLY  TARGETS 

_____ 

GO  TO  lwC: _ _ _ _ 

-----  OBJECTIVE  FUNCTION  (K 1  BETABAR  -  K2  BETA)  - - 


9  j >  CALL  I D  c  N  T  V 
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o  olo  I  i  i  II  '  looln  oj  I  I  i  !  !  olo  ol  S  !  I  oloo 


LCPL=(J-1)*NIT*K 

9 ;  2  A (L^CWjLOCL) =  PREFER  »QPrUX( NN?yr >  _ 

- IF  NO  ATTACH  CM  Dc.FcNSE S  BYPASS  THIS  SECTION 


IF(NC.Lr.l)  GC  TO  91j 


OC  9.F  1 1  -  1 ,  M 

_ 0  0_9r f  K  =  1,^0 

MN?=TlI-l)*Nn  4-  K 

__  LC OL  =  T rX  4-  NN?  _  _  _ _  ___  _ 

9’f  A(L*?OK,LCCL)”s  p^EFEP**  OrfBXO(KN?",NNJ)  '^”1(11,0  ♦  0(11,0" 
91"  CONTINUE 


---  IF  NC  ATTACK  OK'  DcrF.NSrS  BYPASS  THIS  SECTION 


_ IF(Nn.Lr.")  CC  TC  £1F  _  _  _  _ 

nc  91?  I  =  t , N I 

OU  91?  J=l,Nr _ _ _ _ _ _ _ 

JN*>  =  (T-l)  *  MO  +  J* 

_ l  ROW  =  irx  4-  NNg  __  _  _ _  _ 

10  11?  K=  I , NI 

_  10  912  KP=J,NO  _ _ _ 

NN3”=  (K-l)  ~+  Nl";  KR  " 

LCOL  =  TFX  4-  NN3 _ _  _____  _ _ 

TERM  r  '7 

_ If.LL._FO  . K.  AMO i .  J .  KD>  Tjf  91  =  ORB  0  C I ,  J  >*R  1  (I  ,  J )  f*2fC.  _ 

-j--_JHFSF  ARE  SfXCNO  PARTIAL  U  C?OSS  U  DERIVATIVES  WHICH  DIP.  NOT 
---  A F F £  A~R  ~  T N  "  P  A C •  C 9 •  £ ' 

"9 1?  a ( l 90W , lcTcD  ,jr*DTKri^r»lIT!TJ)^ 

_ 1KR)_4-_T-9M) 

'919‘Tfc  ’RELET  .EQ.n  GO  T 0  ic:f 

in_92S  _JJ=1  ,NJ  __  _  _ 

"n io  9? 9"  i^r,KiF  ~ .  ’ 

LPQW=<JJ-1  )*NII4»NI  +  I  _  _  _ _ 

4N?-<JJ-2  )*NT2%I~~ 

_ 00  n? _j=jj*nj  _  _  _  __  _ _ 

LCCw=  (  J  - 1")  ♦NT  I’T'  mi'  4- 1 ' 

917  A_( Lf C W i_L CC L)  -  -  PRfB^T  *_  n^yy  (  NM? ,  J)  _ _ _ 

’""DO  919  J  =  1  » N J?  . .  .  . . . 

LCOL  =  N 1 1 J  »  (J-l)  *  NI?  4-  I  _  __ _ 

919  AtLROPfLCOL)'  =  -  PRESET  *  QBXXl  NN?,  N  J+  j) . 


31  (K  , 


- IF  NC  ATTACK  CN  CCFENSE S  BY  c AS  S  THIS  SrCTICM 


IF(Nn,LE.;)  GO  TO  9?B 
10  9?'  K  =  1 , N 0 


LCOL  =  IfcX  f(Nl4-I-l)  *  Ml  4-  v 

92  v  ACLROW,  LCOLi_=...-FPtfri.  *  DPXD(  NN2,K)  *  ?  1  <  M 1 ♦  I  > K  L*  L  (NI»I,K) _ 


9?3  CONTIM»c 

_ 00  93,  JJ’1,MJ2 _ 

GO  93:-  1=1, NI2 

_ LPQH  =  hIIJv (JJ-1) *NI?»I _ _  .  _ 

^N2=( JJ-l) *NI2fI+Nj 

_ JO  926  JrJJ,NJ2 _  _ _  _ 

LC0L-NIIJMJ-l>n!I2*I 

92  F  A(L?OW,LCCL)=-PR£ObT  *  POXX ( NN2»NJ+J) _ 

c 

C - IF  NC  ATTACK  OK  P>rFFNSF  S  8YPASS_  THIS  S  EjCI  I  CM -  - 

C 

_ I F  (  >IP.L:-.u)  00  TO  9  or _  _ _  _ 

00  923  <=1,N> 

_ LCCL  -  ICX  »  tNI+I-1)  *  N'D  +  K _ _ _ 

928  A(LPOWfLCOL)  =  -PPFlt'T  *  OBXO < NN2  ,< >  ♦  *1  <  NH- 1 ,  K )  *0  (  NI+ 1  ,K  ) 

_ 93  CCMIKU~ _ _ _ _ 

C 

C - IF  NC  ATTACK  QM  r£F£Nf;bS  PYPASS  THIS  StCTJON  -------- _ 

C 

_ IF(Nn>L~^)  GC  T*  lCw-I _ _  „  _ _ 

00  94^  I  =  M1,NIT 

_ PC  94;,  KKsjjNC _  _ _  _  _ 

L°OW  =  TFX  +  ( I- 1 )  *  Vr  f  K< 

_ NN2  -  (1^ J  )_^  N0_4-_KK  -  NO  *  NI  _  _ 

00  9^"<=Xk,N0 

_ l r°L  -  i rx  »  n-n  »  Np  ♦  k  _ _  _ _ 

JFR‘*  =  T. 

_ IFtK.LC.KK)  TEPM  =  3BD(I-NI,K)  *  PltI,K>**2  *  0  OiJ< )  _ _ 

C 

C - T  H  E. c  E  ARh:_S9CCN0  PARTIAL  l  CROSS  U  _CER  T.V  A  T IV  FS  M~JCH_DIC  NOT 

C - t'm  RA C -C°-43 

C _  _ _  _ . 

933  A  (L9r)W,LCCL>  =  -  DRr9rT  ♦  (  }P0 0  ( NN?  , K)  +6  ( I , K  )  *0  <  I , K X )  *8 i ( I , K)  ♦  «1  { I 

_ l,**)  +  T rR ^ )  _ 

9  4  CCN’TINU- 


o  o  o 


SUBROUTINE  I 0 ENT V 


C 


C - THESE.  COMTUTICNS  APE  THE  SAM^  AS  TH OSfLQ ESCRIBED  IN  RAC-CS-43  PART 

C  —  TV-9.  WITH  EXCFPTIONS  AS  NOTEO  IN  THL  COMMENT  CAROS . - 

C _ _ 

COKMCN/SH APE/  X(lCu),  OEL(l*jy>,  A  ( 1 i  C  ,  lu  3  )  ,N  V  ,MC  ,  MN  ,  NP1 ,N*1 

_ COMMOy/nNF/NI«K,JtNII>NJJtNM>N<fNL,NIl,NJltNI2tNJ?tNIIJt  IEX,I£Q,IEN 

i,IEU,IEXl ,IEX£N,Il,I2,I2,l4,I5,I6,I7,N0EFS,NC0FFS,NNfNUl,NU2,NU3 

2, 1 81, NO  ,NU4,NU5,NU6,NU7,NU8,NQ1,N02 _ 

COMMON/TWO/ A  A C6I  ,AA1(2)  ,FP0P<6)  ,FP0P3(2)  ,DL(9,7>  ,  ENL  (13  ) ,  RB  (1C  >  , 
lEKO(li)  ,TITLf2(lw)  J (15  >7)  ,A3(1C  jJ.5)  ,D(e,lJ, 06(2,10  ,FM<15)  , 

2t ( 15) ,FK(15)  ,TITL£1<15)  ,P3(1C ,15) ,AN(15,7) ,SS<5,6> ,tU(15) 

_ COMMON/ THREE /FREE £R, P PI N ,PRE INC, PREP ET, 61(9,5) ,R0 (5) _ 

COMMON/FOUR/  C  (1  c ,  7 ) ,  C0  (1 0 , 15) ,FLNS<5,6> 

_ C0MMQN/FIVE/XX(9,15),C(9,7) ,EN(lw),Q< 9,5) ,U(9,7) _ 

COMMON/ST X/Z (Q,7)  , ZB (2 , lu )  ,F ( 9 , 1 5)  ,S ( it )  , F0  ( 2 , 1 C ) , E^T A  (6 ) ,V (6 )  , 

13ET  AP (2 ) ,  YB(?) ,3FTa,PET _ 

COMMON/SPVFN/CFY  (9,15,15)~,DFD(Q,15,7)  ,PETP(6)  , DYX (6 ,15)  ,OBX (6 , 15)  , 

lQV 0 ( 6, 7 ) ,060(6,7) , PET  BP (2) ,PETPP(6) , PE  TBpP (2  ) _ 

COMMCN/SIGHT/PSX  (5,6)  ,0S0(5’,7)  ,  0F1X(  lu  ,5  ,6)  ,0F2X(  IP  ,5 , 6)  , 

13F1JUC  ,5,7)  ,PFc!P(U,5,7)  , _ _ DYEX  (2,5,6)  ,CY60(?,5,7) _ 

COMMCN/NINE/OEEX  (5,6)  ,PPB0(5,7) 

C  **  THI S  SUBROUTINE  In  FNT  IF  IPS  MODEL  J/_ARI  APLJES_  HITH  SU^T  VARI A3LFS  X  XX 
00  ir  J  =  1  ,MJ 

_ JJ=(J-1) »NII _ 

00  U  1  =  1, Nil 

_ K  =  JJ»I _ 

i'j  XX(I,J)=Y(K) 

_ 00  2 v  J  =  1,NJ2 _ _ _  _  _  _ 

JJ  =  NI T  JM  J- 1)  *NI ? 

_ 00  2.  I  =  M1,NII _ 

K=JJ+T-NI 


L  =  KJ  ♦  J 


?r  xx  (i , l )  =x  (x) 


- IF  NC  ATTACK  ON  TTFENSES'  9Y  PASS  "THIS"  S  FCT  I  C~N 


I F  (  N D  •  L E  •  ‘  )  GC  TO  6r 

_ DO  3:  1  =  1, Nil _ _ _  _ 

C 

C  --  NC  REPLACES  N K  IN  EQUIVALENT  INSTRUCTION  OF  RAC-CR-45  PART  IV-O. 

C  “  . . . 

_ II  =  (1-1)  *  NO  »  IFX _ _ _ 

DO  3 r  K  =  i,~MO 

_ L  =  n»r _ _  _  _  _  _ _ 

3:  U(I,K)  =  X ( L ) 

6:  return _ _  _  __  _ _ _ 

CNO 


o  o  o  oooo  oooo 


SUEROUTINF  FRACTS 


—  THESE  COMPUT I ONS  ARE  THE  SAME  AS  THOSE  DESCRIBEO  IN  RAC-CR-43  PART 

—  IV-B.  WITH  EXCEPTIONS  AS  NOTED  IN  THE  COMMENT  CARCS  . 

CCMMCN/SHARE/  X(luc),  OEL(l?C),  A (100,  lu 0 )  ,NV  ,MC  ,MN  ,NP1  ,NM1 
COMMON/ONF /NI ,NJ  , NI I , NJ J, NM , NK , NL , N II , NJ 1 , NI 2 , N J? ,N II J , IEX , IED , IEN 
1, IEU,IEX1,IEXEN,  11,12,13,  I4,I5,I6,I7,NDEFS,NC0FFS,NN,NU1,NU2,NU7 
2, IB1 , ND  , NU4 ,NU5 ,NU6, NU7 ,NU9,N0i , NQ2 
COMMON/ TWO/ A A (6 ) ,AA3t2) ,F  FOP ( 6 ) ,  FP0p9  (  2) , OL ( 9 , 7 ) , ENl (1C ) , RB (1 3 ) , 
1EK9  ( 1C  )  ,TITLF2(lw)  ,P(15  ,7)  ,AR  (1C  ,15)  ,  0<?  ,1j  )  ,OE  (2, 1G)  ,EM(15)  , 

2E ( 15) ,EK ( 15) , TIT  LEI (15)  ,PB(li,15) ,AN(15, 7) ,SS (5 ,6) , LUC  15) 
COMMON/THREE /PREFER,PRF FIN, ORE INC, PREBET ,81(9,5) ,RO  (5) 

COMMCN/FOUP/  C(15,7) , CB (1 0 , 15) , EL NS (5, 6) 

COMON/riVE/XX  (9,15)  ,0(9,7)  ,EN  (1C  )  ,0(9  ,5  )  ,U<  9 ,7) 

COMMON/S  IX/Z  (9,7)  ,7B(2,U>  ,F(9,15)  ,S(1G)  ,  FB  (  2, 1C  )  ,  BET  A  (6 ) ,  Y  ( 6 )  , 
13ET  AR ( 2) , YP(2) ,BETB,BET 

C0MM0N/SFVEN/CFX(9,15,15)  ,OFO(9,15,7) , EETP (6 ) ,DYX (6 ,15) ,0€X(6,15)  , 
10Y 0(6,7) ,0B0 (6,7), BET  BP (2),BETPP(6)  ,BETBPP(2 ) 

COMMON/EIGHT/CSX (5,6) , PSP (5,7) ,  DF 1X110,5,6) ,0F2X( 1C ,5 , 6) , 

1QF10 ( 10,5,7) , CF2  C ( 10 , 5, 7) ,  OYBX (2,5,6)  ,CYBP(2,5,7) 

COMMON/NINF/DBEX (5 ,6) ,0890 ( 5, 7) 

*  COMPUTE  FRACTIONS  F,S,AND  FBAR  , BET  A , PET AP A R ,  ALSC  Z  ZEAR  Y  YBAR 
— -  IF  NO  ATTACK  ON  DEFENSES  BYPASS  THIS  SECTION  - 

IF(ND.LE.O)  GO  TO  7 
00  5  1=1, Nil 

- NO’  REPLACES  NK  IN  EQUIVALENT  INSTRUCTION  OF  PAC-CP-4?  PART  IV-B. 

00  5  K= 1 , N  C 

5  0(1, K)  =  OL ( I , K)  *  EXP( 91(1 ,K) *U(I,K)  ) 

7  00  30  J  =  1 , N J 
00  30  1  =  1, NI 
SUM=0. 

00  2 u  K  =  1 , NK 
Z(I, K)=7. 

00  IB  L  =  1 , N J 

10  Z ( I  ,K) =  Z ( I , K )  ♦  XX  (I  ,L  )*  AN  (L ,  K ) 

IF  (Z(I,K) .LE.C.) Z(I,K)=1.E-1C0 
2r'  SUM=  SUM+  C(J,K)  *0  ( I ,  K)  /Z(I,K) 

I F (SUM  .L E  .-14C • )SUM  =-140. 

I F ( SUM  .GT .14t.) SUM  =  14". 

30  F ( I , J ) =  EXP ( SUM) 

00  6?  1  =  Ml, Nil 
00  6C  J=1 , N J J 
SUM=0. 

DO  5C  K  =  1 , N K 
Z(I,K)=n. 

00  40  L  =  1  ,NJ J 

40  Z(I,K)=Z(I,K)*  XX(I,L)*AN(L,K) 

IF(Z(I,K)  .LF.C.) Z(I,K)  =  l.E-ltO 
5r-  SUM  =  SUMfC(J,K)*0(I,K)  /Z«I,K) 

I F ( SUM  .LF.-14C.) SUM  =-140. 

IF (SUM  .GT.14u.) SUM  =140. 

6 J  F (I,J)=EXP (SUM) 


B-19 


00  ICt  1=1, NI 
SllM  =  iu . 

00  9;  J=1,K*J 

9**  SLf*=SUM  +  XX  (Tf  J)*F(I,J)*ELNS  (T,  J) 

S.  =SUtf/EN<I> 

':F(SUM  .LE.-14i.)SLM  = •  1 4  n  # 

IF  (SUM  .GT.14C.)  SUV  =  14C. 

1-'  Sil)  =  RR<I)  *  F VP (SUM) 

00  14C  L  =  1,M 
00  14c  M  =  1  ,NM 
SU*2  =  t:  . 

00  l7r  N= 1 , NN 
SUV1=C. 

00  12l  K= 1 ,  NM 

12 J  SUM1=SUM1*S(K)»EN(K)*A*MK,N) 

ZR(L,N)=  FFOFP(L>*  SUM1 
IF(ZB(l,N)  .LE.C.  )ZB(L  ,W)  =1*  E-l?  r* 

1? J  SUM?sSUM2+CB(KfN)*Ce<l,N) /Z9(L ,N) 

IF  (SUM2#LF*-14°#)SUV2=-140* 

IF  (SUM2.GT.14U.) SUM2=14C. 

14;  FB  (L,M) =EXF(SUV?) 

8£T  =  u. 

00  2&w  1  =  1, M2 
Y<I)=r. 

no  15C  J=1,NJJ 

150  YmsYCI)*FCM  +  I,J)*XX(NI*I»J)  *FK(J) 

I F ( Y ( I ) . LF.O • I YC I) =!• E-l^  0 
TERMsAA(I)*SORTCV(I)) 

3ETA ( I) =1  .-<i.*TERM)*EXP(-TEPM) 

2u 0  3 E T r B F T  ♦  FFCPt I ) *  ErT  A ( I) 

B  £  T  9  =  C  • 

00  3„E  L= 1 , NL 
YB (L)=C  , 

00  25i  Ms  1 , NM 

25,  YB (l) =YR  (L  )  ♦  FP(l,V)*S< V) ♦EN(M) *rK9 (M) 
IF(YP(L).LE.?'.)YP(L)=l.F-lw3 
YO  (L )  =  FFOFB(L)  *  Y°(L) 

TFR^=  AAB(L) *SQR  T ( Y B ( L I  ) 

RFTAEi(L)  =l.-(l.fTFRM)  *EXP  (-TERM) 

3l  r  3ET0-PETB  ♦  FFOPB ( L ) *9F YA p< L) 

RETURN 

ENO 
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oooooo  o  o  o  o  o 


SUBROUTINE  DFRB 


---  THESE  COMPUTICNS  ARE  THE  SAME  AS  THOSE  DESCRTBEC  IN  RAC-CR-43  PART 
— -  IV-B.  WITH  EXCEPTIONS  AS  NOTEO  IN  THE  COMMENT  CAROS  - 

CGMMCN/SHARE/  X(1LG),  DELUdC),  A (1 jC , 1C r ) ,N V ,MC ,MN ,NP1 ,NM1 
COMMON/ONE/NI,NJ,NII,NJJ,NM,NK,NL,NIi,  NJ1 ,NI 2 ,N J2 ,NIIJ , IEX, IEO , IEN 
1,IEU,IEX1,IEXEN,I1,I2,T3,I4,I5,I6,I7,NDEFS,NCCFFS,NN,NU1,NU2,NU3 
2,101, NO  ,  NU4 »NU5 »NU6»NU7»NU8»N01»N02 
COMMON/TWO/ AA (6) ,AAB<  2) ,FP0P(6) ,FPOPP( 2) ,0L ( 9 ,7) , EN L (1C ) , RB(1Q ) , 
lEKB(lC)  ,TITLE2(1C)  ,P(15  ,7)  ,AB  ( 10 ,15) ,  0(  2 , 10  > ,  00  (2, 1C)  ,EM(15)  , 

2E ( 15) ,FK(15) ,TITLE1(15) ,PB<10,15) , AN< 15, 7 ) ,SS 15, 6) ,  LU( 15) 
COMMCN/THREE /PREFER, PRt FI N,FRE INC, PREBET ,01(5,5) ,RO (5) 

COMMCN/FCUR/  C<15,7),CR(lti,  15)  ,ELNS(5,6) 

C0MM0N/FIVE/XX(9,15)  » C( 9, 7)  ,EN  (1"  )  ,  0(9  >5  )  ,1(9,7) 

C0M0N/SIX/Z(9,7)  , ZB  (2,10)  ,F(9,15)  ,S ( 10  )  ,F0  ( 2, 1 C ) ,  EET A  (6 ) , Y(6 ) , 
13ETAB ( 2) ,  YE ( 2) , BETPjBET 

CO«MON/SEVEN/DFX (9,15,15) ,DFO( 9,15, 7) , BETF (6) ,0YX(6,15) ,0BX(5,15) , 
10Y0((>,7)  ,rBO(6,7),0ET9P(2)  ,BETPP(f>)  ,BET9PP(2) 

COMMON/ EIGHT/CSX (5,6)  ,p$0  (5,7)  ,  0F1X( 10 , 5 ,6) ,CF2X( 10 ,5 , 6) , 

10F10(1C ,5,7) ,0F2P( 10 , 5, 7) ,  nY0X  (2,5,6)  ,CYBC(2,5,7) 

COMMON/NINE/Oeex (5,6) ,0000(5,7) 

*  COMPUTES  FIPST  CE° IVATI VFS  OF  BETA,  ALSO  DFRS.  OF  F  *** 

90  3'  1=1, NI 
90  3C  J=1,NJ 
00  3l  L=1,NJ 
SUM=0. 

00  2-  K=1,NK 

20  SUM=SUM,C<J,K)*n(I,K) *AN(l,K)/Z(I,K>**2 
3"  0FX< T,J,l)=-SUM*F(I,J) 

00  6C  I  =  M1,NTI 
90  6G  J  =  1  ,NJ J 
90  6C  L  =  1  ,NJJ 
SUM=i . 

00  5L  K=1,NK 

5"  SUM=SUM*C  (J,K)*n(I,K)*AN(L,K)/Z(I,K)**2 
60  OFX(I,J,L)=-SLM*F<I,J) 

-  IF  NO  ATTACK  ON  OEFCNSFS  BYPASS  THIS  SECTION  - 

IFOO.lE.o)  GC  Tr  85 
90  7’J  1  =  1, M 

—  -  NO  REPLACES  NK  IN  ECU IV  AL  ENT  INSTRUCTION  OF  RAC-CR-43  PART  IV-B. 

00  71  J*1,NJ 
00  7.  K= 1 , NO 

7"  9F9(I,J,K)=F(I,J)*C(J,K)/Z(I,K) 

"0  9.  I  =  Ml, Nil 

oo  a:  j = i , m  j  j 

00  8C  K= 1 , NO 

8"  OFB(I,J,x)=f(I,J)*C(J,K)/Z(I,K) 

85  90  1  =  1  .M? 

9G  BET“(I)=AA(I)**2*EXP(-AA(I)*SQRT(Y(D) )/2. 

DO  121  1=1, NI2 
90  12E  J= 1 ,N J J 
SU“  =  C . 
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11" 

120 


DO  Ilf  K=1,NJJ 

SUM=SUN*XX(NI*I,K)*CFX<NI+I,K,J) 
OYX  (I,  J)  =SU^  +F(  M*I*  J)  *EK(  J) 

08V (I , J)=FFOP(I) ♦8£TP (I >VDYX(I , J) 


♦  EK(  K) 


---  IF  NC  attack  on  cefenses  bypass  this  section 

IF  (ND.LE.G )  GC  TO  15C 


—  NO  REPLACES  NK  IN  FCUIVALENT  INSTRUCTION  OF  RAC-CR-1.3  PART  IV-B 

00  1<*C  K  =  1,N0 
DO  14G  1  =  1, NI2 
SUM=;. 

DO  133  J=1 ,NJ J 

13?  DYO(ll'K)+=SU«ItI’J>#CFn,NI  +  I’J’,<>  *  FK(J> 

14c  OBD(I,K)=FPCP(I) ♦8ETD (I ) *0Y0 (I fK) 

150  RETURN  * 

END 


oooooo  oooo 


SUBROUTINE  OFRBB 


—  THESE  CQMPUTIONS  ARE  THE  SAME  AS  THOSE  OESCRIPEC  IN  RAC-CR-43  PART 

—  IV-B.  WITH  EXCFFTI CNS  AS  NOTEO  IN  THE  COMMENT  CAROS  - 

COMMON/SHARE/  X(1CC),  DFL(lOi),  A ( 1 OC , 10 ? ) ,N V ,MC ,MN ,NP1 ,NMi 
COMMON/ONP/NI ,NJ ,NII , NJ J, NM ,NK , NL ,NI1 , NJ1 , NI 2 , NJ? ,N II J , IEX, IEO , IEN 
1,IFU,IEX1,IEXEN,T1,I2,I3, 14 , 15 , IE ,1 7,N DEFS , NCCFFS , N N ,NU t, NU2, NU3 
2f IBl,NO  ,NU4,NU5,NU6,NU7,NU8,NDi,ND2 
COMMCN/TWC/AA (6) , AA9( 2) ,FPOP (6 ) , FFOPR ( 2) , DL< 9 ,7) , FNL (10 ) , R8(10 ) , 
lEKR(li) ,TITLE2C1C) ,P(15  ,7)  ,AB (1C ,15) , B(2 ,1C ) ,06 (2, It) ,FM(15) , 

2E ( 15) , EK ( 15) ,TITLE1(15) ,*>B ( U , 15) , ANC 1 5, 7) ,SS (5 , 6 ) , LU( 15) 
COMMON/THREE/FREFEP, PREFIN, PREINC ,PREP^T , 91(9, 5)  ,RQ  (5) 

COMMON/FOUR/  0(15,7) , CP (lu, 15) ,ELNS<5,6) 

C0MM0N/FIVE/XX(9,15) ,0(9,7) ,EN (1C) ,Q(9 ,5) ,1(9,7) 

C0MM0N/SIX/Z(9,7)  ,ZB  (2 ,1 w) ,F( 3,15) ,S( 10 ) ,FP (2, 1C) , EFTA  (6? , Y(6) , 
1BFTAEH2 )  ,YP  (2)  ,DFT6,BET 

COMMON/SEVEN/CFX (9,15,15)  ,0F0(9,15,7) , *ET F (6 ) , OY X  (6 , 15 )  ,0  6X(6, 15)  , 
10Y0(6,7) , CRD (6, 7), BET 6P (2) ,8ETPP(6) ,BETBPP(2) 

COMMON/ EIGHT /CSX (5,6)  ,OSO(5,7)  ,  DF1X(  lti  ,5  ,6)  ,0F2X( 10 ,5 , 6) , 

10F10(10,5,7) ,CF20(10,5,7) ,  DYBX (2,5,6) ,0YBD(2,5,7) 

C0MM0N/NINF/0P8X (5,6) ,0890(5,7) 

00  8C  M=1,NI 
00  J=1,NJ 
SUMrj. 

00  10  L=1 ,NJ 

1?  SUM=SUM+XX(M;L)*CFX(M,L  ,  J  )  *ELNS  (M  ,L  ) 

SUMrSl!M  +  F(M,  J)*ELNS(M,J) 

20  OSX(M,J)=S(M) *SUM/EN(M) 

---  IF  NO  ATTACK  CN  CEFENSES  BYPASS  THIS  SECTION  - . - 

IF(ND.LE.O)  GC  TO  60 

—  -  MO  REPLACES  NK  IN  ECU IV  AL ENT  INSTRUCTION  OF  RAC-CR-43  PART  IV-9. 

00  4C  K= 1 , NO 
SUM=Q. 

DO  3c  L  =  1 ,  N  J 

30  SUM=SUM4-XX(M,L)#CFC(M,L  ,K)*ELNS(M,L  ) 

4G  OSC(M,K)=S(M)»SUM/FN(V) 

80  CONTINUF 

00  12c  L=1,NL 
no  12C  M=1,NM 
00  120  1  =  1, NI 
00  12C  J= 1 , N J 
SUM=W. 

00  110  N  - 1 , NN 

11"  SUM=SUM^CB(M,N)*CB(L, N) *AB(I,N)/Z9(L,N)**2 
SUM=-FP0P6(L) *FR(L,M) *DSX (I , J) *EN (I )  #SUM 
GO  T0(111,112),L 

111  DFlX(M,I,J)=Sl)M 
GO  TO  1 2D 

112  0F2X(M,I,J)=SUM 
12 j  CONTINUE 

IF(ND.LE.C)  GC  TO  195 
00  14C  L  =  1 , NL 


ooo  ooo  ooo 


00  141  K=1,NM 
no  14*  1  =  1, NX 

---  NO  REFLACFS  NK  IN  EQUIVALENT  INSTRUCTION  OF  PAC-CR-43  PART  IV-E. 
* 

DC  140  K=1,ND 
SUM=C. 

DC  13C  N=1,NN 

13:  SUM=SUM+CB(M,N)*CB(L,N)  *AB(I,N)/ZB(L,N)**2 
SUM=-FP0PP(U *FR(L,M> *PSD (I,K> *FN(I)  ♦  SUM 
GC  T0(13l,13?) ,L 

131  QF10(M,I,K)=SUW 
GO  TO  14* 

132  DF20(M,I,K)=SUM 

14^  continue 

195  DO  231  1=1, Nl 
00  2"*'-  1  =  1, NI 
DO  2?  t  J=1,NJ 
SUM  =  '-. 

DC  22C  M=1,N* 

GO  TO  (211, 212), L 

211  TFMP=0F1X(M,I,J) 

GO  TO  22P 

212  TE.MP  =  DF2X  (M,  I ,  J) 

220  SUM=SUM+TFMP*S(M)*EN( *) *PKR(M) 

SUM=SUMfFF(L, I)*0SXCI ,J )*ENU)  *EKB(I) 

23 j  0Yex<l,I,J>  =  FP0Pn<L)*SUM 

-  IF  NC  ATTACK  CN  CEFENS^S  BYPASS  THIS  SECTION  . - 

IF(NC.LE.C)  GO  TO  295 
GO  26C  1  =  1, NL 
DO  26 J  1  =  1, M 
DO  260  K=1,NC 
SUM  =  i, . 

DO  25:  M=1,NM 
GO  TO  (241, 242), l 

241  TEM?  =  DF1D(M,I  ,K) 

GC  TC  25C 

242  TFM?=DF?0(M,I,K> 

25:  SUK=SUM*TCMP*S(M)»EN(M) *FKB  (M) 

SUM=SUM*FB(L,I)#DSC(I,K)*EN(I) *FKB< I) 

26*  OYBD(L,l,K)=FFCPF(L)*SUM 
295  DO  3  j (  L  =  1,NL 

3 )  u  DET3P(L)=AA9(L)**2#FXP(  -AAB  (L)  *  SORT  ( YB  (L  ) )  )/c% 

DO  35 C  1  =  1, NI 
DO  350  J= 1 , N J 
SUM=-. 

DO  34C  L  =  1 , NL 

34 J  SUM  =  SUMfFFOF0(U  *BETDP( L) ♦DYDX (L , I, J) 

35;  0BRX(I, J)=SUM 

---  IF  NC  ATTACK  CM  CFFENSr S  BYPASS  THIS  SECTION . -  — 

IF(NO*LC»0)  GC  TC  IjjO 
no  370  1  =  1, NI 
DO  3 7 C  K  =  1 , Nn 
SUM=C. 


363  SU"=SUIHFP0P9 <D  *BET9PI  t»  *!mO  ( L»  I»  K) 
378  OaeO(I»K>sSUM 
1000  RETURN 
ENC 
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subroutine  oerzb 


-  THESE  COMFUTICNS  ARE  THE  SAME  AS  THOSE  DESCRIBED  IN  RAC-CR-43  PART 

— -  IV-3.  WITH  EXCEPTIONS  AS  NOTED  IN  THE  COMMENT  CAROS  - - 

COMMON/SHAPE/  X<1009,  OEL(lu3>,  A(10O,  lui>  )  ,NV ,MC, MN,NP1  fNMl 
COMHON/CNE/NI,NJ,NII,NJJ,NM,NK,NL,NI1,NJ1,NI2,NJ2,NIIJ,  IEX,IEO,IEN 
I, IEUfIEXl,IEXENf I1,I2,I3,I4,I5,I6,I7,N0EFS,NC0FFS,NN,NU1,NU2,NU3 
2,I61,NC  ,NU4,NU5 ,NU6,NU7, NU8,N01,N0? 

COMHCN/THC/AA  (6)  ,AA*f  29  ,F«>0F<69  ,FPOPB<  21  ,DL<  5,79  ,ENL(lo  9,RB(11  9  , 
1£KB(1C9 ,TITLE2(1C9 , Ff 15  f7i  ,AB CIO ,15) , B(2,1G ) , 06 (2, 1C) , EH (15)  , 
2E(159,EKI15),TITLE1C159 ,PB< 1C , 15) , ANC15,7> ,SS (5 ,69 , LU(159 
CC*MCN/TFREE/FREFER ,PREFIN,PREINC , PRESET ,01(5,5)  ,R0<5> 

COMMON/FOUR/  C<15,7) , CB (10,15) ,ELNS<5,69 
C0MHCN/FIVE/XX<9, 159 ,0(9,79  ,EN (1C  9 , C(9 ,5) ,U(9,7) 

COMMON/SI X/Z (9,7  9  , ZB (2,109 ,F( 9,159 ,S( 10 9 ,FB (2,109 ,BETA (6) ,Y(6) , 
1BETAB(?9 , YB(29 ,BETB,BET 

COMMON/SEVEN/OFX <9,15,159 ,0F0(9,15, 79 , BETP (6 9 ,0YX (6 ,159 ,DBX<6,159 , 
10Y0(6,79 ,OBO(6,79 , BET BP (2 9 ,  BETPP(6> ,BETBPP<29 
CQKMCN/EIGHT/CSX (5*69 ,OSO(5,79 ,  OFiX<lU,5,69 ,CF2X( 1C ,5 , 69 , 

1DF10(10  ,5,79  ,DF20(13, 5,  79  ,  OYBX  <2,5,69  ,CYB0(2,5,79 

COMMON/NINE/DBBX (5,69 ,OBBO(5,79 
COMMON/ ELEVEN /OB XX (9a ,1 59 ,06X0 <9C ,79 ,0B0P C 30 ,59 
COMMON/THFLVE/DSXX (34  ,6  9,  OSXO(  10  ,79  ,0S00(6,5) 

COMMON/FCUPTN/OBBXX (3C, 3u ) ,300X0(30 ,35 9 , DPBDO (3C ,20 9 ,02YB (2) 

00  5  1=1, NI2 
TERM  =  SORT (Y  (1 9  9 
IF(T€RM9 A, ?, 4 

3  9FTPP(I9=-1.E41UC 
GO  TO  5 

4  9ETPP(I9=-AA(I)**3*£XP<-AA(I9*TERM  /  <4.*TER*9 

5  CONTINUE 

00  7C  I=NIi,NJI 
00  J-l , N J J 

NN2  s  (  J  -  1  9  *NI2  ♦  I  -  NI 

00  3C  KR=1,NJJ 

SUM=n. 

00  2C  K=1,NJ 
NN1= (K-19  *NI I+I 

2 J  SUM=SUMfOFXX (NN1 , J,KR>*  XX(I,K)  ♦  EK(K ) 

IF (NJ.EQ  *NJJ9  GO  TC  26 

00  25  K=NJ1,NJJ 

NN1  =  NIIJMK-1-NJ>  *NI2*i-NI 

25  SUM=SUMfDFXX<NNl,J,KR)*XX(I,K) *EK (K ) 

26  O2YXX=SUM  +  0FX(I,KR,J) *EK< KR9+DFX(I, J,KR9  *EK< J9 
I II=I-NT 

30  0BXX(NN2,KR9  =  FP0P(III 9* <BETPP< III9 *DYX  <IIIfJ)*DYX<III,KR9 ♦BETP (III 
19*02YXX) 

—  IF  NO  ATTACK  ON  OFFENSES  BYPASS  THIS  SFCTIGN  . 

IF  (NO*  LE • C  9  GO  TC  6? 

■*—  NO  REPLACES  NK  IN  EQUIVALENT  INSTRUCTION  OF  RAC-CP-43  PART  IV-B. 

00  50  KR  =  1 , NO 
SU*=l. 
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00  4C  K-i , NJ  , 

NN1=(K-1)*MI»X  , 

40  SUH=SUM*CFXO(NNl,JJ,KR>*XX<IfK)*EK(K) 

IF(NJ.FQ.NJJ)  GO  TO  46 
00  45  K=NJlfNJJ 
NNlsNIIJ^CK-l-NJ)*Nl2*I-NI 

45  SUK=SUH+OFXO(KNl,JfKR)*Xy (I,K)*EK(K) 

46  02YXD=SUrUDF0(I,  J,KR>*EMJ>  . 

K=I~NI 

5«  D9X0(NN2>KR>-FP0e(K)MBETPP(KIY0YX(K?J)»DY0(KfKR)*8ETP(KI*02YXOI  : 
60  CONTINUE  ; 

!  J 

—  IF  NO  ATTACK,  ON  DEFENSES  BYPASS’  THIS  SECTION - 

'•  \ 

IF(NC.UE.u)  gc  to  ?c 
00  56  K=1,N0 

i 

—  -  THESE  ARE  SECOND  PARTIAL  ,0  CROSS  D  DERIVATIVES  WHICH  OIC  NOT 

—  APPEAR  IN  RAC-CR-43  BUT  ARE  OESCRIBEO  IN  THE  PRESENT  REPORT  — - - 

NN?=  CI-1)  *  NC  *  K  -  NG  *  NI 

00  56  KR*1,N0  ,  f 

SU*=J. 

DO  54  J=i,NJj  1 

NNi  *  (J  -  i)  *  Nil  ♦  I: 

54  SUH=SUH  ♦  OFDC  (NNI  ,K>  KP )  *  XX(I,J)  *  EKCJ). 

IF(NJ.Ed.NJJ)  GO  TC  53  ■  * 

00  52  J=NJ1 fNJJ 

NNi  =  NIIJ  ♦  (J-i-NJ)  *  NI2  ♦  I  -  NI  i  , 

52  S|JH=SUM  f  DF0CCNN1  >  Kf  KP I  ♦  XX(I,J)  *  EK(J)  i 

53  02Y0D=SUM  \  i 

III-  I  -  MI 

56  DBDD(NN2,KR)  ~  FFOF  (I  mMBETPPJIII  )*DYO  (III  ,‘KJ  *CYO  (III  ,KR>  *BETP<  I 
III)  *D2YQn)  «* 

70  CONTINUE  ■  » 

RETURN  *  ' 

END 
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! SUBROUTINE  0ER2S 

—  THESE  COHPUTIONS  ARt  THE  SAME  AS  THOSE  DESCRIBEC  IN  RAC-CR-43  PART 

—  IV-0.  WITH  EXCEPTIONS  AS  NOTED  IN  THE  CCHNENT  CAROS  — - — 

COHNON/SNARE-'  XdCO),  OEL(13‘),  A (1.3 ,  It  t )  ,NV ,MC, FN ,N'P1  ,NP1 
C£W'MCN/CNE/NI,NJ,NII,NJJ,NM,NK>NL,NI1,NJ1,NI?,NJ2,NIIJ1IEX,IED,IEN 
l,IEU,IEXi,IEXEN,U,I2,I3,I4,I5,IF ,I7,N0EFS,NC0FFS  ,NN,NU1, NU2, NU3 
2,IB1,N0  »NU4,NU?  »NL6,  NU7,  NU8»NC1»ND2 
COMMCN/TNO/A A  (6)  ,AAB(2)  ,FF0P(6 ) ,FP0P8( 2)  ,OL(«,?>  ,ENl (10 ) ,RB(11 ) , 
lEKFdO  ,TITLE2(1C>  ,P(15,7)  ,  AB  (lu  ,15)  ,  B(  2 , 10  )  ,DE  (2,  1C)  ,  EH  (15)  , 
2£d5),EK(15),TITLEl(15)  ,P8(ll>,15)  ,AN(1  5,7)  ,SS (5,6) , LU(15> 
C0*'MCM/THFEE/FREFtR,PRFFIN,PREINC,PPEBET,ei{<3,5),R0(5> 

CO^NCN/FOUP/  C(15,7),Cr<(13,15>  ,ELNS(5,6) 
i  COMMON/FI VE/XX (9 ,15), 0(9, 7) ,EN (If) , C(9 ,5) ,U(9,7) 

C0MN0N/SIX/7t9,7)  ,Z0(2,lt)  ,F(9,15>  ,S(lli)  »FB(2»10)»  EETA  (6)  ,Y(6) « 
13ETA8(?),Y8(2)»0ETB,BFT 

COHrtON/SEVEN/CFX  <9,15 ,15) ,0F0( 9, 15,7) , BETP(6 > ,0YX (6  ,15) ,DBX(6,15> , 
10YO(6,7) ,080 <6, 7) , BETBP (2) »9ETPP(6) ,BETBPP(2) 

C0P80N/FICHT7CSX (5,6) ,0  SO  (5*7) ,  0F1X( 10 , 5 ,6) ,0F2X( 10 ,5 , 6) , 

lOFlOdC , 5,7)  ,CF2D(12,  5,  7)  ,  OYBX  (2,5,6)  ,0YB0(2,5,7) 

C0N*10N/NINF/0BnX  (5,6)  ,0880(5,7) 

C0*N0N/tLEVEN/0BXX  (90  ,15)  ,08X0  (9t  ,7)  ,OP0D  ( 30  ,5) 
O0*,N0N/TWELVE/0SXX(3ii  ,6 ) ,  OSXO(  30 ,7)  ,0S00(6,5) 
COMHON/FCURTN/OBEXX(30»  3u ) »OBBXO(33 ,35 ) ,08600(30 ,30 ) ,0 2 Y8 (2) 

00  4C  Hsl,NI 
00  4C  K=1,NJ 
MN?=(K-1)*NI*M 
00  4C  L-l ,NJ 
IF (S (H) -C  • ) 2C  »2T  ,25 
2‘;  DSXX (NN2, L) -2  . 

’  GO  TO  4J 
25  Sl)Y=t. 

Dd  32  J*1 »NJ 
NN1=(J-1)*NII+V 

3w  SUM=?UM+XX(M,J)*DFXX(NN1,K,L)*ELNS(M,J> 

OSXX()IN?,L)=OSX(M,t)*CSX(M,K)/S(M)+S(H)*(OFX(H,K,L) *ELNS(N,K)» 
10FX(F,L,K)*ELNS(F,L)»SUP)  /EMM) 

4C  CONTINUE 

— -  IF  NC  ATTACK  CN  CEFENSES  BYPASS  THIS  SECTION  - 

I:F(ND.LE.D  GO  TO  100 
00  9c  M=1,NI 
00  7C  K=1,NJ 

— -  NO  REPLACES  NK  IN  ECUIVALENT  INSTRUCTION  OF  RAC-CR-43  PART  IV-B. 

MN2* (K-l) *NI  +  F 
00  7u  L*1,N0 
IF(S(M)-C  .)45,45,5t 
45  OSXO  (NN2  ,L)  =  !*. 

GO, TO  73 
50  SUMsO. 

00  6C  J*1,NJ 
NN1*(J-1)*NII+M 

63  SUM=SUK*XX(M,J)*0FX0(NN1,K,L>*ELNS(M,J) 


b-28 


o  o  o  o 


OS>DCNN0,LMOSO(MfL)*OSXCMf  K)/S(H)fS<H)*CDFO  <MfK #LI *£LNS < M,KI  ♦ 
1$UM)/EN<M> 

70  CONTINUE 

—  THESE  ARE  SECCNO  PARTIAL  U  CROSS  U  DERIVATIVES  NHICH  010  NOT 

—  APPEAR  IN  RAC-CR-43  BUT  ARE  DESCRIBED  IN  THE  PRESENT  REPORT  — 

DO  85  K=1,ND 
NN2  -  <M-1)  *  NO  ♦  K 
00  8C  L=l|Nn 
IF(S(H|- 1.172,72, 75 
72  DSO0 (NN2 » L)  =  C. 

60  TO  80 
75  SUM  s  t  . 

HO  78  J=1 >NJ 
^Ni-(J-l)  *  Nil  ♦  P 

78  SUM  =  SUM  ♦  XX  (M ,  J )  *  OFOO  (NNi  ,K,L>  4  ELNS(M,J) 

3S09<NN2,L>  =  OSD(M,L)  *  DSO<M,K)  /  S(M)  ♦  SCM)  *  SUM  /  EN(M) 
8*  CONTINUE 
95  CONTINUE 
100  RETURN 
END 


SUEROUTINF  OFRgPp 


C 

C THESE  CQHPUTIONS  ARE  THE  SAME  AS  THOSE  OESCRIBFP  IN  RAC-CR-43  PART 

C  —  IV-3.  KITH  EXCEPTIONS  AS  NOTEO  IN  THE  COnfENT  CAROS  - 

c _ ; _ 

COMMON/SHARc/  VtlCu).  D€L(l)si>,  A  (1  lie  )  ,NV,HC,**N,NP1  ,N*»1 

_ CC)H0N/CN£/NI.NJ.NII.NJJ.NM.NK.NL.NI1.NJ1.NI2.NJ?.NIIJ.IEX.IE9.IEN 

1,ISU,IEX1,IEXEN,U,I2,I3,I4,I5,I6,I7,N0EFS,NC0FFS,NK,NU1,NU2,NU3 

2. 1 91. NO  .NU4.NUS.NU6.NU7.NU3.ND1.N02 _ 

COHMON/TWC/AA (6) ,AA8(2> ,FP0P(6> ,FP0PB(2> ,DL(9,7> ,ENK1G > . RB(13 » , 
1EKS (10 )  ,TITLt2(ir.)  ,P(15.7>  ,A3<  1C. L 5) .9(2. lu). 09 (2.10)  .EM (15)  . 
2E(15>,E<(15),TITIE1(15> ,PQ(lg,15> , AN(15, f ) ,SS(5,6) , IU(15) 

_ COMMON/ THREE /PREFER .prffI N, °REINC .PREOET ,91( 9.5) <RD  (5) _ 

COMMON/FOUR/  C(15,7),CB(19,15) ,ELNS(5,6) 

_ C0K*<0N/FivE/xx(q.io.r;t°. 7)  .en(  10.0(9.5 >  .u(9.7) _ 

COMMOM/SIX/Z (9,7 )  .7°  (2*16)  .F (9,15)  ,S  ( 1C  )  ,FP  (2, li. ) , BETA  (6 ) ,  Y(6 )  , 

ldFTAB(2) .Y9(?).3;T3,9ET _ 

CON1ON/SF.VEN/0FX  (9,15,15)  ,DCQ(9,15»7),  BETP (6  ),0YX(6,15)  ,0EX(6,15)  , 

19Yn(6,7),0e0(6.7).9ETBp(2).d£TPP(6)  ,9ET9PP(2) _ 

COMMON/EIGHT/CSX (5,6) ,0SD(5,7) ,  DF1X( If , 5 ,6) ,0F2X( 1C ,5 , 6) , 

lOFlOd"  ,5  ,7)  ,  CF?P  ( IT .  5,  7)  . _ 0Y3X(2,5.6  >  .CYBC(2,5.7) 

COMMON/NINF/OBPX  (5,6) ,0690(5,7) 

_ COHHCN/ELEYEN/CBXX (9 j ,15) ,00X0 (9- .7) ,0dPn(3Q ,5) _ 

COMMON/TWiiL  VE/DS  XX (3g ,6 > , OSXO( 30 ,7) ,DSCD<6,5> 

_ C0H90N/F0URTN/03FXX(3^.  30 .009X0(33, 35). OOPOCO;  .3:>,02YB(2) _ 

00  5  L  =  l.K'L 

_ TERM=SQRT (Y8(L) ) _ _ _ 

IF(TERM-C .>4,3,4 

3  iET9PP(L)=-l.E*l-i  _  _  _ 

GO  TO  5 

4  0ET'3PP  ( L  )  a  -AAiHL>»»3»EXP(-AA3(L)»TrRM)/(4.«TERH) _ 

5  CONTINUE 

_ 00  65  1=1, NI _ _ _ 

90  65  J=1,NJ 

_ 4N?=(J-1)»NI  FT _ _ _ _ _ _ 

00  6?  K=1,NI 

_ 00  3:  KR=1,NJ _ 

3N3  ■  (KR  -  1)  *  NI  ♦  K 

_ DO  21  L  =  1.NL _ _ 

SUH  =  v. . 

00  1"  K=1.NM _ _ _ 

)N1= (L“l ) *NM*M 

13  SUHgSUH»PFexX(NNl,NN?,NN3i»S(H)* E N(  M)  »  FK8  ( M) _ 

GO  TC(12,14) , L  ” 

12  TFRH1L=0F1X(K.I.  J)  _ _  _ _ _ _ _ 

TERH2L=0F1X(I,K, KR) 

GO  TO  15 _ _ _ _ _ 

14  TEFilL=0F2X(K,I, j) 

_ TERM2L=0F2X(I,K,KR) _ 

15  TERM?=Tcr«il*OSX (K,KR)*EN(K)*CKB(K) 

_  TERH5  =  TEPK2L*CSX  <I,J)*EN(I) »£Ke(I)  _ 

T£RH4=v . 

_  IF(I.EC.K)TEPH4=FB(L.I) *EN(I)*EKP(I)*OSXX (NN2, KR)  _ 

2,*  O2Y3a>=FFOP0a>*(SUM*T  tRKZ'*TERM2>'TEP*4)  “ 

_ SUM=t. _ } _ _ 

00  25  L=t,NL 

25  SUK=S»K»FFOFB(L) *(PETOFP(L) *0Y9 X ( l,I,J)»p YFX ( L ,K ,KR) +EETBP (L) » _ 


102VB (L) ) 

30  3RfiXX(NN2,NN3)=SUM _ 

IFCNO.LS*'.)  GC  TO  6? 

_ 3SL5LJ&=A*hn _ 

Ulc  =  (K-l)  ♦  NO  ♦  KR 

_ OP.-5*  l  =  _ 

$UH  =  % 

no  4i.  Msj  i  NM _ 

•JNt=(L-l>*NKf* 

40  SUM=SU*+0F5XCCKN1,KN2 «NN3 ) »S (M) »EN( M) ♦ EKE <N)  _ 

GO  TO (42 f 44)  ,1 

4?  TFfmL  =  QFlX(KfIt  J) _ _ _ 

TFPN2L=0F1C(I,K,KR) 

GO  TO  45 _ 

44  TFR*ll=0F2XCK,I, Jl 

TFRM2L=3F2H(I,KtKR) _ _ 

45  TER*2=TERN1L  *r)SP<  K,KR)  ♦ENC  i<)  ♦FKBOO 

TFRN3rTEP^2L  »CSX(I,J)*£N(I)»EKe W) _ 

TFRM4=:  . 

IF(I»EG«K)  T£RM4-sFB  CL«  I>  »£N  (  II  *  EKB  (I)  »D  SXO  <NN2 » KR ) 

50  02 Y 3 ( L ) =FF0PQCL) *(SUM+TFRM24TFRN?4TERM4) 

SUM  =  C  . _ _ _ 

OC  55  L-1,NL 

55  SUM=SUN»FPQPP(L) » ( RET 0PP( L) »OY 6X (L , I » J ) »0Y80 (L , K , KR) »9E TOP (L) 
132Y9CU) 

ec  055X0 (NN2tNN3)=SU* _ 

5?  CONTINUE 

55  CONTINUF _ _ 

i 

;  ---  IF  NC  ATTACK  CN  C£P£NSF 3  BYPASS  THIS  SFCTION . . 

IF (ND«LF • " )  GO  TO  1 JO 


, THESE  ARE  SECCNO  p ARTIAL  0  CROSS  0  _0  ERIVATI VFS  WHICH  QIC  NOT 

J  -  APPEAR  IN  RAC-CR-43 

bOR'  lTiTNI  . 

00  9i  Js 1 1 NO _ 

NN2  r’  (I- 1)  *  NO  ♦  J 

_ 00  9’.  K=1  ,NI _  _ _ 

00  9*  ~KR  =  1,NP 

IN3  =  ( K-l)  »  NO  »  KR _ _ 

00  8r  1=1 ,NL 

SU*»C  * _ 

00  7t  M=lfNM 

NNlsg-1)  »NM  +  M _ 

7r  SUV  =  SUH fGFPOn (NN 1 f NN2 |NN3 ) ♦S (H ) *EN( hV^tkb 

_ GC  TO  (72 » 74)  ,L _ _  _ 

72  TEFN1L=0F1C<K,I, J> 

TFPH21  =  DFID  Cl  t  K,  KR) _ _ _ 

GC  TO  75 

74  T£RilL*CF 20 ( K  »I t  J) _  _  _ 

TERM2L=0F20(I,K,KR) 

75  TFRN2=TERM1L  *OSO< K tKR) *EN< K I »EKR(K)  _  _ 

TFRH3«TPPH2L  •OS«MI|J)*eMCn*rKB(li 

TFRN4si. _ 

IF(I.FC*K)  TERM4=FB (t*  I)  *EN(I)*EKr(I)#nS0P(KN2fKf<) 

8?  0? YR (L) sFFQPH  (LI MSUM»TERM2»T£RH3+T£RM4) _ 


8*  5SM»SU«’FFCP8a)*t9P.T0PP<L>*nY9OCL,I,J)*OY0OtL,K,KR)*BEr^CO* 

in2Y3»L»  > 

gj  aaai9(NK?,KN3)=suM 
IJ..  PETURt* 

PNO 
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FUNCTION  OFXO(NNl,t,KS) 


—  THFSE  COMPUTIQNS  ARE  THE  SANE  AS  THOSE  OESCRIBEO  IN  RAC-CR-43  PART 

—  IV-8.  WITH  EXCEPTIONS  AS  NOTEO  IN  THE  COMMENT  CAROS  - 

COMMON/SHARE/  X(100),  DEL(130),  A(10o, IOC) ,NV,MC,MN,NP1 ,NM1 
C0MM0N/0NE/NI,NJ,NII,NJJ,NM,NK,NL?NI1,NJ1,NI2,NJ2,NIIJ,IEX,IE0, IEN 
1,IEU,IEX1,IEXEN, I1,I2,I3,I4,I5,I6,I7,NOEFS,NCOFFS,NN,NU1,NU2,NU3 
2,IB1,ND  ,NU4,NU5,NU6,NU7,NU8,ND1,N02 
C0MH0N/TW0/AAC6) ,AA9(2) ,FF 0P(6) ,FP0P8(2)  ,01(5,7) , ENL (10 ) , RP(10 ) , 
1EKB(1€),TITLE2(10) ,P(15,7>  ,AB(10 ,15) , B(2,1G ) ,06(2, 10) , EM (15)  , 

2E ( 15) ,EK (15) , TITLE1 (1 5) ,PP( 10,15) ,AN(15,7) ,SS(5,6 ) , IU(15) 
COMMON/THREE/PREFER, PREFIN, PREINC ,PRE8ET , 61(9,5) ,R0(5) 

COMMON/FOUR/  C (15 , 7) , CB (1C , 15) , FINS (5,6) 

COMMON/FIVE/XX (9 , 15) ,0(9,7) ,EN (1C ) , C(9 ,5) ,U(5,7) 

COMMON/SI X/Z (9,7 )  ,Z9(2,10) ,F(9,15) ,S ( IQ )  ,FB (2 , 10 ) , BETA (6) ,Y(6> , 
1BETA6(2) ,Y6( 2) ,8ET6 ,BET 

COMMON/ SEVEN/ CFX (9,15,15) ,OFO( 5,15, 7) , BFTP(6 ) ,OYX (6 ,15 ) ,06X(6,15) , 
1DY0(6,7) ,060(6,7 ),8ET9P(2), BETP°( 6) ,BETBPP(2) 

COMMON/EIGHT/OSX (5,6) ,0S0(5,7) ,  OF1X(10,5,6) ,0F2X( 10 ,5, 6) , 

IDF 10 (1C ,5,7) ,OF2D(10,5,7) ,  OYBX (2,5,6) ,QYB0(2,5,7) 

C0MM0N/NINE/C8BX (5,6) ,0BB0(5,7) 

COMMON/ ELF VEN/OB XX (9u ,15) ,OFXO (9C ,7),DBOO(30 ,5) 

COMMON/TWELVE/OSXX (30 ,6>, DSXH( 30 , 7) ,0S00(6,5) 
COMMON/FCURTN/nBcXX(30,33) ,CBBXD(3J ,35 > ,08B0C (3C ,20 ) ,02YB (21 

-  IF  NO  ATTACK  CN  DEFENSES  BYPASS  THIS  SECTION  - 

IF (NO. LE «C )  GC  TO  1?" 

IF(NNl-NIIJ) 8, 8, 7 
7  J=NJ  ♦  (NN1-NII J-l) /NI2  ♦  1 

I =NI  ♦  (NN1-NIIJ)  -  ( J- NJ-1 ) *N 12 
GO  TO  9 

3  J=(NN1-1) /NII+1 
I=NN1-(J-1)*NII 
9  SUH=0. 

00  7C  K=1,NK 

73  SUP=SUMfC(J,K)*0(I,K) *AN(L,K)/Z(I,K)#*2 

OFXOs-OFO(I, J,KSI*SUM-F(I , J ) *C ( J, KS ) *A N( L ,KS ) /Z ( I ,KS) 

103  RETURN 
END 


o  o  o  00 no 


FUNCTION  CFDP(NN1,L,KS> 

—  THESE  ARE  SECCNO  PARTIAL  0  CROSS  0  DERIVATIVES  WHICH  010  NOT 

—  APPEAR  IN  RAC-CR-43  8UT  ARE  DESCRIBED  IN  THE  PRESENT  REPORT  — - 

COMMON/ SHARE/  XIICG),  OEL(IQO),  A <103, 10P *  ,NV,MC,MN ,NP1  ,NH1 
C0MM0N/0NE/NI,NJ,NII,WJ,NM,NK,NI,NI1,NJ1,NI2,NJ2,NIIJ,IEX,IE0,IEN 
1,IEU,IEX1,IEXEN,I1,I2,I3,I4,I5,I6,I7,N0EFS,NC0FF$,NK,NU1,NU2,NU3 
2, 191, NO  ,NU4,NU5,NU6,NU7,NU8,N0i,N02 
COMMON/TWO/AA  (6) ,AAB(2> ,FP0P(6> ,FPOP8(2) ,OL( 9*7) ,ENK1Q ) , RB(10 » , 
lEKB(lC) ,TITLE2(1C) ,P(15,7>  ,AB(10 ,15) , B( 2 , 10 ) ,08 (2, 10) , EM <151 , 
2E(15)  ,EK(15)  ,TITLE1  (15)  ><>6(10,15)  ,AN(15,7)  ,SS(5,6>  ,UU<15) 
CCMMON/THPEE/PREFER, PREFIN, PREINC, PREBET, 81<9»5) ,RO(5) 

COMMON/FOUR/  C(15,7),CBC10,15) , FUNS (5,6) 

COMMON/FtV£/XX (9 ,15) ,0(9,7) ,EN (IS ) ,0(9 ,5 ) ,U(9,7) 

COMMON/ SIX/Z (9,7 )  ,78(2,10) ,F(9,15) ,S(10) ,F8 (2, 101 , BET A (6) ,Y(6 ) , 
1BETA9(2),Y8(2),BETB,RET 

COMMON/SEVEN/CFX (9,15,15) ,OFO(9,15, 7) ,BETP(6> ,0YX(6,15> ,08X(6,15) , 
1I)YO(6,7),C0O(6,7),BETBP(2)  ,0ETPP(6)  ,3ETBPP(2) 

COHMON/EIGHT/CSX (5,6) , OSD (5,7) ,  OF1X(10,5,6) ,0F2X( 10 ,5 ,6) , 

lDFlO (10 ,5 ,7) , OF2D ( 10, 5, 7) ,  OY0X (2,5,6)  ,0YB0(2,5,7) 

C0MM0N/NINE/0B6X (5,6) ,0000(5,7) 

COMMON/ ELEVEN/OBXX (90,1 5) ,08X0 (90 ,7 ) ,OBOO ( 30 ,5) 

COMMON/TWELVE/OSXX (30 ,6) , OSXO( 30 ,7) ,OSOO<6,5) 

COMMON/ F0URTN/0BPXX(3v,,  33)  ,088X0(30  ,35) ,  08000  (30 ,30 )  ,02YB(2) 

— -  IF  NO  ATTACK  ON  OEFENSES  BYPASS  THIS  SECTION  - 

IF(NO.LE.C)  GC  TO  IOC 
IF(NN1-NIIJ)3,8,7 

7  J=N J  ♦  (KN1-NIIJ-D/NI2  +  1 

I  =  NI  ♦  (NN1-MIJ)  -  (J-NJ-1)*NI2 
GO  TO  9 

8  Js(NNl-l)/NII*l 
I=NN1-(J-1)*NII 

9  OFDO  =  F(I,J)4C(J,L)*C(J,KS)  /  (Z( I,L> *Z( I,KS )) 

100  RETURN 

END 


o  o  o  o 


FUNCTION  DF8XX(NN1,NN2,NN3) 


—  THESE  COMPUTIONS  ARE  THE  SAME  AS  THOSE  DESCRIPEO  IN  RAC-CR-43  PART 

—  IV-R.  WITH  EXCEPTIONS  AS  NOTEO  IN  THE  COMMENT  CAROS  — — 

COMMON/SHARE/  X(1GU),  OEKiOC),  A(100, 110) ,NV ,MC ,MN ,NP1  ,NM1 
COMNON/ONE/NI,NJ,NII,NJJ, NM,NK,NL,NI1,NJ1 ,NI2,NJ2,NIIJ, IEX,IEO,IEN 
1, IEUtlEXl, IEXEN, 11,12,13, I4,I5,I6,I7,N0EFS,NCGFFS,NN,NU1,NU2,NU3 
2, 181, NO  ,NU4,NU5,NU6,NU7,NUQ,NOi,NO? 

COMMON/TWC/AA (6) ,AAB(2)  ,FP0P<6) ,FP OPB(2) ,OL ( 5 ,7 ) ,ENL <10 ) , R0C1O ) , 
lEKB(U)  ,TITLE2U0)  ,P(15,7)  ,AB(1C  ,151 ,  B<  2 ,10  )  ,  OB  (2,  lu)  ,  EM  (15)  , 

2E  ( 15)  ,EK(  151  ,  TITLEM15)  ,PB(U.,15)  ,AN(15,7)  ,SSf 5 ,6) ,IUC 15) 
COMMON/THREE /PREFE 8, P RE FI N,PRE INC, PREbET ,81(9,5) ,RO(5) 

COMMON/FOUR/  C (1 5, 7) , C8 t 10 , 15) ,EINS(5,6> 

COMMON/FI VE/XX(9, 15), OC 9, 7) ,EN(10) , Q<9,5> ,U(9,7) 

COMMCN/SIX/Z  f  9,7 )  ,28(2,10) ,F<9,15) ,S(  10  )  ,F9  (2,10  ,  EETA  (6),  YC6)  , 
19ETAB (2) , Y6 (2 ) ,BF  T8,BET 

COMMON/SEVEN/CFX (9,15,15)  ,OFO(  9 , 15,7)  ,  BETP  (6)  ,OYX  (6  ,15) ,0BX(6,15) , 
10Y0(6,7) ,0B0(6,7) , BET BP (2) , 9£TPP(6)  ,BETBPP(2) 

COMMON/EIGHT/CSX (5,6) ,OSO(5,7) ,  DF1X( 10,5,6)  ,0F2X( 10 ,5 , 6) , 

19F10(10 ,5,7) ,OF2P(10,5,7) ,  OYBX (2,5,6) ,OYBO(2,5,7) 

COMMON/NINE/OPBX (5,6) ,0880(5,7) 

CO MM ON/ EL EVEN /OB XX (9u ,15) ,06X0 (9C ,7) ,0800(30 ,5) 

COMMON/TWELVE /OS XX (30 ,6 ) , 0SXO(30 ,7) ,OSOO (6 ,5) 

C0MM0N/F0URTN/0B8XX ( 30,  33 ), 088X0(30,35 ), 08800(30, 30), 02YB (2) 

L* (NN1-1) /NM+1 

M=NN1-(L-1)»NM 

J=(NN2-1)/NI*1 

I-NN2-(J-i)*NI 

KR- (NN3-1) /NI*1 

K=NN3- (KR-1) *NI 

SUMl=r. 

SUM2=C* 

00  10  N*1,NN 

TERM*Cfi (M,N)*0B(L,N)*A8(I,N) 

SUM1=SUM1*TERM/ZR(L,N)**2 
10  SUM2-SUM2+TERM*AB(K,N)/ZB (L ,N) 

TEPML=0F1X(M,K,KR) 

IF(L*E0.2)TERML=DF2X(M,K,XR) 

TEPM1=-TFRML*0SX (I , J) *EN( I) *SUM1 
TERM2«C • 

IF(I.E0.K)TERM2r-F0(L,M)*DSXX(NN2,KR)*EN(I)*$UMl 
TERM3ss2»*FP0PB(L)#fcN(  I)  ♦EN  (K)  *OSX  (I  ,J)  *PSX(K,KR)*F8  (l,W)#$UN2 
0FBXX*FPCFB(t)MTERMl*TERM2*TERM3) 

RETURN 

ENO 
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FUNCTION  0F0X0(NN1,NN2, NN3) 


—  THESE  CCMPUTICNS  ARE  THF  SANE  as  THOSE  DESCRIBED  IN  RAC-CR-43  PART 

—  iv-R.  KITH  EXCEPTIONS  AS  NOTEO  IN  THE  COHMENT  CAROS  - 

CCMKCN/SHARE/  X(100>,  OEL(IjC),  A (100, 100 ) ,NV ,MC ,NN ,NP1 ,NN1 
COMHON/ONf /NI ,NJ,NII,NJJ,NM,NK,NL»NI1,NJ1,NI2 ,NJ?,NII J, I£X,IEO , IEN 
1,IEU,IEX1,IEXEN, 11,12, 13, 14 ,T5 , IF, I7,NDEFS ,NC0FFS ,NK ,NU1, NU2,NU3 
2,  I'd, NO  , NU4 ,NU5 ,NU6,NU7,NU9,N01,N02 
CONNON/THC/AA (6) ,AA9(2) ,FP0P(6) ,FPOP0(2) ,0L(9,7> ,ENL (10 ) , RB(10 1 , 
lFKB(10),TITLE2(lo) ,P(15,7)  ,AR(1C ,15) , B( 2 , 1U ) ,DB (2, 10) , EH (15) , 

2£< 15) ,FK( 15) ,TITLE1(15) ,P6(10,15) , AN(15,7 ) ,SS (5 ,6) , LU<15) 
COHNON/THRtE/FEEFFR  .PfiFFIN,  PRE  INC,  PREBET  ,B1(9,5),R0(5) 

C0M50N/F0UR/  C(1 5, 7) , C8 (10 , 15) ,ELNS(5,6> 

COHMON/FI VE/XX (9 ,15) , D( 9, 7) ,EN ( 10 ) , 0(9 ,5 ) ,U( 9, 7) 

C0KMCN/SIX/Z(9,7)  , 79  (2 ,1C ) ,F ( 9 , 15) ,S(  10  >  , FB  (2, 10) , eET A  (6) , Y(6) , 
1BETAB(2) ,YB(2),BET9,BET 

COMMONXSE VEN/BFX (9,15,15) ,0FD ( 9 , 15, 7) , «ETP (6 ) ,OYX (6 ,15) ,0eX(6,15> , 
10YC (6 ,7 ) ,0 B0 (6, 7), RET  BP (2), 9ETPP(6) ,8ETBPP(2) 

COMNON/EIGHT/CSX (5,6) ,0S0(5,7) ,  OF1X(10,5,6) ,0F2X( 10 ,5, 6) , 

lOFlQdC1 ,5  ,7) ,  0F20(  10 , 5,  7)  ,  0YBX(2,5,6)  ,OYe0(2,5,7) 

COMHCN/NINE/OBPX (5,6) ,0BB0(5,7) 

CCMMCN/EL£VEN/CBXX(9Q ,15) ,OeXO (90 ,7) ,0  BOD (30 ,5) 
C0H«)0N7THELVE/nSXX<3B  ,6 ) ,  DSXO(  30 , 7)  , OS  00(6, 5) 

COMMON/FOURTN/OBBXX (3ii»33)»OBBXD(3Q,35)»OBBOO(3U  ,30),02YB(2) 

—  IF  NO  ATTACK  ON  OEFENSFS  BYPASS  THIS  SECTION  - 

IF(NC.LE.o)  GO  TO  100 
La (NN1-1) /NH  +  1 
M=NN1-(L-1)*NN 
J= (NN2-1) 7NI+1 
I=NN2- ( J-l ) *NI 

—  NO  REPLACES  NK  TN  EQUIVALENT  INSTRUCTION  OF  RAC-CR-43  PART  IV-B. 

K=  (NN3  -  1)  /  NO  ♦  1 

...  NO  REPLACES  NK  IN  EQUIVALENT  INSTRUCTION  CF  RAC-CP-43  PART  IV-0. 

KR  =  NN3  -  (K-l)  *  NO 
SUMt*3. 

SUM2=C  . 

00  50  N*1,NN 

TERN*C9(M,N)*0E(L,N)*AP(I,N) 

SUH1«SUM1*TERH/ZB(L,N)**2 
3J  SUN2*SUH2+TERN*AB(K,N)/ZB(L,N) *•? 

TERNL*0F10(M,K,KR) 

IF(L.E0.2)TFPI*L*OF20(N,K,KR) 

TERN1-TFRHL*0SX (I,J)*EN(I ) *SUM1 
TERH2*C. 

IF(I.EO.K)T£RH2*FB(L,H) *0SX0(NN2,KR)*EMI)*SUH1 
TERM3a-2.*FB(L,M)aEN(I) *EN (K) * OSX (I , J) *0S0 (K ,KR) *FPOPB ( L) *SUM2 
OFOXC*-FPCPB<L)*(TERHl*TERH2»TERM3) 

ICO  RETURN 
ENO 


O  l) 


FUNCTION  DF800(NN1,NN2,NN3> 


:  —  THESE  ARE  SFCCNO  PART IAt  0  CROSS  0  DERIVATIVES  NHICH  DIC  NOT 
:  —  APPEAR  IN  RAC-CR-43  BUT  ARE  OFCRIBEO  IN  THE  PRESENT  REPOOT  - 

COMMON/SHARE/  XfluO),  DEL  (ICC ) >  A  (lCu ,  10  G )  ,N  V  ,HC ,  PN  ,NP1  ,NM1 
C0M10N/0NE7NI,NJ,NII,NJJ,NM,NK,NI,NI1,NJ1,NI2,NJ2,NIIJ,IEX,IE1,IEN 
1,IFU,IEX1,I£XEN,I1,I2,I3,I<*,I5,I6,I7,NOEFS,NCOFFS,NN,NU1,NU2,NU3 
2, 181, NO  « NU4 ,NU5 ,NU6,NU7,NU8,N01,N02 
COPMON/TWO/A A (6) ,AA9(2) ,FP0P(6> ,FPOPB(2> ,CL(9,7> ,ENL<1« ) , RB(10 > , 
lEKB(lO) , T ITL£2< ID ,P(15,7)  , AB (1C ,151 , 8(2 ,10 » ,DE <2, 10) ,EH(15> , 
2E(15),EK(15> ,TITLE1(15> ,PP(lw,15) ,AN(15,7) ,SS(5,6),UU(15» 
COPPON/ThREE/PRFFER*PREFIN, PREINC , PRESET ,61(9,5) ,R0(5) 

C0MM0N/F0UR7  C(15 , 7) ,  CB  (lu,  15)  , FLNS  (5,6) 

COKMCN/FIVE/XX (9 ,15) ,0(9,7) ,EN (1C ) , C(9 ,5 ) ,U( 9, 7) 

COHMON/SIX/Z (9,7)  ,Z8(2,10) ,F(9,15) ,S( 19 ) , FB (2, 10) , EETA (6 ) , YC6 ) , 
19ETAB(2) » Y8(2) ,BFTE  »BET 

COPMON/SFVFN/CFX (9,15,15) ,OFO (9 ,15, 7) , PETP (6 > , OYX (6 ,15) ,0ex(6,15) , 
10Y0(6»7), CBO (6,7 ) , BET  BP (2 ) , BETPP(  6) ,BETBPP(2) 

COHNON/EIGHT/OSX (5,6) ,OSO(5,7) ,  0F1X(10,5,6)  ,0F2X( 10 ,5, 6) , 

10F 10(1C , 5,7)  ,0F20 ( 10 , 5, 7) ,  0YBX(2,5,6) ,CYBC(2,5,7) 

COMMON/NINE/CBBX (5,6) ,0080(5,7) 

COPHON/ELEVEN/OBXX (90 ,15) ,08X0(90 ,7),08DO(30 ,5) 
C0“N0N/THFLVE/nSXX(3l,  ,6) ,  OSXO(  30 ,7)  ,0S00(6,5) 
CONMON/F0URTN/C9BXX(30,3u ) ,C9BXD(39 ,35) , 08E0C ( 3C ,30 > ,02 YB (2) 

- IF  NO  ATTACK  ON  DEFENSES  BYFASS  THIS  SECTION - - 

IF(NC.LE.C)  GO  TO  1JC 
Ls(NNl-l)/N**l 
M=NN1-(L-1)*NM 
I  *  (  NN2  -  1)  /  NC  ♦  1 
J  *  NN2  -  (  I  -  1  >  *  NO 
K=  (NN3  -  1)  /  NP  ♦  1 
<R  =  NN3  -  (K-l)  *  NO 
SUH1=7. 

SUP2=0. 

90  3 j  N*1 ,NN 

TFRHsCB(N,N)*08(L,N)*A8(I,N) 

SUM1*SUH1+TERH/ZP(L,N)**2 
30  SU-2*SUM2»TERP*AE(K,N)/ZB(L,N) **3 
TERML*0F1P(M,K,KP) 

IF  (L.EC.2)TERMl=OF2D(H*K,KR) 

TERH1«TERML*0S0(I,J)*EN(I)*SUM1 
TFR*2*C. 

IF(I.F0.K)TERH2=FB(L,M) *OSOO(NN2,KR)*EN'T)*SUM1 
TERM3*-2.*F0(L,H)*EN(I) *EN(K)*OSO(I,J) *DSO(K,KR) *FPOPB(L>  *SUH2 
OFROO^-FPOPBtDMTERHl+TERMZ+TFRPZ) 

100  RETURN 
FNH 
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SUBPOUIME.  T  ap  t  p 

S^^UjB,egyT JE^jrs.  a  SWUNG  E.EA.SI3LIL»QI»^ 

o 

. . MWM05IYSMAPE/  JC(ltO),  DEI  (10  G) ,  A  ( it)  tf  ,104)  »  N  V.#MC.»..H!L»BPl.*.M.Hi . . 

COW#W/'WT/*>Jr,  N'J'f  WT  I  fMJJ,  NM  ,  NK  t  NL  ,  Ni  l  ♦  N  J1  ,  N I  2  ,  N  J?  ,  N 1 1 J,  I FX  ,  I  ED,  I EN 
.JL^IEUaEXUIEK£NfIi,I2.n*I4,I5fI6>I7fNDEFS,.ND0F.E^NN^NUl*NU2^1U3 
2,TRifM0  ,NU4f NU5fNU6,NU7,NUP,N01,Nn2 

_ ,  A  A  31 2)  .t F.90 *  (6  >  , FPOP  B ( 2 ) , DL  ( 9 ,T)  *£NL.  CIO  )-*PJB.lli3)^ _ 

1FKO (10) ,TITLF?(1C>,P(15,7)  , AB(iw,15) ,B (2, 1 0 ) ,0P(2, 1 3) ,  EM (15) , 

2E_(.L5)  *  EK.il  5.)  , TITLE! (15)  ,PB( 1C , 15) , AN (15,7)  ,SS (5,6) ,  UU(15> . . 

COMVON/TMOEE/PPE^rP^oPEETN, PPEINC,oPEBET,31  (9,5)  ,Pf)  (5) 

.  C  Q.HOTN/  EQUIPS  •.  C  ( 15  ,  7 )  ,  CB  ( 1 G ,  1 5 ) ,  ELNS  (5,6)  . . . . . . . 

COMMON /F TV F7 X X (9,15) ,0(9,7) ,FN(1C) ,0  (9,5) ,U (9,7) 

-M**Q*mX/7  ( 9,_7)  ,79(2,11)  ,F  (9,15)  ,5  ( 1 0)  ,  F  B  ( 2 , 1 C )  ,  BETA  (  6)  ,  Y  (6)  , . 

19ET  AP( 2) , YB(  2 ) , BET9 , 9ET 

. C0NNON/SEVEN/0FX  (°,  1^,15)  ,0FD(9 ,15, 7) , 9ETD  (6)  ,.CYX (6 ,153 ,03X16,1.5  )  # 

10YO (6,?), 000(6,7) ,9ETPP(2),PFTD°(6) ,9ETB99(2) 

. .  .  C0^WQN/EIGKt/aSX(5,6)  ,DS0(5,7)  ,  OF1X  (10 ,5, 6)  ,0F2X  IU+5+6J  , 

10F10(1Q,5,7)  , OF 20 (10, 5, 7) ,  OY9X (2,5,6) ,OYBO (2,5,7) 

_ _ C!)MpQj|/niHE/llSlBXCrf6)  ,D9?C(5,7)  . . . . . . 

COMMON  /OUTIN/  P9P(?" ),7N1(2‘) , 7N2(2G) ,  RAT (1C) 

..  . COMMON  /EQAL/  H,H1,M7  . . 

DTMEN9T0N  ALL (  lb) 
flY  =.  PP99ET/ (PPEFFP  ♦  ORFR^T) 

>yviv’-~cv“ 

_ -S3A  LL=  •  31 . _ _  _ _ 

FL1  =  FLOaT(NTI)/rLnflT{NlT4-l) 

. FL_?  r  fLOAT(MT2)  /  FtPAT (NI?  +  1) 

po  l»j  J s  1 , N  J 

. 11  ALL  ( J)  =  ELI  *  EM  (J)  . .  ... 

00  ?0  TsK'IlfMH 

_  30_  10.  J  =  N;J  ij  NJJ  .  . . 

2U  XX(I.j)  =  FM(j)  ♦  FL?  *  Fdod(I-NI)  +  SMALL 

TF  NO  ATTACK  ON  OFFENSES  BYPASS  THTS  SECTION  —  —  —  — 


T F ( NO ,L F • 0 )  00  To  SGE 

_ J E1H02 1 E,0.t.2JL  GO.  TP  ICO  . 

00  to  (7.1,40,50,60  ,MPl 

_ 3*  SUM  =  1.1 

QO  7*  KslfN0 

. .  ALUK)  =  ALt(K)  /?,0 

no  7?  1=1, NT 

_ J.2..  UlI.iKl.  =..  CE_T__._  aiL(K).  /  FLOAT(.NT)  ♦  SMALL 

no  74  T  =  NI  1 ,  N T I 

. ?4  JHIfK)  =.  rv  *  ALL  (K>  ♦  FPno(j.Nli  ♦  SMALL 

79  CONTINUE 
„  .  CO  .JO  503 

43  ALL  (1)  *  ALL ( 1)  /  7.0 

_ A ILL  ( ?J„ .*  ALL  (  2).  „ /.  2 .  *J.  . 

ALL (7)  =  ALL ( 7)  /  2.' 

. . 00_4F_,J  =  lf  ? 

00  4?  T  =  1 ,  NT 

. .  bZ  U(IjJ)  -  .  ALL  ( 1 )  ♦  Cr  /  FLOAT(NI)  +  sm\ll 

00  44  I =NI 1 , N IT 

*  cv  *  fpop(i-ni) 
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46  rOMTIMM- 

DO  47  1=1, NT 


U(T,?> 

=  all  (?) 

* 

n  F 

/ 

FLOAT(MT) 

♦  small 

47 

U(T,4) 

=  ALL  (3) 

* 

rF 

/ 

FLOAT (NT) 

♦  SMALL 

po  4* 

I -NT  If  MIX 

U(I,S) 

=  .ALL(?) 

♦ 

rv 

♦ 

FPOP(I-NT) 

♦  SMALL 

4« 

U  ( T  4) 

=  ALL(^) 

* 

rv 

FPOP(I-NT) 

♦  SMALL 

GO  Tn 

**  j; 

5J 

ALL (1) 

=  ALL(l)/ 

4 

.? 

ALL  (?) 

=  ALL(?)/ 

2 

on  gg 

4=1.^ 

no  c? 

T  =  1  ♦  K,T 

G? 

U  (  T  ,  J  ) 

=  ALL(l) 

* 

CF 

/ 

FLnAT (NT ) 

♦  small 

DO  54 

I=NI1,mii 

04 

U(I,J) 

=  ®  LL  ( 1) 

* 

nv 

» 

CO OP (I- MX) 

♦SMALL 

56 

CONTINUE 

00  G7 

1=1,  NT 

57 

•J(I,4) 

=  ALL (2) 

♦ 

CF 

/ 

cLnAT  (NI ) 

♦  SMALL 

DO  G« 

T  =NT 1, NTT 

G  ^ 

U  (  I  ,  4) 

=  ALL ( 2) 

* 

rv 

♦ 

ropp ( T-NI) 

♦  SMALL 

Gn  TO  go: 

6'j.  AUL  f  1>  =  Al L  <  1 )  / 

DO  66  J=1,ND 
DO  6?  T=t,MT 

6?  U  (T ,  J)  =  flLL(l)  *  oc  /  FLOAT(KT)  ♦  SMALL 
nn  64  I  =L!1 1 ,  M  1 1 

54  U(T,J)=  ALL ( 1 )  *  rV  *  FonP(I-MX)  +  SWALL 
6*5  CONTINUE 
GO  To  G * 

20*3  ALL  <1>  =  ALL  (  1)  /  ^.C 
ALL  (?)  =  ALL  f  ,?>  /  7.  : 

K  =  j 

DO  ?19  J=l,4 
tp(j.gt.?)  <<  =  ? 

nn  202  T  =  t ,NT 

2)2  U(T,J)  r  ALL(K)  *  OF  /  CL  OA  T  (NT )  +  SMALL 

nn  204  T=MT1 , NT 2 

2y 4  U ( T  9  J )  =  ALL(K>  *  CV  *  F°f  ° ( J-M I)  +  SMALL 
21 O  CONTINUE 
*>;. :  on  .  J=  i*MJ 
no  i=  itMv 

5?s  xxa,j>  =  op  *  alkj)  /  floatcni)  ♦  small 

no  c-r,  T=NT1,NTT 

VX(TfJ)=  CV  *  ALL  ( J )  *  F°0°(T-MT)  *  SMALL 

6  j  0  OONTTMIjF 

on  7  0  •  T =1  j  NIT 
no  721  j=i f n j 

Ms ( J-l ) *MT I  * 

7)1  V(N)=  XX  (T , J) 

On  7l«J  J=1,NJ2 
On  71]  T=1,MI2 
M=  N.TTJ  ♦  (4~  1 ) *  MT  ?  ♦  T 
71;  Y  (N) sXX (NT+T  «  NJ  ♦  J ) 

-  JF  NO  ATTACK  ON  GFFFMcrs  ?Y°ASS  THIS  SECTION 

C 

TF(Nn.LF.O.G)  On  to  l c J 3 
HO  7.23..  1=1  ,NTT 


no  7 21  K=ltN* 

K'=  Try  +  fT-l)*h»n  +  if 
7?»;  XCM|=  U(T,K) 

3  P.5TUP*I 
ENO 


<5*S> 
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SUBROUTINE  REAOOU  (KT) 


i 


147  5<J 


X 


—  THIS  SUBROUTINE  SUMMARIZES  THE  RESULTS  OF ' AN  EXCHANGE  — - 

•  ! 

COMMON/SHARE/  X(ICC),  OFL(l)C),  A  (ICC  ,  l9ty  ,NV,MC  ,MN  ,NP1  ,NN1 
COMMON/ ONE/NI,NJ  ,  NII ,  NJ  J,  NH  ,NK, NL  »N II  ,  NJ1 ,  \I  2, N.J2,N II  J,  IE,X,  IEO,  itN 
li IEU,IEX1,IEXEN, 11,12,13, 14 ,15 , 16 ,I7,N0EFS,NC0FFS ,NN ,Nlfl,NU2,NU3 
2, 101, NO  ,NU4,NU5,NU6:»NU7,NU8,N01,N02 
COMMON/ THO/ A A (6)  ,AAB(2) ,FF0P(6 ) ,FF0F8(  2) ,0L(B,7I ,ENL(1G) , RB(19 ) , 
1EK8(10),TITLE2(11) ,P<15,7)  , AB (10 ,151 , BC2 , 10 > ,OE C2, 10) ,EM<15)  ,  1 

2E(15),EK(15)  ,  TITLEK15)  ,°B(  15  ,15)  ,AN(15,7) ,SS (5 ,6) , LU(15) 
C0MM0N/THPEF/FREFER,PRFFIN,FREINC,PREBFT,81(9,5),R0(5) 

COHMON/FOUR/  C(15,7),tB(10, 15) ,  FLNS  (5, 6)  , 

COMMON/FI VE/XX(9» 15)  ,0(9,7)  ,FN(10)  ,0(9,5)  ,U(9,'7> 

COMMON/SIX/Z  (9,7)  ,Z8(2,1«)  ,F( 9,15)  ,S(  1C  >  ,FB  (2, 1C  f  ,  EETA '(6 >  ,Y(6 )  ,  , 
19ETA0(?)  ,YP(.2),8FTP,BET  1 

COMMON /SEVEN/CFX (9,15,15) ,0F0(9,15, 7) , DETPC6) , OVX (6 ,15) iDBX(6,15> , 
10Y0(S, ’) , 080 (6,7) , 8ETBP (2 ) , PETPP(6) ,BETPPP(2 ) 

COMMON/EIGHT/CSX  (5,5)  ,0S0(5,7)  ,  nFlX(13,S,6>  ,0R2X(  10,'5, 6) , 

10F10(1G,5,7) ,0P2O(13,?,7) ,  1 0 VBX (2, 5 ,6)  ,0YB0( 2,5, 7) 

COMMON/NINF/OBBX (5,6) ,OBd0(5,7)  1 

COMMON  /OUTIN/  RBB(2o),ZNl(2u) ,ZN2(2c) ,RAT(lt) 

DIMENSION  Xl(2u) ,X2(2C) ,X3(2') ,X4(20),X5(20»,X6(2C) 

WRITE  ( 6 ,655) 

WRITE  (6,669)  '  i 

CALL  IOENTV  , 

CALL  FRACTS 

95j  FORMAT(  //1CX,  ?ch*****  9AMAGF,  TO  SIDE  ONE  ♦,**** ,F5 .3)  , 

955  F0RMAT(  //IjX,  3I.H*****  OAM AGE  TO  SIDE  TWC  **»**, F5. 3) 

WRITE  (6,495) 

WRITE  (6,95.0  )  BETP  '  i 

WRITE(6 ,955) PET 

WRITE  (6,41i)  1 

WRITE  (6,493)  ,i  ‘ 

WRITE(6,82E) (TITLEl(J), J=1,NJ> 

00  955  1=1, NI  ’  , 

00  840  J=1,NJ 

X1(J)  =  XV(I,J)  /  R83*J>  1  '  , 

X2(J)  =  XX(I,J) 

XT ( J)  =  X2(J)  *  7N1(J)  .  1 

X4(J)  =  X3(J)  *  F( I , J ) 

X5  ( J)  =  X4(J)  /  EN (I) 

IF(X5(J)  .GT.14C.)  X5('J)  =14C  • 

84J  X6(J)  =  SS(I,‘J)  **  X5  (J ) 

WRITE(6»8C5) (X1(J) ,J=1,NJ)  ■ 

WRITE(6,8i:) (X2( J) ,J=1,NJ)  , 

WRITE(6,815ITITLF2(I) ,(X3(J),J=1,NJ) • 

WRITE  (6,8.K)  (X4(  J)  ,  J=  1,  NJ) 

WRTTE(6,82i.)  (X5(J)  ,J=1,NJ) 

WRITE (6, 8?.) (X6( J) ,J=t,NJ) 

855  CONTINUE  ' 

WRITE (6 ,655) . 

WRITE  (6,54?)  1 

IF(NJJ,LE.12)  Ml  =  NJJ  1 

IF(NJJ.GT.12j  Ml  *  12  i  * 

WRITE(6,8JC)  (TITLEKJ),  J=1,M1) 

00  896  1=1, NI2  S  ,  1  ■ 
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I 


154C0 

i 

154C5 
,  15905 


16195 


TI  =  NI  ♦  I 

DC  985  J=1,M1 

X1CJI  =  XXCII.J)  /  R8B(J) 

X2CJJ  a  XXCII,J) 

X3(J)  a  XXUI,J)*FCII,J) 

885  X4CJ)  a  X3C  J)  *CKtJ) 

HRIim.mi  CX1(J)  ,Jal,Hl) 
WRITE<6,895)  I,(X2(J) 
WRITE<6,810) CX3C  J)  fJaiyr<i) 

WRITE (6, 8 li)  (X4C  J)  »J=1,M1) 

896  CONTINUE 

IF(NJJ.LE.12)  GO  TC  905 
WRITFC6 , 925) C TITLF1 CJ ) , J=13 ,NJ J) 
00  9' t  1=1, NI2 
00  989  J=13,NJJ 
II  =  NI  ♦  I 

X1(J)  a  XX(II,J)  /  R80(J) 

.  X2(J)  a  XX ( II , J) 

X?(J)  a  XX(II,J)*F(II,J) 

889  X4CJ)  i  X3(J)  *EK( J) 

WRITE(6,83C) (X1(J) ,Jal7,NJJ) 
WRITE (6 , 895)  I, < X2 l J) , J=13, NJ J) 

'  WRITE(6,81C) ( X3( J) ,  J=  13 ,N JJ) 
WRITEC6',  810)  <X4(J)  fJ=13,NJJ) 

900  CONTINUE 
9C5  CONTINUE 

WRITE (6 ,655) 


— *  IF  NC  ATTACK  CN  CEFENSrS  BY»ASS  THIS  SECTION 


IMNC.LE.O)  GC  TC  787 


1  WRITE (6, 575) 

WRITE (6 , 835)  (K,K=1,NQ) 

00  781  1=1, Nil 

DO  77f  M=l,NO 

XHN)  =  U  ( I  ,N)  /  RBB  ( M) 

X3(M)  a  01  (I,N)/(?0(W) 

X4  M)  a  C(I,M)/RQ(M) 

77;  X2 (4)  =  U(I,»j 

, WRITE  (6, 830  (X1(M> ,M=1,ND) 

WR ITE C6 ,895)  I,  (  X2CM)  ,M=i  ,ND) 
WRITE  <6 , 91u ) ( X  3 ( N ) , Ha  1 , NO  I 
W9ITE(6,81i.)  SX4CH)  ,M=l,NO) 

780  CONTINUE 
787  CONTINUE 

WRITE (6,655) 

WRITE  (6,545) 


00  740  L=1,NL 
00  750  4=1, NM 
X1(M>  a  EMM) 

X2(4)  a  EMN)*  S  CM)  /  PP(M) 
>3(4)  a  X? (M )  *  R«MH) 

XMM)  =  X3(N)  *  FB(L,M) 

750  X5H)  a  X  4  (  4  )  *  FKB(M) 

WRITE(6,825) CTITLE2(M),M=1,nm) 
WRITE(6,83C) CXl(H) , Me  1, NM) 
WRITE  <6 ,895)  L, ( X2 (M) ,m=1 ,NM) 
WRITE<6,810 (X3CM) ,M=1,N4) 
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-n»^cet^  *■ r?****. 


WRITE (6, 8101 <X4<H| ,M=i,NM) 

WRITE(6,81U) (X5(M) ,M=1,NM) 

760  CONTINUE 

NU8  =  NU8  ♦  1 

IMNU8.GT.NU4)  KT  =  2 

IF (NU8 . GT .NU4)  RETURN 

IF (RAT (NU8  ) .GT. l.E+7)  PREBET=10C. 

IF  (RAT (NU8  ) .GT, l.E+7)  PREFER=0. 

IF  (RAT  (NUB  >.GT.l.F«-7)  GO  TO  100C 
PREFER  =  ICO. 

PRESET  =  PREFFR  *  RAT (NUB) 

u5  FORMAT  (47X,30HSTRATEGIC  ANALYSIS  OF  EXCHANGE)  167700  S 

ilL  FORMAT  (/5iiX,3CHNEW  SIOE  CNF  ARSENAL  ) 

935  FORMAT (  //li  X  ,12  (4X,I  2,  4X )  > 

575  FORMAT <  35X,  47HSICE  CNE  COUNTER  DEFENSE  STRIKE  CHARACTERISTICS/ 

1  23X,  9JHSICE  ONF  ALLOCATIONS)  ON  DEFENSES  OF  RESOURCE  TYPE  (ROW 

2S)  DEFENDED  BY  OFDR.  TYPE  (COLS. >/29X ,1C H(ARRIV ALS ) /20 X, 


3 

90HSXCE  TWO  ( 

ORIG. 

NO. 

)  OF 

OEFOPS. 

OF 

RESOURCE 

TYPF 

(ROW 

4S) 

ATTACKED  BY  OFFS. 

/20X 

«5 

91HSI0E  TWO  ( 

SURV  . 

NO. 

>  OF 

OEFOPS. 

OF 

RESOURCE 

TYPE 

(ROW 

6S)  ATT ACKEO  EY  OFFS.  ) 

50  FORMAT  (/35X,57H AGGREGATE  COUNTERVALUE  OAMAGF  CHARACTERISTICS  OF  E  168350 
1XCHANGE//45X, 50 H TOT  AL  NUMBER  OF  ALL  CCA  TFO  (COLUMN)  CN  (ROW)  CITIES  1684C0 
2/19X,18HINT£RPPETATION  OF  , 17X , 44HNUM8FR  OF  ARRIVING  (COLUMN)  ON  (  168450 

3RON)  CITIES  /ll X ,20HE ACH  SET  OF  FNTPIES  , 15X , 38HNUMEER  OF  IMPACTIN  169500 
4G  (COLUMN)  ON  (ROW)  / 45 X, 58HNUMPER  OF  1  NT.  EQUIVALENTS  ON  (ROW)  C  168550 
5ITIES  FRCM  (COLUMN)  )  168600 

>55  FORMAT  (/20X , 8HM ISSILES,2UX ,7H BOMBERS, 2t X , 5HSLOMS ,20X , 11HPERCENTAG  168650 
IE  /75X, 13HSI0E  TWO  ONLY ,12X ,6H0AMAGE)  168700 

-  i9C  FORMAT  (45X,46HSIDE  ONE  COUNTERFORCE  STPIKE  CHARACTERISTICS  //30X  169050 

1, 56HMISSILES  CF  TYFE  (COLUMN)  ALLOCATED  AGAINST  TARGET  (RCW)/30X,5  169100 

22HMISSILE  CF  TYPE  (COLUMN)  ARRIVING  OVFR  TARGET  (ROM  /5X.19HINTER  169150 
3PRFTATI0N  CF  ,7X ,55HWAPHEAOS  OF  MISSILE  TYPE  (COLUMN)  ARRIVING  OVE  169210 
4R  TARGET  /5X,?(,HEACH  SET  OF  ENTRIES  ,5 X, 60HW ARHE ACS  OF  MISSILC  TYP  169259 
5E  (COLUMN)  IMPACTING  CN  TARGET  (ROW)  /30X ,66HHAPHEACS  OF  MISSILE  T  169300 
6YPE  (COLUMN)  IMPACTING  ON  EACH  MISSILE  (ROW)  /3C X .73HSURVI VAL  PROB  169350 
7.  OF  EACH  MISSILE  TYPE  (ROW)  FROM  ATTACK  PY  MISSILE  (COLUMN))  169400 

535  FORMAT  (//25X ,67HNUM9ER  OF  EACH  SIOE  TWO  MISSILE  TYFE  SURVIVING  CO  169850 
1UNTERFORCF  STRIKE  )  169900 

540  FORMAT  (4GX,46HSI0E  ONE  COUNTERVALUE  STRIKE  CHARACTERISTICS  //30X  169950 

1, 61HWEA°0NS  CF  TYPE  (COLUMN)  ALLOCATED  AGAINST  CITY  CLASS  (ROW)  /  173CU9 

26X,1BHINTERPRETATI0N  CF  , 6X ,56HHEAP0NS  OF  TYFE  (COLUMN)  ARRIVING  0  17J050 

3VER  CITY  CLASS  (ROW)  /6X,2CHEACH  SET  OF  ENTRIES  , 4X , 48HWE APONS  OFT  170100 
4YPE  (COLUMN)  IMPACTING  ON  CITY  CLASS  /3CX,74HN0.  OF  1  MT.  EOS.  FRO  170159 


5M  WEAPON  TYFF  (COLUMN)  IMPACTING  ON  CITY  CLASS  (ROW)  )  170209 

545  FORMAT  (41X,39HSI0E  TWO  SECOND  STRIKE  CHARACTERISTICS  /30X,68HNO.  173259 

10F  SIOE  TWO  WEAPONS  OF  TYPE  (COLUMN)  BEFORE  COUNTfRFORCE  STRIKE  /3  170300 

25X  ,71HNO,  OF  SIOE  TWO  WEAPONS  OF  TYPE  (COLUMN)  SURVIVING  COUNT  ERF  0  170359 

3RCE  STRIKE  /6X,1 8HINT  ERPRET  ATI  ON  OF  ,6X,62HNC.  CF  SICE  TWC  WEAPONS  179400 
4  OF  TYPE  (COLUMN)  ARRIVING  OVER  CITIES  /6X,20HEACH  SET  CF  ENTRIES  170450 
5,4X,7LHN0.  CF  SICE  TWO  WEAPONS  OF  TYPE  (COLUMN)  IMPACTING  ON  SIOE  170500 
60NE  CITIES  /33X » 70HNO .  OF  1  MT.  EOS •  FROM  EACH  WEAPCN  TYPE  (COLUMN  179550 
7)  IMPACTING  CN  CITIES  )  170609 

655  FORMAT  (1H1) 

66''  FORMAT  (29X,62H**********STRATEGIC  WEAPONS  EXCHANGE  MOOELS  (SWEM) 
1***********,//45X,34H  ALLOCATION  OPTIMIZATION  MOCEL  ) 

8n'j  FORMAT  //IX, 4HCIT Y,5 X, 12 (4X, A  6) /IX ,5HCL ASS) 

805  FORMAT  ».'1CX ,  12(  2X ,F6 .1 ,2X) ) 


910  FORMAT  (1CX,12I2X,F6.1,2X)> 

915  FORMAT  (2XtA6»2X»12(2X»F6tl»2X)) 

920  FORMAT  (1»'X,12(2X,F6.  3,  2X ) » 

925  FO<<MAT<  //11X,12 <4X,A6) > 

933  FORMAT  </li!X  ,  12<  2X  ,F6.1 ,2X)  > 

995  FORMAT  <4X,I2,4X ,12 (2X, F6 .1 ,2X ) ) 

iOi."  IFINU1.EQ.1)  CALL  STARTS 
RETURN 
ENO 
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FORTRAN  NOTATION 

wt 

DYBD(L,X,K) 

betbp(l) 

SBBX(l,j) 

#- 

dbbd(i,k) 

8din 

a*FtJ 

DF30C(NN1,L,KS) 

axj* 

NN1  depends  on  1  and  j 

a*F, , 

'  BF3CD(.m,l-,KS) 

axu 

NN1  depends  on  i  and  j 

a\ _ 

DBXX(NN2,KR) 

»■♦!,  KTTT 

NN2  depends  on  1  end  J 

a8s 

DBXD(NN2,KR) 

M|*it  3^1*1} 

s  NN2  depends  on  1  and  j 

a8s. 

DSX)C(NN2  T,) 

ax.A  ax.„ 

NN2  depends  on  m  and  k 

a8s. 

DSXD(NN2,L) 

w.A  aclM 

NN2  depends  on  m  and  k 

BFtl 

bfx(i,j,l) 

sf,, 

^7 

DFD(I,J,K) 
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*L 


b**r  dxU 


dfbxd(nni,m2,kn3) 

Ml  depends  on  l  and  m 
NN2  depends  on  i  and  j 
M3  depends  on  k  and  r 

DBBXX(M1,M2) 

Ml  depends  on  i  and  J 
M2  depends  on  k  and  r 


&!jL— 


.DBBXD(M1,NN2) 

Ml  depends  on  i  and  j 
M2  depends  on  k  and  r 


^1*  ^li 


DFDD  (NN1,L,KS) 

NKL  depends  on  i  and  j 


d33 

ddl  ddT, 

1+  u  1+  ij 


DBDD  (M2,KR) 

NIG  depends  on  i  and  j 


d*S. 

&d.i  K* 


DSDD  (M2,L) 

NN2  depends  on  m  and  k 


ddkP  dd14 


DFBDD  (Ml,  M2,  NN3) 
Ml  depends  on  l  and  m 
NIC  depends  on  i  and  j 
M3  depends  on  k  and  r 


daI 

3dkr  ddt7 


DBBDD  (Ml,  M2) 

NN1  depends  on  Jt  and  m 
NIC  depends  on  k  and  r 


D-fc 


n^n 


